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Theories for consideration divide themselves into two 
categories, those in which the primaries function as the 
measured entities, and those in which secondaries or higher 
order radiations play a significant or even a dominant role. 

Primaries as measured entities. 
such that the rays which 


There exists a curve 
relating altitude and latitude 
can just penetrate the earth’s magnetic field at an assigned 
latitude can penetrate the atmosphere down to the altitude 
given by the curve. The curve cuts the sea-level axis at a 
with that at which the 


latitude customarily associated 


increase of intensity with latitude ceases for further in- 


crease of latitude. For all altitudes above the curve in 
question, and in particular for all altitudes above sea level, 
with latitude less than the sea-level critical latitude, there 
can be no variation of primary intensity with altitude; 
and, appeal must be made to secondaries for variation of 
intensity with altitude in these regions. In regions where 
the primary intensity varies with altitude, that variation 
must depend upon the energy spectrum of the incoming 
rays. The type of energy distribution necessary to secure 
an exponential law is determined. Its form depends upon 
the law of energy absorption in the atmosphere. Two cases 


are considered. They correspond to 


—dF, dx=a 1); —dk,/dx=at hk 2) 


where £, is the energy at the distance x from the top of 
the homogeneous water equivalent atmosphere, a@ corre- 
sponds to loss by ionization, and \£, to loss by secondary 
production. Case (1) leads to an energy distribution which 
predicts at high altitudes a latitude variation which is 
enormous and out of all possible harmony with the facts. 
Case (2) leads to a reasonable high altitude latitude varia- 
tion. However, the other difficulties concerned with a 
predicted absence of variation of intensity with altitude 
in equatorial regions remain. 

Combination of primaries and secondaries as measured 
entities. The writer generalizes his former theory accord- 


* A paper delivered at the Tercentenary Conference of 
Arts and Sciences at Harvard University, September, 1936. 


ing to which the primaries produce secondaries in numbers 


which increase with increase of primary energy. The theory 


of the soft component is first considered, with (2) as a 


basis. The energy distribution function, F(/) necessary to 
give an exponential law for the primaries is of the form 
F(k) B, k+a/X)”: and 
IS us =(p 
the quantity W=(£+ea/\A), which is nearly the same as / 


with this form the absorption 


coethcient u 1)A. For mathematical convenience 


is used throughout for purposes of discussion instead of 


the true energy £. It is shown that if the number of 


secondaries m, accompanying a primary is given b 


n,=7W,*, then below the critical latitude u, = As, and above 
the critical latitude yu p—1)a. If, in 


Rk. A. Millikan’s observations we are to have yz, the same 


accordance with 


and below the critical latitude we must have 


p—1). If, 


meter of water equivalent we must have s=1.5 


abov 4 
moreover, u, has to have the value 0.5 per 
and so 
p=2.5 when X is assigned the value 0.33, which value is 
determined from the energy for penetration of the earth's 
field The 
equality of u, above and below the critical latitude leads 


magneti 610° ev, at the critical latitude 


inevitably to an expression for the secondary intensity 


whose evaluation involves the integration of dW/W be 
tween a lower limit which is finite and an upper limit 
which is infinite unless we assign a finite upper limit Fy, 
beyond which F(£) is zero. The concept of an upper limit 
Ey, and so an upper limit of Wy of W is consequently 
idmitted; and, for purposes of illustration the limit taken 
rhe 


* log 


intensity of the secondaries is 


if Vv W. below 
tude, where W, is the value of W for penetration of the 


earth's field at latitude ¢ 


is £ v= 5x10! ev. 


proportional to e* the critical lati 


It thus shows a latitude varia 


tion. There is also a stronger latitude variation of the 


primary intensity which, however, does not vary with 
altitude below the critical latitude. The combined sea-level 
latitude variation is too great to correspond to the obser- 
vations and it is necessary to look for the hard component 
to iron it out. 

The hard component. A hard component showing no 


latitude effect when combined with the foregoing contri- 
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butions to the intensity can be made to provide a satis- 
factory sea-level latitude variation, a satisfactory ratio of 
the intensities of the hard and soft components, and a 
satisfactory increase of latitude effect with altitude through 
the logarithmic term referred to above. A photon com- 
ponent can therefore fulfill the requirements; but, it is also 
shown that a charged particle component, operating mainly 
through its secondaries, will do equally well. Were it not 
for requirements of variation of latitude effect with alti- 
tude, it would be possible to tell the story of both the hard 
and the soft components as the outcome of the suitably 
chosen energy distribution in the incoming rays, which 
rays could then be all of one type. The increase of latitude 
effect with altitude demands, for the two types of pri- 
maries, properties which are different as regards production 
of secondaries. 

Miscellaneous phenomena. It appears that the average 
energy of the hard component primaries is determined 
mainly by the upper limit of their energies while that of 
the soft component primaries is determined mainly by the 
lower limit for penetration of the earth’s magnetic field. 
The conditions provided for the hard component second- 
aries are such that those secondaries completely outweigh 
the primaries in number; and, this explains why but few 
very high energy primaries are found in cloud chambers. 
The higher energy members of the soft component pri- 
maries are outweighed by their secondaries but this is not 
so for the lower energy members. The maximum in the 
intensity altitude curves shown by Regener’s results and 
by those of the Bartol Research Foundation of the Frank- 
lin Institute is attributed to a range phenomena of the 
secondaries and results from the limitations of the atmos- 


1. INTRODUCTION 

HE evidence drawn from the variation of 

cosmic-ray intensity, with latitude at high 
altitudes, has created a school strongly in favor 
of the view that the entities which enter our 
atmosphere from outer space and are directly or 
indirectly responsible for cosmic-ray phenomena 
are charged particles.! It is true that the facts 
are unable to deny that an appreciable portion 
of the cosmic rays observed at sea level may 
result from uncharged entities entering the at- 
mosphere, but they do not require this; and, 
according to A. H. Compton, practically a hun- 
dred percent of the rays observed at high alti- 
tudes owe their origin to charged particles. It is 
true that some of the high altitude data upon 
which Compton’s conclusions were based have 
become subject to question as the result of obser- 
vations secured since the time when he expressed 
the foregoing opinion. However, the likelihood 


1A. H. Compton, Rev. Sci. Inst. 7, 71 (1936). 
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¥phere. The so-called secondary hump in the said curves 


represents the transition which occurs on crossing the 


critical-altitude critical-latitude curve. The theory pro- 
vides for the fact that the east-west asymmetry increases 
with altitude more rapidly than does the total radiation. 
It provides for such apparently paradoxical conclusions as 
those resulting from the fact that there is very little 
asymmetry in shower production in lead, suggesting that 
showers are produced by hard rays, while shower produc 
tion increases with altitude more rapidly than the cosmic 
radiation, suggesting that it is produced by soft rays. The 
theory in fact makes shower production in lead increase 
with the energy of the primaries, but provides neverthe 
less for its rapid increase with altitude in such a manner 
as to give, at first sight, the impression that it is produced 
by a soft radiation. 

An alternative choice of constants. If 


value of \ as large as 0.5, it is possible to provide for 


one permits a 
us=0.5 by the assumption that the number of secondaries 
accompanying a soft component primary is proportional 
to W,, or approximately to E, as assumed in the writer's 
theory published earlier. This is attended with several 
advantages; but, necessitates relinquishment of the asso 
ciation between the sea-level critical latitude at 41 degrees 
and the critical latitude for the soft component. It becomes 
necessary to assume a second critical latitude phenomenon 
for the soft which critical latitude 
phenomenon does not show itself at all until an altitude 
becomes 


component primary, 


of 0.6 meter above sea level is reached. It now 


necessary to look to the hard component for the sea-level 
latitude effect. The soft component, however, provides the 
variation of latitude effect with altitude. 


of a purely corpuscular origin seems sufficiently 
great to make worth while an attempt to account 
for all of the phenomena in terms of charged 
particles. 

There are three possibilities as to the relation- 
ship between the primary rays and the phenom- 
ena observed. The first visions the observed rays 
as the primary particles themselves; the second 
regards them for the most part as long range 
secondaries, initiated directly or indirectly by 
the primaries, and perpetuating in large measure 
the directional characteristics of the primaries. 
The third possibility represents a combination of 
the first and second. Presumably nobody would 
deny that even if the first possibility represented 
the main story, there would be, in the observed 
phenomena, many secondaries, since such are 
actually observed in cloud chamber and other 
experiments. In the second or third type of hy- 
pothesis cited above, however, I do not imply a 
condition in which the role of the secondaries is 
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limited to a mere complicating function which 
tends to obscure the but rather one in 
which the law of their production, as a function 


story, 


of the energies of the primaries, etc., plays a 
definite part in determining the main features 
characteristic of cosmic-ray phenomena and, 

particular, the law of absorption. 
Neither is it implied that the secondaries are 
necessarily produced from the primaries in single 
Photons, may function as 


intermediaries. 


apparent 


acts. for example, 


2. NATURE OF PRIMARY Rays IN RELATION TO 
THE HypotHeEsis THAT THEY COoNsTI- 
TUTE THE OBSERVED ENTITIES 


_ Electrons (positive and negative), alpha-par- 
ticles, protons and other heavy particles have 
been suggested as primaries.' The investigations 
Bartol Foundation of the 
both in the laboratory? and 


have not given any evidence 


performed by the 
Franklin Institute 
in the stratosphere,* 
for the occurrence of protons, or other heavy 
particles functioning as the measured entities, 
however much such 
through 
may 


with observable absorption,‘ 
particles may participate in “‘silento”’ 
the action of the 
Thus, 


Maries as measured entities, 


secondaries which they 
in considering the case of pri- 
it would seem that 
we must confine our attention to positive and 


produce. 


negative electrons, although much of the analy- 
sis which follows is applicable to other entities 
as well. 

We shall confine attention throughout to the 
vertical intensity. The energy for vertical entry 
through the earth’s magnetic field is 17 K 10° ev 
at the equator, 6X 10° ev at 41 degrees magnetic 
latitude, and about 310° ev at 50° magnetic 
latitude for example. 


?W. F.G. 


Swann, Phys. Rev. 49, 478 (1936); also C. ¢ 
and D. D. Montgomery, W. E. Ramsey and W. F. ( 
Swann, Phys. Rev. 50, 403 (1936). 


*>L. H. Rumbaugh and G. L. Locher, Phys. Rev. 49, 
855 (1936). 

*The hypothesis that the primaries are the observed 
entities implies that a measurable fraction of them termi- 
nate their paths in the region of observation. The experi- 
ments cited in reference 2, were based upon an attempt to 
detect the protons near the ends of their ranges. High 
energy protons which passed right through the atmosphere 
without being stopped therein in measurable amount 


would be permissible in the light of reference 2. Such 


protons might be the basis of initiating secondary rz adiation 
which functioned as a large part of the observed radiation. 


ORY OF ¢ 
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OSMIC RAYS 


At the equator, the rays which can penetrate 
the earth’s magnetic field can also penetrate the 
lati 
energv of 


, however, a magnetic 
tude ¢ that the 
penetration of the field is only just sufficient to 
provide for penetration of the atmosphere. This 
latitude taken latitude at 
which intensity fails to increase 
further with increase of latitude. It is 41° and 
corresponds, as has been stated, to an energy of 
penetration of 6X 10° ev. 

For each altitude there is a critical magnetic 


atmosphere. There is 


sufficiently great 


is customarily as the 


the sea-level 


latitude ¢.., such that the rays which can pene- 
trate the magnetic field can also penetrate down 
to the altitude concerned. The form of the curve 
giving the relation between critical latitude and 
critical altitude is a function of the law governing 
the absorption of energy in the atmosphere. Two 
laws present themselves for our consideration in 


view of what will follow. They are 


—dkE, dx=a, 1 


~ 


—dE,/dx=\E,+a, 2 


where E£, is the energy of a primary ray at a 
depth x below the top of the homogeneous water 
equivalent atmosphere, a is the loss per unit of 
path in ionization, and } is a constant.° 

On the basis of the first law we cannot account 
for an energy of atmospheric penetration any- 
thing like as high as 6X 10° ev. The second form 
of law with \=0.333 and a=0.77 X 10° per meter 
of water equivalent® provides a satisfactory solu- 
tion, giving a range equal to that of the vertical 


atmosphere for an energy of 610° ev. The so- 
lution of (2) is 
E.t+ar=(E+a A)je, (3) 


where E is the energy of the ray at vertical entry 


has already been introduced 


48, 641 


5 The second of these laws 
by the writer in a former publication, Phys. Rev. 
(1935). See also W. F. G. Swann, Phys. Rev. 46, 828 (1934). 

° This value of a is based upon A. Eisl’s value 32.2 ev 
as the energy necessary to produce a pair of ions (see 
Ann. d. Physik 3, 277 (1929)), and upon an assumption 
of 38 ion pairs per centimeter of path at N.T.P., as found 
by W. E. Ramsey and the writer, in an investigation as 
yet unpublished. C. D. Anderson and S. H. Neddermeyer 
(Int. Conf. Phys. 1, 171 (1935)) quote data corresponding 
to values of @ ranging between 10° ev and 1.5 X10® ev 
but it would seem that their values require division by 
the density of air at N.T.P. in relation to the density of 
water X 10 
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Fic. 1. The curve gives for any magnetic latitude the 
depth to which primary electrons, which can just penetrate 
the earth’s magnetic field at that latitude, can reach in 
the homogeneous water equivalent atmosphere. The 
curve is calculated on the law of absorption specified in 
Eq. (2) with A=0.33 and a=0.77X108 per meter of 
water equivalent. 


into the atmosphere. The range R of the rays is 
given by putting E,=0, so that 


eR=1+4+ XE /a. (4) 


Fig. 1 shows the relation between the critical 
latitude and critical altitude as drawn from (4), 
with the foregoing quoted values of \ and a. 

Now it is manifest that if we regard the pri- 
mary rays as the measured entities there should 
be absolutely no variation of intensity with alti- 
tude in the regions lying above the curve in 
Fig. 1; and, in particular, there should be no 
variation with altitude for all latitudes below the 
sea-level critical latitude. Such a conclusion is 
contrary to the facts. In order therefore to secure 
an increase of intensity with altitude below the 
critical latitude it is necessary to suppose that the 
primaries are accompanied by secondaries whose 
number increases with the altitude and so with the 
energy of the primary.’ 

Quite apart from the foregoing consideration, 
however, other difficulties arise as a result of the 
law of distribution of incoming primary energies 
which it is necessary to assume in order to ac- 
count for the variation of intensity with altitude 
in the regions where such variation of intensity 
is to be expected. For the elucidation of this 
matter, and also to provide the material for sub- 
sequent developments we shall proceed to deter- 


71 exclude here considerations founded upon ‘‘strag- 
gling,”’ as such lines of escape do not seem to have been 
developed in sufficient detail to result in the correlation 
with all of the facts concerned. 
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mine the energy distribution necessary to provide 
for an exponential law or a combination of con- 


tributions obeying exponential laws. 


3. SPECTRAL ENERGY DISTRIBUTION NECESSARY 
TO PROVIDE EXPONENTIAL LAW OF 
ABSORPTION OF PRIMARIES 


In the remainder of this paper, whenever we 
speak of the ‘intensity of the rays,”’ we shall refer 
to the coefficient of dsdw in the expression for the 
number of rays passing per second through the 
horizontal element of area ds within the limits 
specified by the vertical solid angle dw. 

Let us suppose that the range R in the homo- 
geneous atmosphere is given as a function of its 
initial energy E by the relation 


Let N be the intensity at a depth x below the 
top of the homogeneous atmosphere. If dN is the 
change of N in a distance dx, then —dJN is the 
contribution to the intensity by those rays which, 
at the top of the atmosphere have energies be- 
tween E and E+dE, where 

E=y(x); dE=(d/dx)[ (x) jdx. 6 


Thus, if the energy distribution of the vertical 
rays at the top of the atmosphere is given by 
dN,=F(E)dE (7) 
at the poles, for example, where all of the rays 
can enter the atmosphere, we have 
—dN= F(y(x) ](dy/dx)dx. (8) 
If NV is to be exponential in x, with an absorp- 
tion coefficient u, then 
N=WN ve 9) 
and 


—dN=puNoe**dx. (10) 


It is to be observed that NV is given by 


Ne= [ F(E)dE. (11) 
0 


It is further to be observed that, in accordance 
with the conclusions reached in Section 2, rela- 
tion (9) holds only for depths x greater than 
X», Where 

¥(x,) = Ey. (12) 
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In other words, it holds for g>¢.z. By compari- 
son of (8) and (10), we thus have 


FL Y(x) ](dy/dx) =p Noe. 
The solution of this equation for F is 
F(&) =u Noe**® ((dx(€)) /dé), (13) 


where £ is any variable, and x(é) is the function 
which expresses R as a function of E through (5) 
in the form 

R= x(E). (14) 


Thus F(E) =uNoe*® (dR/dE), (15) 


where R is given as a function of E by (14), or 
if we prefer by (5). The energy spectrum thus 
becomes determined explicitly, directly we have 
assigned the range as a function of the energy 
through (5). It is to be observed, however, that 
the experimental existence of an exponential law 
of intensity up to any altitude 4 only determines 
the energy distribution for energies greater than 
a minimum value determined by E=y(h), so that 
at any latitude ¢, the energy distribution is only 
determined for values of E greater than the 
minimum energy of entry E,, as indeed might 
be concluded from fundamental considerations. 

Case of constant energy loss per unit of path. 
This corresponds to Eq. (1) and yields 


F(E) = (w/a) Noe*##'2, (16) 


or, starting with F(£) as fundamentally assigned 
we see that if 

F(E) = Ae** (17) 
an exponential law will result for the intensity 
for g>¢z with 


p=aB; No=(a/p)A=A/B. (18) 


With a=0.77X10°; and, w=0.5 for the soft 
component of the radiation, both expressed in 
terms of 1-meter water equivalent, we have 
8=0.65X10-*. The vertical intensity at the 
top of the atmosphere at magnetic latitude ¢ 
is given by integrating (17) from the entrance 
energy £, to infinity, and is proportional to exp 
(—0.65X10-* £). The ratio of the correspond- 
ing intensities at latitudes 41° and the equator, 
for example, is therefore exp [0.65 X 10~-8(17 —6) 
xX 10° ]=exp (71.5), an enormous quantity out of 
all harmony with any possibility in relation to 
the facts. 


Case where the energy loss per unit of path is 
composed of a constant part a, plus a part \E, 
proportional to the energy of the primary. Using 
(15) in combination with (4), which was derived 
from (2) we find 


pa” No a (itp A), 
F(E)= - | e+ ] (19) 
Ute A) X 


In terms of F(Z) regarded as fundamentally 
assigned, we conclude that provided that 


F(E) =B,/(E+a)?, (20) 


where B, and p are constants, an exponential 
law of the type (9) will result for the intensity, 
for ¢>¢z, with 


u=(p—1)Ar, (21) 
(p—1)No=(A/a)?"B,, (22) 


provided that a=a/X. 

It is of significance to note that both for this case 
and for that founded upon constant energy loss per 
unit of path, uw is not determined uniquely by the 
energy loss per centimeter of path, but depends also 
upon the law of energy distribution in the spectrum 
of the entering rays. In principle, the condition 
a=a/d is illogical since it makes F(Z) depend 
upon processes of energy loss in the atmosphere. 
The condition a=0 would remove the logical 
difficulty, but would destroy to some extent the 
exponential form for the intensity inherent in 
(20). However, since a/\ is of the order 4X 105, 
and so is much smaller than even the energy 
necessary for penetration of the earth’s magnetic 
field in all but the highest latitudes, the neglect 
of a/d is not of fundamental importance. In 
other words, the a/\ contribution is not of funda- 
mental importance one way or another. For these 
reasons we shall retain it, since its presence adds 
to the mathematical simplicity of our operations. 
In fact, in all that follows, we shall make our 
discussions in terms of E+a/X instead of E, and 
shall use, in fact, a new quantity W, defined by 


W=E+a/x. (23) 


Again, the vertical intensity at latitude ¢ is ob- 
tained by integrating F(£) from E,, to infinity, 
and by use of (20) and (21) is seen to be pro- 
portional to (E+a/\)“-”), where (p—1) =y/A. 
Thus, in terms of the values already quoted, we 
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find that the ratio of the intensities at 41° and 
at the equator is (173/63)'*=4.3. This value is 
quite reasonable in relation to the possibilities 
suggested by such experimental data as exist. 

However, as already pointed out, such varia- 
tion in the vertical intensity as is found at the 
top of the atmosphere should, on the present 
form of theory, be found also at all altitudes, 
including sea level, for ¢<¢». Moreover, no 
variation of any kind with altitude is provided 
for the regions ¢ <¢,, and no explanation is pro- 
vided for the so-called hard component. We shall 
attack these matters in the next section. 


4. DEVELOPMENT OF THE VIEW THAT THE OB- 
SERVED ENTITIES ARE COMBINATIONS OF 
PRIMARIES AND SECONDARIES 


Following the first analysis made by R. A. 
Millikan, various analyses of intensity-altitude 
cosmic-ray data have been made from time to 
time. W. Kramer® specified four components for 
the radiation with separate absorption coeffi- 
cients, and with corresponding intensity coeffi- 
cients J which correspond to the relative inten- 
sities for the different components when x=0. 
Kramer's values are as follows: 

iu 0.52 0.16 0.075 0.021 
J 92 4.3 3.3 0.4 


We shall omit the small component with u = 0.021. 
We shall assign »=0.5 for the soft component, 
and we shall lump the second and third compo- 
nents into one hard component with ~=0.12. 
We shall assign the intensities of the hard and 
soft components so as to result in the same ratio 
for the sea-level values as is given by Kramer's 
soft component, and the combined value of his 
second and third components. Thus we have 


BL 0.12 > 0.50 
a 9.6 ; 90.4 


In former publications, the writer has empha- 
sized the role played by secondaries and even 
rays from cosmic-ray bursts in contributing to 
the observed cosmic radiation ;° and in particular 





8 W. Kramer, Zeits. f. Physik 85, 411 (1933). 

*W. F. G. Swann, Phys. Rev. 43, 945 (1933); Phys. 
Rev. 44, 1025 (1933); Phys. Rev. 46, 432 (1934); also 
J. Frank. Inst. 220, 373 (1935), in which the property of 
high degree of conservation of direction by secondaries is 
demonstrated. 
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he has developed a theory” in which it was sup- 
posed that the primaries passed, for the most 
part, right through the atmosphere and pro- 
duced along their paths secondaries which per- 
petuated the directions of flight of the primaries. 
It was supposed, to a first approximation, that 
the product of the number of secondaries pro- 
duced, per unit of path, and the range of the 
secondaries was proportional to the energy of the 
primary. These assumptions led to an exponen- 
tial law for the measured intensity of the sec- 
ondaries. It was not supposed that the primaries 
of necessity produced their secondaries in single 
acts. Photons, for example, might intervene as 
intermediaries. In the more detailed development 
of the theory, the law of loss of energy was 
generalized to the type (2), and it was shown 
that, in this form the theory possessed the power 
to correlate a large number of experimental phe- 
nomena. The primaries paid to the measured 
radiation a direct contribution whose story was 
not developed in any great detail. In the present 
section we propose to amplify the story in greater 
detail, particularly with regard to the mutual 
relationship of the primaries and higher order 
radiations, and we shall develop the theory to 
the point of harmonizing its numerical predic- 
tions with experiment. 

The primary background. It will turn out that 
we shall have to use a primary background com- 
posed of contributions which, in the case of the 
background of the soft component, at any rate, 
obeys an exponential law. For this purpose, there- 
fore, we have the form of F(£) already developed 
in Eq. (20), which, as will readily be seen, is not 
limited to a single component but, in the form 


F(E) =B,/(E+a/\,)?+B,/(E+a/d,)%, (24) 


for example, is capable of giving rise to two pri- 
mary absorption coefficients of the form (p—1)A, 
and (g—1)A, as the result of a single type of 
primary entity. Unfortunately the desirable con- 
summation inherent in such a hope will have to 
be modified for reasons which will develop later ; 
but, in the meantime, it is of importance to 
realize that the different contributions to F(£) 
pay independent contributions to the primary 
intensity, and therefore to the intensities of the 





10 W. F. G. Swann, Phys. Rev. 46, 828 (1934); Phys. 
Rev. 48, 641 (1935). 
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secondaries arising from those primary con- 
tributions. 

In what follows, it will become necessary to 
introduce an upper limit Ey to the energy of the 
primaries. This will destroy the rigorous reali- 
zation of the exponential form for the component 
concerned, but it will turn out that no serious 
consequences result therefrom. 

We confine our attention to a single compo- 
nent. We utilize the quantity W defined in (23), 
which quantity is approximately the energy of 
a ray; and, bearing in mind (4), we obtain the 
vertical intensity N, of the primaries for ¢> ¢.,. 
by integrating F(£) up to Wy=Eyt+a/X, from 
a lower limit ae“, corresponding to E,=0. 
For ¢< ¢z, the lower limit of Wis W,=E,+a/X, 
where E, is the energy for vertical entry at the 
latitude ¢. Thus 


B, (*) ir AWau\ ~-? 
ms c lé ; | . a ( ) | 
(p—1)\a a 


for g>¢er, (25) 





B, 
N,= [We — Wy 7 
(p—1) 


for ¢<¢ez, (26) 


where ? is related to the absorption coefficient yu, 
by 
(p—1)A=u. (27) 


For latitudes below the sea-level critical latitude, 
there is no variation with altitude; and for 
greater latitudes NV, changes from the form speci- 
fied by (25) to that specified by (26) at the 
altitude given by Fig. 1. When Wy, is infinite, 
(25) is purely exponential. 

For the case where p is equal to or less than 
unity, an upper limit to Ey is necessary for the 
realization of a finite value for N,. The case 
where p=1 is of interest; and, in that case 


N,=B, log. [AWu/ae™] for ¢>¢2, (28) 
N,=B, log, [Wa W, | for ia) <= Pex: (29) 
When JN, does not obey a true exponential 
law, it is convenient to define a variable coeffi- 
cient of absorption yu, as 


us= —(1/N,)dN,/dx. 


We then have 


uz=(p—1)d/[1—(ae*/AW yy)? ‘i, (30) 
and, for p=1 
us=/log. [AWar/ae*]. (31) 


Both (30) and (31) approach their ideal values 
uw=(p—1)dA, when Wy=@. 


Theory of soft component 

Certain guiding considerations control us in 
our choice of a theory at this stage. The first 
concerns the striking fact, which has been care- 
fully established by R. A. Millikan and his col- 
laborators," that the coefficient of absorption for 
the soft component is 0.5 everywhere, i.e., for 
regions above and below the critical latitude. 

Suppose now that the number », of sec- 
ondaries accompanying a primary at the place 
where it has an energy £, is a function of E,, 
or rather of W, of the form n,=f(W,). Then, in 
view of (3) we have 


n,=f(We™), (32) 


where W is the value of W, at x=0. The num- 
ber of primaries lying within the range W, to 
W.+dW, is equal to the number which at x=0 
lie within the range W to W+dW, where W is 
related to W, by (3). Thus, at places for which 
¢<¢cz, and in particular for ¢ less than the sea- 
level critical latitude, we must have, for the 
intensity of the secondaries 


In order that N, shall be exponential in x, with 


Wm 
F(W)f(We-™*)dW. (33) 


We 
a coefficient y,, it is necessary and sufficient that 
f(We~*), that is m, shall be of the form 
n,=y(We-*)*, 

where sA=uy,, and y is a constant. 

In this case, with F(W) =B,/W? (33) becomes 

Ww 
N,=vB,e nef W°-»)dWw. (34) 
“We 


Incidentally, with the value A\=0.33 and 
us=0.5 already assigned, the condition s\=uz, 


gives s=1.5. 


uy, S. Bowen, R. A. Millikan and H. V. Neher, Phys. 
Rev. 44, 246 (1933); also R. A. Millikan, H. V. Neher 
and S. Korff, Phys. Rev. 49, 871 (1936). 
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For ¢> ¢cz, i.e., for latitudes greater than the 
critical latitude, the primaries also suffer absorp- 
tion. Now it has been shown elsewhere by the 
writer” that if the primaries obey a true expo- 
nential law, then over the range in which they 
obey that law there is a constancy of ‘‘quality”’ 
with altitude in the sense that, at all altitudes, 
the energy spectrum must be the same. All ener- 
gies change in the same proportion with altitude. 
Under these conditions, regardless of the law of 
production of secondaries as a function of the 
primary energy, the factor of proportionality 
between intensities of secondaries and of pri- 
maries must be the same at all altitudes. Under 
such conditions the secondary intensity must obey 
an exponential law with the same absorption coeffi- 
cient as that of the primaries. Thus if Ey is 
sufficiently large to guarantee a true exponential 
relation for N,, in order to obtain a coefficient 
s\ for the secondaries we must have, in accord- 
ance with (27) for Ey infinite, 


(p—1)A=syX. 


Hence (p—1)=s. (35) 
Thus (34) yields 


N,=7B,e-*** log. (Wu /W,) for ¢e<¢.z. (36) 


For ¢>¢-z, the only change necessitated is the 
replacement of the lower limit W, by ae/X. 


Thus 
N,=7B ye" log. (AW /ae®”) for ¢>¢ez. (37) 


It is to be noted that regardless of the values 
of s and of \ the logarithmic form in (36) and 
(37) is inevitable. 

With s=1.5, necessitated by the values of u, 
and \ concerned, (35) gives p=2.5, so that 


F(£) = B,/W?-. (38) 
For the general case 
p=1+n/d. (39) 


It will readily be seen that with the value 
p=2.5, and with an upper limit as low even as 
5X10" ev the departure of (25) from the true 
exponential form is negligible. 

It is of interest to observe that the equality of 
us for regions above and below the critical latitude 
leads inevitably to the form of dependence of sec- 





2 W. F. G. Swann, Phys. Rev. 47, 575 (1935). 
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ondary production upon primary energy and (o the 
value of p for the primaries. 

Note on the significance of the relation between 
secondary production and primary energy. We 
have regarded the number , of secondaries 
accompanying a primary as a function of W,, i.e., 


n,=f(W.). (40) 


If r is the range of a secondary, and o is the 
number of secondaries produced per unit of path, 
we have, approximately 


n,=or=f(W.). (41) 


If « is the energy of a secondary, and if we 
regard the term AZ, in (2) as representative of 
the loss encountered per unit of path in produc- 
ing these secondaries, then 


oe =E,=)W,, approximately. (42) 
Thus, from (41) and (42) 
e/r=W./ (f(W.)). (43) 


For the special case to which we have been led, 
where n,=7W,’, we have 


e/r=i/yW,°-”. (44) 


The content of our assumption is comprised in 
the following : Suppose that theory or experiment 
has given us some relation between 7 and e. Then 
(43) becomes an expression of an assumption as 
to how e depends upon £,. Next (42) becomes 
a statement of how a is to depend upon £, in 
order that the energy loss per unit of path shall 
be proportional to E,. With (42) and (43) hold- 
ing, (40) becomes provided for. In subsequent 
determination of the constants concerned, the 
constant y, for example, our main care must be 
to see that e/r as given by (44) does not come 
out smaller than would be reasonable for the 
average energy per unit range of a secondary. 
Thus, ¢/r should not be less than 10° ev per 
meter of water equivalent, so that if the upper 
limit of the energies concerned is Ey (or Wy 
approximately), we must have 


<A (Wx,° > 10%), (45) 
for the case where s=1.5. 


The hard component 
At this stage our theory has given us a satis- 
factory preliminary account of the soft compo- 
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Fic. 2. Data, quoted from A. H. Compton, on the 
variation of cosmic-ray intensity with latitude at sea level 
and 4360 meters altitude. 


nent. We must, however, examine it in relation 
to the latitude effect, and here the hard compo- 
nent will become involved. We shall denote by 
Ry the ratio of the vertical sea-level intensity at 
the sea-level critical latitude to the corresponding 
sea-level intensity at the equator. R, shall be the 
corresponding ratio for the same latitudes at the 
depth x below the top of the atmosphere. 

For the secondaries, R is obtainable from (36). 
For an upper energy limit equal to 5X10" ev 
R,=Ro=1.92, and for an upper limit of 10" ev 
it is 1.56. 

The latitude ratio for the primaries is given 
by (26) ; and, for the above upper limits the ratios 
are approximately the same and are 5.5 and 4.9, 
respectively, since ~=2.5. The resultant sea- 
level latitude ratio for the soft component cannot 
be less than the value given by the logarithmic 
term. It is then necessary to look to the hard 
component to iron out the sea-level ratio to the 
observed value, say 1.14," leaving the soft com- 
ponent dominant at high altitudes, however. 

An easy way to dispose of the hard component 
is to attribute it to photons, in which case it 
would have no latitude effect. On a preliminary 
survey such a procedure presents difficulties 
which tend to evaporate upon further scrutiny, 
however. The situation is worthy of discussion. 

The latitude effect in the light of a hard compo- 
nent which itself shows no latitude effect. If we 
take Kramer’s data, we find that the intensities 
of the hard and soft components at sea level are 
in the ratio 4.8 to unity. Taking these figures as 
applying to some place above the sea-level criti- 


cal latitude, we find that if y is the ratio of the 


18 P. Auger and L. Leprince-Ringuet, Int. Cong. Phys. 
1, 188 (1935). 


vertical intensity for the soft component at 
¢=¢- to that at ¢=0, and if 1.14 is the corre- 
sponding ratio for the total intensity, and if the 
hard component has no latitude effect 


Rey = (N,+4.8N;,)/(Ns/y+4.8N,) = 1.14, 


where JN, is the intensity of the soft component. 
Thus y=3.5. Incidentally this represents the 
maximum latitude effect which can be produced 
between ¢.o and ¢=0 at any altitude apart from 
the invocation of circumstances of a type extra- 
neous to those involved in the expression of the 
intensity in the form of the sum of two expo- 
nential contributions. Of course, the latitude 
effect extends to greater latitudes at higher alti- 
tudes. In fact, at any altitude, it extends to the 
latitudes given by Fig. 1, and naturally the 
intensity ratio goes on increasing up to this 
latitude. The value y=3.5 while theoretically 
attained only at the top of the atmosphere would 
be approached at altitudes for which the soft 
component dominates the hard component. Thus, 
for example, at an altitude corresponding to x =6 
meters of water, we should have for the ratio of 
the intensity at ¢=42° to that at ¢=0, a value 
Rx given by 


Re = (4.52 +4.67) /(4.58/3.5+4.8) =1.55, 


for the case where R,,=1.14 for sea level. The 
numbers 4.58 and 4.8 represent the relative ver- 
tical intensities at x=6 meters of water, calcu- 
lated from Kramer's data cited earlier, in which 
the relative intensities of the soft and hard com- 
ponents at x=0 are 90.4 and 9.6, respectively, 
and in which the corresponding absorption co- 
efficients are 0.5 and 0.12. 

Now at an altitude of 4360 meters, which is 
not far from the equivalent of x=6 meters of 
water, A. H. Compton’s data shown in Fig. 2 
give a value of Rs equal to 1.30. The excess of 
this ratio over unity, which is the thing which 
counts, is thus only a little more than one-half 
of the excess calculated on the assumption that 
the hard component shows no latitude effect. A 
latitude effect for the hard component would 
reduce the foregoing discrepancies, but would 
call for a charged corpuscular origin for the hard 
component, or at any rate for a portion of it. 

While the foregoing considerations appear to 
mitigate against the claim of photons for the 
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position of the hard component, their force be- 
comes annulled in the light of the conditions to 
which our theory of the soft component has led 

conditions in which for ¢< ¢.» there is a primary 
contribution to the intensity, a contribution in- 
dependent of altitude but dependent on latitude, 
and which possesses the power to provide at sea 
level the necessary latitude variation which a 
photon component would lack. In the light of 
these considerations, let us therefore examine the 
possibilities inherent in the assumption of a ‘‘ non- 
field sensitive,’’ a noncharged hard component, 
remembering of course that if such a component 
is to affect Geiger counters, it must act through 
its secondaries. Using (25), (26), (36) and (37) 
and recalling that p= 2.5 we find for the vertical 


cosmic-ray intensity /, 


2/r\ 15 Wr 
1=B (2) +78. (>) 
3\a ae** 


2 1 \'s 
+ A e-Har — B,( ) for ¢> ez, (46) 
3 Ww 
2 1 1 
|= a r 
3 We Wy! 
Wa 
+B] log ( )} usrt. 4 e-uaz 
w,/. 
for ¢< Pez, (47) 


where A is the intensity of the hard component 
whose absorption coefficient is pa. 

Only for g > ¢.z does our theory give a primary 
as well as a secondary variation with altitude for 
the soft component ; and, it must be in this region 
that we regard Kramer's analysis as applicable. 
The contribution in (46) which is proportional to 
y is not truly exponential, but involves a slowly 
varying logarithmic coefficient. For the purpose 
of comparison with Kramer's analysis we shall 
regard the comparison as made at sea level, and 


14 No apology is needed for our departing from the true 
exponential forms. As a matter of fact the correct forms 
for the intensities above and below the critical latitude 
must fit together on the curve represented in Fig. 1. The 
sum of two exponentials, with constant absorption coeffi- 
cients could not possibly fit together on this curve even 
with the best adjustment of the intensity coefficients. 
Thus the form of theory leading to true exponentials 
would be illogical. The departure from the true exponential 
forms in our theory provide for just those elements neces- 
sary to secure the fit along the curve Fig. 1. 
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shall evaluate the logarithmic coefficient in 
question there. 

We assume A=0.33 a=0.77X10%, as 
before, write youyX10", Ay=AX10", and 
assume, for purposes of trial that Ey =5 X10" ev, 
which corresponds to Wy=502X10° ev. Then 
Kramer’s value 9.4 for the ratio of the intensities 
of the soft and hard components causes (46) 


to yield 


and 


0.190B,+2.041y.B,=9.4A 0- (48) 


Eq. (47), with u,=0.5 and uw, =0.12 gives, for the 
sea-level critical-latitude ratio 


1.27 X10-*B,,+13.76 X10-*yoB,+0.301A 5 


_ 0» 


2.88 104B,+7.12X10-y70B,+0.301A49 (49) 
(48) and (49) yield 
(7.35 +79.1y0) /(6.37+72.4y0)=Ro. (50) 


If we assign Ry=1.14, we find yo=0.043 so that 
y = 0.043 X 107", 

This value of y leads to results which are 
not unreasonable. Thus, for the highest en- 
ergy assigned, Wy,=502X10° ev (44) gives 
e/r=0.3X10*. This value, for the average energy 
per unit range of the secondary is dangerously 
small. However, in all places where observations 
are made which would have any significance in 
relation to our formulae, W would be much 
smaller than 5X 10'° ev. Thus at x =4 meters, the 
maximum energy would be 5 X10!°Xexp [ —1.3] 
= 1.210", and the corresponding value of ¢/r 
would be about 0.75 X 10° ev. This represents the 
worst case, and for all lower values of W, the 
value of ¢€/r is greater. 

The foregoing value of y means that a 10'° ev 
primary would be accompanied by about 40 
secondaries. 

The ratio of the total number of secondaries to 
primaries for ¢>¢,z is given by the ratio of the 
two contributions to the soft component in (46). 
It does not vary greatly with altitude (for 
¢>¢-z) and amounts at sea level and at mag- 
netic latitude 50° to 1 to 1.9. However, as already 
implied, the ratio is greater for the high energy 
primaries, being, as we have already shown, 40 for 
a 10'° ev primary. It is thus easy to see why cloud 
chamber photographs would not show many high 
energy particles. Such particles would be repre- 
sented mainly by their low energy secondaries. 
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Fic. 3. Relative cosmic-ray intensities as a function of 
depth below the top of the homogeneous water equivalent 
atmosphere. The curve is drawn from the theoretical 
formulae given in Eqs. (46) and (47). The point P is 
representative of the secondary hump in the intensity 
altitude curves. 


Variation of latitude effect with altitude. On ac- 
count of the logarithmic term in (47), the latitude 
effect varies with altitude. With the value of y 
chosen to give a sea-level latitude ratio of 1.14, 
the value of R, at 6 meters of water equivalent 
depth comes out as 1.19; and, at x =4 meters it is 
1.29. The maximum vertical current ratio for any 
altitude is the ratio of the intensity at the latitude 
and altitude determined by Fig. 1, to the inten- 
sity at the same altitude at the equator. The 
ratio in question comes out as 1.79 for x=6 
meters and 2.66 for x =4 meters. 

The increase of latitude effect with altitude is 
thus of the same order of magnitude as to cor- 
respond to Compton's data but is somewhat 
smaller. If the hard component possessed a sea- 
level latitude effect, the effect would be to weaken 
still further the increase of latitude effect with 
altitude. The increase of latitude effect with 
altitude is bound up with the logarithmic term 
in (47) ; and, as we have seen the logarithmic form 
is inevitable if u, is to have the same value above 
and below the critical latitude. The maximum 
latitude effect attainable at high altitudes is 
limited by the logarithmic form of this term. The 
increase of latitude ratio with altitude is very 
sensitive to the choice of the upper limit. Only 
when that upper limit approaches close to the 
limit necessary for vertical entry at the equator 
does one obtain any increase of latitude ratio with 
altitude at all. In fact, our theory makes the 
increase with altitude a symbol of the existence 


of an upper limit. If Ey =1.7 X10" ev, which is 
the energy for vertical entry at the equator, the 
latitude effect of the logarithmic term will be 
infinite. We are limited, however, in the degree of 
smallness of the value which we may assign to 
Ew, since we must permit rays to’come in from 
the due east and due west ; and, to secure entry in 
the equatorial plane this would necessitate an 
energy of about 6.010" ev. Even if we should 
adopt as low a value as 2.5 X10" ev for Ey, the 
maximum latitude ratio for the logarithmic term 
for x=4 meters would only amount to 8.75, and 
the complete latitude ratio, including the hard 
component would be 3.0. All of the published 
data on variations of intensity with latitude give 
ratios less than 2, when the observations have 
been made with a single instrument at different 
latitudes. Such latitude ratios of the order of 2 to 
1, and even 10 to 1 and more as have been quoted 
as applicable to the top of the atmosphere, depend 
upon comparison of observations taken with dif- 
ferent instruments in different flights. Such com- 
parisons are attended by the greatest difficulty. 
In the first place, the whole significance of the in- 
tensity is modified by the existence of a max- 
imum. This maximum, on the theory developed 
in this paper, would simply be indicative of the 
fact that near the top of the atmosphere the 
range of the secondary electrons is greater than 
that corresponding to the whole of the air above 
the point of observation, so that the intensity 
suffers in relation to the value given by the 
theory on account of the absence of enough air to 
produce all of the secondaries required by theory. 
On this basis the observations of the Bartol 
Foundation lead to a secondary range of about 82 
centimeters of water equivalent. However, our 
theory provides for a variation of the energy and 
so of the range of the secondaries with primary 
energy. There should, therefore, on this account be 
some variations in the altitude of the maximum 
with latitude. Then, matters are complicated 
considerably by the fact that, as the Bartol Re- 
search Foundation observations have shown," 


15 See, for example, W. F. G. Swann and G. L. Locher, 
National Geographic Society, Contributed Technical 
Papers, Stratosphere Series, No. 1, 7 (1935), in which 
publication it appears that the horizontal intensity forms 
20 percent of the vertical at 40,000 feet ; see also W. F. G. 
Swann, J. Frank. Inst. 222, 1-11 (1936) in which from 
our observations in the flight by Professor and Mrs. Jean 
Piccard we obtained more than 50 percent for the ratio 
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the horizontal intensity is comparable with the 
vertical intensity in the stratosphere on account 
of the bending of the paths of the secondary rays 
by the earth’s magnetic field. Undoubtedly many 
of the secondary rays actually approach the ap- 
paratus from a direction below horizontal, so that 
the variation of shielding conditions in different 
experiments is apt to have an important bearing 
upon the results, particularly as it was formerly 
supposed that shields were only effective when 
placed above the ionization chambers. 

The secondary hump in the intensity-altitude 
curve. Intimately bound up with the latitude 
effect is another having to do, possibly, with the 
secondary hump—the hump below the maximum 
—characteristic of intensity-altitude curves. At 
any latitude, there is a value of x determined by 
Fig. 1, for which the expression for the vertical 
intensity changes from the type characterized by 
(46), to that characterized by (47). Above this 
altitude the chief change is that resulting from 
the absence of the direct contribution to the soft 
component by its primaries. We should thus ex- 
pect a hump in the ionization altitude curve at 
the value of x concerned. This hump should start 
at sea level for ¢= ¢.. =41°, and should gradually 
climb up the curve with increasing latitude. In- 
tensity-altitude curves obtained with ionization 
chambers must, of necessity, tend to obscure this 
hump, since there is a different critical latitude 
corresponding to each direction from the zenith. 
Fig. 3 represents the vertical intensity-altitude 
curve plotted from Eqs. (46) and (47), with the 
constants already determined and for magnetic 
latitude 55°. The curve does not, of course, show 
a maximum, since the maximum represents an 
added feature depending upon the limitations of 
the atmosphere as already stated. It does show 
the secondary hump at x=6.6 meters of water. 
This value of x is considerably greater than the 
value suggested by high altitude data. It is for 
this reason moreover that, in spite of our adjust- 
ment of the absorption coefficients and the rela- 
tive sea-level intensities of the hard and soft 
components to their experimentally established 


of the horizontal to the vertical intensity. In the last 
flight made by Captain A. W. Stevens and O. Anderson 
our apparatus gave practically as much horizontal as 
vertical intensity. It is of course to be observed that 
horizontal intensities as measured by Geiger counters 
correspond to the sum of the intensities for two diametri- 
cally opposite directions. 
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values, the ratio of the total intensity even at 
x =0 to that at sea level is given by the curve as 
only about 17 instead of three or four times that 
value, as it should be. Above the secondary 
hump, the radiation becomes robbed not of its 
absorption coefficients but of part of its intensity 
factor; and, if the hump comes at too low an 
altitude, the ratios of the high altitude intensities 
to those at sea level became seriously affected. 
While the exact position of the secondary hump 
may be regarded at first sight as a matter of 
refinement which an approximate theory cannot 
be expected to reveal, it must be admitted that 
the defect inherent in a failure to make the right 
prediction becomes more vital on mature con- 
sideration. Fortunately, a way of overcoming the 
difficulty exists; but, since it involves a rather 
radical change in our views as to the significance 
of the critical latitude, and since it has been 
developed since this paper was presented at the 
Harvard Tercentenary Conference, it is here 
relegated to an Addendum at the end of the paper. 

Possibilities inherent in a field sensitive hard 
component. If the hard component is of a charged 
particle type, and is to show an increase with 
altitude in equatorial regions, it also must be 
represented by secondaries. The coefficient of 
absorption is one-quarter that of the soft com- 
ponent, so that if the number of secondaries 
accompanying a hard component primary, which 
we shall call a g primary is to be proportional to 
(E+a/X,)*9, i.e., We, we must again have 
Mh=S qd. 

The most elegant solution of the problem 
would be one in which we had simply one form of 
primary, with an energy spectrum given by 


F(W) =B,/W°+B,/W%, (51) 


and with a law of production of secondaries such 
that the number ”, of secondaries accompanying 
a primary was 


n,=yW'>+kRW'?", (52) 


The working out of the theory would then tell the 
story of both the hard component and the soft 
component simply as the outcome of the form of 
the energy distribution in the incoming primaries. 
The constants \ and a@ would, of course, be the 
same for the g primaries as for the p primaries. As 
a matter of fact, it is perfectly possible to fit up a 
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theory along these lines, and provide for the 
right ratio of the intensities of the components 
and for the sea-level latitude ratio. However, 
difficulty is encountered in the matter of the 
variation of latitude effect with altitude. The 
significant elements are these: If for the g 
primaries we proceed as we have done for 
the p primaries, we should conclude that 
(g—1) =u, \=0.38; and, in line with our former 
demonstration, the contribution of the g second- 
aries to the intensity through the term kW'*/4 
would come out proportional to B, log, (Wa/ W,) 
below the critical latitude. However, there would 
be two additional “‘cross product terms” in the 
contribution to the secondaries, a term propor- 
tional to the integral of (1/ W?*)(W'!*4)kB,, and 
a term proportional to (1/W'*5)(W'+)yB,. Below 
the critical latitude, for example, these two terms 
are proportional, respectively, to B,k(W,-':*8 
—Wy**) and Byy(Wx "—W,°-”). Of these 
two terms, the former contributes to the hard 
component and the latter to the soft com- 
ponent. There are thus two contributions to 
the soft component secondaries, one propor- 
tional to B,y log. Wy W,, and one proportional 
to Boy(Wa ” — W,2-”). Now when we make use 
of the intensity ratio of the soft and hard com- 
ponents, it becomes possible to show that the 
latter contribution to the secondaries dominates 
the former, and so determines in large measure 
the latitude effect at high altitudes since it is part 
of the soft component. However, the latitude 
effect of this term is weaker than that of the 
logarithmic term. It is, in fact weaker in this 
respect than all of the other contributions from 
the p primaries or secondaries or the g primaries 
or secondaries. Thus, its dominance at high alti- 
tudes tends to cause a decrease rather than an 
increase of latitude effect with altitude. Jt is for 
this reason that it appears necessary to separate the 
pb primaries and the q primaries into two distinct 
classes, one whose law of production is given by 
n,=yW,'*, and the other whose law of produc- 
tion is given by kW,'*"*. As already stated, apart 
from the requirements of the increase of latitude 
effect with altitude, it would be possible to fit up 
a theory based upon (52), although even then it 
would be necessary, as indeed it is in any case, to 
assume for the g primaries a value of g less than 
1.38, or even a negative value of g so as to avoid 


having too large a sea-level latitude effect for the 
hard component. 

Let us then take for the g primaries a contribu- 
tion to F(W) of the form 


F(W)=B,W" (53) 


and assume that the number of secondaries » 
accompanying a primary is 


n,=kW, 
where Sq\q= Un. 


We then have for the primary intensity 


B, aeary ntl 
N,= om | wess—( ) 
(n+1) do 


for g>ecz, (54) 


N, L(Wa,)""*—(We)"*] 


(n+1) 
for g< Sem (9) 
and for the secondary intensity 


kB, 
N.,.= | (wr, hg 
(n+1+4r/dA,) 


— (aee7/Q,) "ttn se “e® for g>Pecz, (56) 


kB, 
Nu= — | (Wynter 
(n+1+pr/dq) 


wn 
~I 
— 


—(We,)"titn sf “et for o< Goce 


Here we have recognized values of \Wy and 
critical latitudes which may be different for the q¢ 
primaries than for the p primaries. 

The essential thing to observe is that the 
greater , the more does NV, approach a radiation 
with zero absorption coefficient, the more nearly 
does N,, approach a true exponential with ab- 
sorption coefficient u,, and the smaller the lati- 
tude effects of both NV, and N,,, as shown by (55) 
and (57). Thus, for example, if m=2, and 
Wu ,=5X10" ev, the ratio of N,, at g=41 
degrees to the value at the equator is less than 
[1—(6/50)*]/[1—(17/50)*], ie., less than 4 per- 
cent. For »=3, it is only about 1.3 percent. 
Moreover, there is no limitation upon k, and so 








~ 
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upon the ratio of the secondaries to the primaries 
other than that imposed by the necessity of (44) 
leading to a reasonable value for e/r where y is 
replaced by k, \ by A,, and s by s,. Remembering 
that ”,, the number of secondaries per primary 
is equal to kW*, (44) becomes ¢«/r=\d,W/n,. 
Hence, if «7 is to be greater than 10° and if A, 
is comparable with \ for the primaries, we have 
W/n,>3X10*. Thus for a 5X10" ev g primary, 
we can have ”, comparable with 200. This would 
correspond to k=200/(5X10!")°**; and a 10° 
volt g primary would be accompanied by 
200/(50)°**=50 secondaries. Thus, under such 
conditions, the secondaries would dominate the 
primaries. 

With a very small latitude effect provided for 
the hard component, there would be no marked 
critical latitude for it ; and, its relationship to the 
soft component would be exactly the same as that 
already discussed in the case of the nonfield 
sensitive hard component. The net latitude effect 
would be calculated as before, and the variation 
of net latitude effect with altitude would follow 
as before. 

Thus, it appears possible to supplement the 
story of the soft component equally well with a 
hard component which is field sensitive or with 
one which is not—a photon component, for ex- 
ample. The ratio of intensities of hard component 
to soft component can be provided for. The cor- 
rect latitude effect can be provided for as can also 
the variation of latitude effect with altitude. 


Relation of theory to shower production and to 
asymmetry 


Shower production. It is a well established ex- 
perimental fact that the rate of production of 
small showers from lead increases with altitude 
at a rate greater than the rate of increase of the 
total radiation. In fact, the shower production 
increases approximately with a coefficient equal 
to that of the soft component of the radiation. 
This has given rise to the statement that the 
shower producing radiation is a soft radiation. 
According to our theory, if the showers in lead 
are produced according to a law anything like 
that for the air secondaries, it is the most energetic 
primary rays which are the most efficient in the 
production of showers; but, in spite of this, they 
endow the shower production with a coefficient 
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of increase with altitude which is large in propor- 
tion to the power of the energy which determines 
the shower production. It is true that, in the case 
of a photon hard component we could not predict 
the net increase of shower production with alti- 
tude without knowledge of the relative efficiency 
of the photons and the charged particles in pro- 
ducing showers (directly or through photon 
intermediaries) ; but, in the case of a hard com- 
ponent composed of particles which lose energy 
more or less continuously along their paths, so 
long as the shower production in lead by such 
particles depends upon a higher power of the 
energy than does the air secondary production by 
the hard component, shower production will 
increase with altitude more rapidly than does the 
cosmic-ray intensity. In the foregoing connec- 
tions, it may be well to emphasize again, that the 
“softness” of the component, in the case of the p 
primaries, at any rate depends not merely upon 
the energy absorption but also upon the energy 
distribution of the incoming radiation. 

East-west asymmetry. The view here taken re- 
moves a paradox which has existed between the 
variation of shower production with altitude, and 
the east-west asymmetry in shower production. 
According to the naive view that showers are 
produced by a soft radiation we might expect that 
they would show great asymmetry, since the 
east-west asymmetry is greater for the primary 
rays of smaller energy than for those of larger 
energy. Experimentally, however, very little 
asymmetry has been found in shower production. 
The reason on our present view, is very clear. It 
is the rays of high energy which produce the 
showers, so the showers show very little asym- 
metry. On the other hand, for the reasons stated, 
the intensity of shower production increases with 
altitude at a rate determined by a large co- 
efficient. 

Dependence of asymmetry upon altitude. The 
story of the dependence of asymmetry upon 
altitude is of interest. To fix our ideas let us 
imagine two similar cosmic-ray telescopes in a 
plane parallel to the equatorial plane and 
pointed at equal angles on the east and on the 
west side of the zenith. The asymmetry observed 
by these two instruments depends upon the rays 
which enter the atmosphere within certain energy 
limits, and which also’ succeed in penetrating 
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down to the instruments. If the apparatus be now 
taken to a higher altitude, the asymmetry will 
alter mainly through the effect of decreased 
absorption permitting rays at the lower energy 
limit to enter, which rays would not have pene- 
trated to the apparatus at the lower depths. 
Now, according to the form of theory we have 
developed the average entrance energy of the p 
primaries is approximately 3£, and is deter- 
mined by the Jower limit of the energy for entry, 
while even for the case of a charged particle hard 
component the average entrance energy of the g 
primaries is approximately (n+1)Ew,/(m+2) 
and is determined by the upper limit of the 
energy. Thus, it is the » primaries which are 
mainly concerned in the variations of asymmetry 
with altitude, which variations are determined 
by variations in the numbers of rays entering the 
instruments at the lower energy limit. We should 
thus expect that variation of asymmetry with 
altitude would be determined by the soft component 
absorption coefficient, rather than by the varia- 
tions with altitude of the cosmic radiation as a 
whole. A further conclusion derivable from the 
theory, and one in harmony with the results of T. 
H. Johnson, is the conclusion to the effect that the 


latitude effect for shower production is smaller 
than for the total radiation. Since, according to 
the theory, shower production depends upon the 
higher energy rays, its latitude effect will be 
determined largely by the latitude effect for 
those rays ; but, it is only the rays lying below the 
upper limit for equatorial entry that show any 
latitude effect at all. 


The nature of the hard component primaries and 
the soft component primaries 


While the present theory does not lead to any 
very specific conclusions as to the actual nature 
of the particles concerned, the fact that the east- 
west asymmetry is determined by the low energy 
rays, points to its becoming associated with the 
p primaries ; and, the fact that asymmetry meas- 
urements require positively charged rays suggest 
that the p primaries may be protons, while the g 
primaries are electrons. An assumption of protons 
for the » primaries would not be inconsistent 
with our failure to find such protons at sea level,” 
since the p primaries and their secondaries are 
associated with the soft component, and so play 
but little role at sea level. 


ADDENDUM 


(Made since delivery of the paper at the Harvard Tercentenary) 


The chief weakness in the theory as given in 
the foregoing lies in its prediction of too low an 
altitude for the ‘‘secondary hump” in the 
intensity-altitude curve. This difficulty could be 
alleviated if we could permit for the ~ primaries a 
larger value of X. If, in fact we could take \=0.5 
instead of 0.32, the modified curve corresponding 
to Fig. 1 would show an altitude corresponding to 
x=5.3 meters for latitude 55°, so that the second- 
ary hump would come at x =5.3 for this latitude, 
and would be better in agreement with such 
experimental data as exist. A still larger value of 
\ would improve matters still further. The value 
4=0.5, would result in s=1, so that the number 
of secondaries accompanying a primary would be 
proportional to the first power of the energy as 
assumed in the first form of the theory which I 
published in 1934.!° With s=1, Eq. (35) requires 
p=2, so that F(W,) would be inversely propor- 
tional to W,?. With s=1, moreover, requirement 


(44) now reverts to e/r=X/y so that y can be as 
large as \X10~-° viz. 0.5 10-5 ev without caus- 
ing €/r to be less than 10° ev per meter of water. 
With this value of y, a 10"° ev electron would be 
accompanied by 50 secondaries and there is no 
value of W,, no matter how large, which would 
give trouble as regards condition (44). However, 
with A=0.5, and a=0.77 X 10°, Eq. (4) gives for a 
1.7 X10" ev electron a range of only 9.4 meters 
of water equivalent. Thus, there would be no 
sea-level critical latitude for the soft component. 
It would be necessary to leave the sea-level 
latitude variation entirely to the hard com- 
ponent. Eqs. (55) and (57) provide the story. For 
n=0.5 in (53), and with u,/A,=0.36, which cor- 
responds to 4, =0.12, and the value \,=0.33 ap- 
propriate to the sea-level critical latitude we find 
1.13 for R, the ratio of the sea-level vertical 
intensity of the secondaries at 41° to the cor- 
responding value at the equator. This value of R 
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is quite satisfactory. The g primaries have a 
larger latitude effect than the secondaries as 
shown by (55); but, with the value of & already 
taken in our former illustrations, the g primaries 
are entirely masked by their secondaries. 

Under the foregoing conditions there would be 
no variation of latitude effect with altitude until 
the altitude corresponding to x=9.4 already 
cited, i.e., the altitude 0.6 meter of water, was 
attained. At this altitude, a critical latitude 
would appear for the soft component, and would 
gradually work its way towards increasing lati- 
tude with increasing altitude. Thus, for altitudes 
above 0.6 meter of water, the latitude effect 
would increase with further increase of altitude as 
the soft component assumed more and more of 
the dominant role. 

If only one is justified in assuming for the soft 
component a critical latitude different from the 
known sea-level critical latitudes, but one which 
only makes its appearance at altitudes greater 
than 0.6 meter, then I favor the story just 
given in preference to that developed in the previ- 
ous part of this paper, although the only differ- 
ence is in the choice of the constants. The funda- 
mental theory remains the same for both; and is 
exhibited better in the general form in which we 
have presented it in the body of this paper. 

A recent developmeni by George Pfotzer. Since 
the foregoing paper was presented at the Harvard 
Tercentenary Conference, a paper by George 
Pfotzer in the August, 1936 issue of the Zeit- 
schrift fiir Physik has come to my attention." In 
this paper, the author following the lines indi- 
cated by B. Gross, develops the consequences of 
my 1934 theory,’ ' which theory, as already 
stated, resulted in an exponential law of apparent 
absorption and is based upon the law of energy 
loss specified in Eq. (2), combined with the 
assumption that the number of secondaries ac- 
companying a primary increases with the energy 
of the primary. It is convenient to take this 
opportunity to discuss some of the points in 
Pfotzer’s paper for purposes of comparison. In 
the development which I formerly published, the 
law of the increase was one of proportionality 
between secondaries and the first power of the 
primary energy; and, this is the form used by 
Pfotzer. Incidentally Pfotzer also gives the form 


16 G. Pfotzer, Zeits. f. Physik 102, 41 (1936). 
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of energy distribution stated here in Eq. (19) 
which I have derived as a particular case of the 
more general form developed in (13).!7 

In the development of the present paper I had 
started with the hope of utilizng the law in 
which secondary production was proportional to 
the primary energy, but was driven from the 
simplicity of this assumption by the necessity of 
providing for a soft component absorption co- 
efficient equal to 0.5 per meter of water both 
above and below the critical latitude as de- 
manded by Millikan’s observations.'' As already 
demonstrated, the necessity for this provision 
demanded that s\=4u,, where s is the exponent of 
the primary energy in the expression of the 
number of secondaries accompanying a primary. 
With the value \=0.25 which was derived from 
the energy 6X10° ev for penetration at the sea- 
level critical latitude and from a value of a@ falling 
within the limits specified by C. Anderson,® s 
became equal to 2; and, in accordance with (39) p 
became equal to 3 which gave rise to a law of 
approximately the inversed cube of the energy 
for the energy distribution function. Such a law 
is used by Pfotzer, but purely on an empirical 
basis. In this form the theory gave considerable 
difficulty in providing for the proper value of the 
latitude effect without invoking a situation in 
which the calculated energy per unit range of the 
secondaries was too small. I felt the necessity of 
assuming a smaller value of s, but the value of a 
stood in the way, since it controlled the value of 
\. The situation became eased when it was 
pointed out to me by Dr. C. G. Montgomery that 
the values quoted by Anderson for the equivalent 
of a should really be divided by the density of air. 
It was now possible to choose between the ex- 
perimentally allowable limits a value of a which, 
with 6X10° ev for the energy of penetration of 
the atmosphere gave s=1.5. This gave a much 
more satisfactory story as regards the secondary 
production in relation to the latitude effect and, 
incidentally, it led to an energy distribution 
function inversely proportional, approximately 
to the power 2.5 of the energy (see Eq. 38). The 


17 The special form (19), but not the general form (13) 
quoted by Pfotzer from B. Gross, Physik Zeits. 37, 12 
(1936). In Gross’ reference to my paper he quotes J. 
Frank. Inst. 219, 97 (1935). A more complete account of 
the theory concerned is, however, given in Phys. Rev. 48, 
641 (1935). 
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more radical step involving a return to s=1 
through the choice of a value of \ as high as 0.5 
has been referred to above, and necessitates the 
assumption of two critical latitudes, one of which 
does not appear, however, at sea level. 

Pfotzer, on the other hand, uses a very large 
value for u,, a value equal to 0.85 per meter of 
water, and he chooses a very large value, 3 X 10° 
for a, which value he obtains by empirical 
adjustment to fit the experimental data. 

Pfotzer does not place upon his development 
the requirement of equality of a soft component 
absorption coefficient above and below the 
critical latitude. Had he done so, then, quite 
apart from the magnitude of all of the constants 
concerned, the logarithmic forms (36) and (37) 
would have resulted inevitably, and would have 
demanded the assumption of an upper limit to 
the primary energy. 

Another point of difference between the present 
paper and that of Pfotzer lies in the fact that in 
the former we have provided for a variation of 
intensity with latitude and for a variation of the 
latitude effect with altitude, which phenomena 
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placed very stringent requirements upon the 
theory. Pfotzer obtains the sea-level altitude 
effect from the hard component, but does not 
provide for an increase of latitude effect with 
altitude. 

Pfotzer attributes the maximum in the in- 
tensity-altitude curve to a range phenomenon of 
the secondaries and determines from it the 
secondary range at these altitudes. In this, his 
conclusions are substantially similar to those 
which I have reported here and on various former 
occasions,'* the calculated range being equivalent 
to about 82 centimeters of water. 

In conclusion, I wish to express my apprecia- 
tion of the services of Dr. and Mrs. C. G. Mont- 
gomery who have given me considerable assist- 
ance, necessitating careful discrimination, in the 
numerical calculations, and who have checked 
the mathematical calculations. 

BW, F. G. Swann, “Cosmic-Ray Measurements.” 
Presented at Washington, May 1, 1935, as part of a 
symposium on the 1935 National Geographic U. S. Army 
Air Corps Stratosphere Flight. See also W. F. G. Swann, 
G. L. Locher, W. E. Danforth, C. G. and D. D. Mont- 


gomery, National Geographic Society Contributed Techni- 
cal Papers, Stratosphere Series, No. 2 (1936). 
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Cosmic Rays as Electrical Particles* 


ARTHUR H. Compton, University of Chicago, Chicago, Illinois 
(Received October 12, 1936) 


Positive evidence that the primary cosmic rays consist 
of electrical particles is drawn from three types of experi- 
ments: 

1. Latitude and directional asymmetry effects. Clay 
finds 83 percent as intense ionization at the equatorial 
minimum as in high latitudes. Of the remainder, Rossi's 
directional experiments show that about 12 percent at 
least is due to positively charged particles. Corresponding 
to the 73 percent “nonfield sensitive’’ remainder at sea 
level, at high altitudes the remainder must be less than 
20 percent, perhaps no more than 2 percent, of that ob- 
servable in polar regions. An energy distribution analysis, 
following the method of Zanstra, but using new latitude 
effect data collected on the Pacific Ocean in collaboration 
with R. N. Turner, shows a continuous energy distribution 
of the primary cosmic-ray particles between 0.9 and 
1.9 10" ev, and indicates the electrical particle origin of 
a large part, very possibly the whole of the ionization. 


*Based upon a paper delivered at the Tercentenary 
Conference of Arts and Sciences at Harvard University, 
September 8, 1936. 


2. Coincidence experiments. Auger, Street and their 
collaborators have proved that most of the multiple coin- 
cidences observed with counter tubes are caused by single 
high energy ionizing (hence electrical) particles. Experi- 
ments by Rossi and Hsiung show that these coincidence 
producing particles are not secondaries, but originate be- 
yond the atmosphere. Likewise latitude effect experiments 
by Johnson and absorption experiments by Rossi indicate 
that the shower producing radiation is produced by elec- 
trical primaries. Cloud chamber studies show that almost 
all of the observed cosmic-ray ionization is due either to 
particles of the coincidence type or to showers, and is 
hence ascribable to electrical primaries. 

3- Galactic rotation effect. A directional asymmetry 
of cosmic rays ascribable to the motion of the earth with 
the rotation of the galaxy, seems to be established by 
sidereal time variations very recently reported for the 
northern hemisphere by Illing and for the southern hemi- 
sphere by Schonland, Delatizky and Gaskell, and by a new 
provisional observation by Compton and Turner of a 


~~ difference of 0.6 percent between the northern and southern 
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hemispheres. These results are in quantitative agreement and his collaborators have been ascribed to primary 
with the approximate predictions of Compton and Getting photons. A more detailed examination of the photographs 
based upon electrically charged primaries; but the observed in question, however, seems to show that these shower 
diurnal effect is only about }7th as great as is to be ex- _ producing rays are definitely of local and hence secondary 
pected if the primary cosmic rays are photons. The data origin. The arguments that have been raised by Bowen, 
are thus difficult to reconcile with more than a few percent Millikan and Neher against the explanation of cosmic- 
of primary photons. ; ray ionization as due to electrically charged primaries 
Two lines of positive evidence for photon primaries . : - 
4 es or photon primaries are all found to have a straightforward answer if we 
have been put forward. 1. The flattening of the altitude *s e i ‘ : 
atic et : ; : assume the correctness of current theories which make 
vs. ionization curves at very high altitudes, once inter- . . - : 
ema : : photon production the chief mechanism whereby a cosmic- 
preted as indicating the gradual growth of secondary a . Apgat 
gig : ; ray electron loses its energy. The conclusion is drawn that 
ionization as primary photons enter the atmosphere, has 
been shown by later balloon flights at different latitudes 
to be a geomagnetic phenomenon, indicating that the 


primary photons are not responsible for any of the major 
features of the cosmic rays as we now know them. It seems 


corresponding portion of the radiation is electrical in "e€cessary to assume the presence of both positive and 


character. 2. Certain high energy nonionizing shower negative electrons among the primaries. There is also some 


producing rays shown in cloud photographs by Anderson _ suggestion of the existence of protons. 


N several recent addresses' I have presented Permit me first to present the positive evidence 


evidence which seemed to justify the con- that electrically charged particles constitute the 
clusion that cosmic rays consist almost exclu- major part of the primary cosmic rays: 
sively of electrically charged particles, and that 

electrically neutral primary rays, if present, are 1. 
responsible for so small a fraction of the observed It is well recognized that the geographical 
cosmic-ray effects as to be indistinguishable by variations in cosmic rays follow the earth’s 


EVIDENCE FROM THE LATITUDE EFFECT 


our present experiments. A similar view has been magnetic field so closely that the effect must be a 
expressed also by many others.* More recent data magnetic one. This means that a significant 
have seemed to me only to reinforce this poasee. component of the cosmic rays observed near sea 
There have however been put forward from time level is due to electrically charged primaries ; for 

es — x 2 ae > ails : apap, sii, gi 
to time weighty arguments” for the contrary — it has been shown’~? that the geographical effect 
view that the major part of the cosmic rays is is quantitatively explicable on this basis ac- 
electrically neutral. As far as I am aware, nO cording to the theories of Stormer! and Lemaitre- 

g 
attempt has been made to answer these argu- Vallarta," whereas it cannot be explained in 
“ . . . < da, : c >, c 

~ 5 c 3 : se Ss c - oO . ’ ° » 

ments in detail. It is ~ or = this paper to terms of the effect of the earth’s magnetic field 
, r > Sta ance y Se . . . 

review the evidence in light of our present on secondary electrical particles produced in the 
knowledge, and to attempt to reconcile the data atmosphere.” !: "3 The question remains as to 

« - . . c Ds ae 


that have been presented. what fraction of the observed cosmic-ray effects 
1a. Guthrie Lecture of the Physical Society, Feb. 1, owes its origin to such electrically charged rays. 
1935. Proc. Phys. Soc. London 47, 747 (1935). b. American Three typical values of the magnitude a the 


Philosophical Society Address, Apr. 20, 1935. Proc. Am. 
Phil. Soc. 75, 251 (1935). c. Address before American sea-level latitude effect are Clay’s datum”4 of 17 
Physical Society and Section B of AAAS, Jan. 1, 1936. a eee F : on ee 
Rev. Sci. Inst. 7, 71 (1936). percent reduction with an unshielded chamber 
2? E.g., (a) W. Bothe and W. Kolhdrster, Zeits. f. Physik between 45° and 0° magnetic latitude, my own 
56, 751 (1929); (b) B. Rossi, Zeits. f. Physik 82, 151 (1933); fi f . Pig : Se hy 
Proc. London Conference on Nuclear Physics 1, 233 ‘figure of 14 percent with a 6 cm lead shield," and 
(1935); (c) W. Kolhérster, Physik Zeits. 34, 809 (1933); Millikan and Neher’s value® of 12.5 percent with 
(d) J. Clay, Proc. Roy. Soc. A151, 202 (1935); (e) T. H. ———_— 
Johnson, Carnegie Inst. Supp. Pub. No. 13 (1935); (f) P. 7 A. H. Compton, Phys. Rev. 43, 87 (1933). 
Auger, Comptes rendus 200, 739 (1935); (g) P. M. S. 8 J. Clay, Physica 1, 363 and 829 (1934). 
Blackett, La Radiation Cosmique (Hermann & Cie., 1935). *R. A. Millikan and H. V. Neher, Phys. Rev. 47, 205 
3H. Kulenkampff, Physik Zeits. 30, 561 (1929). (1935). 
4 Regener, Kramer and Lenz, Zeits. f. Physik 85, 411 10 C, Stérmer, Zeits. f. Astrophysik 1, 237 (1930). 
and 435 (1933). 1G. Lemaitre and M. S. Vallarta, Phys. Rev. 43, 87 
51. S. Bowen, R. A. Millikan and H. V. Neher, Address (1933). 
before London Conference on Physics, October, 1934. 21. S. Bowen, Phys. Rev. 45, 349 (1934). 
(London, 1935), p. 206. Also in Phys. Rev. 46, 641 (1934). 13 P. M. S. Blackett, reference 2 (g), p. III-18. 
®R. A. Millikan, paper before National Academy of 14 A. H. Compton, Trans. Am. Geophys. Union (1933), 


Sciences, Apr. 28, 1936 (unpublished). p. 154. 
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a 12 cm shield. Of these, Clay’s value of 17 
percent is to be preferred for our present purpose. 
This fraction is known to be due to electrically 
charged primaries, leaving 83 percent to be 
accounted for. 

That a part of the remainder is electrical is 
established by directional experiments, made 
with coincidence counting tubes, which show 
that the rays reaching the earth at the equator 
come mostly from the west.'® This, like the 
latitude effect, is due to the action of the earth's 
magnetic field, and shows an excess of positively 
charged particles. At a zenith angle of 45°, 
having filtered out the soft secondaries with 8 cm 
of lead, Rossi'®* finds near the equator 0.76 as 
many rays from the east as from the west. For 
these equatorial rays, therefore, at least 24 
percent of those incident at 45° west, or about 12 
percent of those incident vertically, are due to 
positively charged primaries. Taking this 12 
percent as a reasonable average, there thus 
remain 0.83 X 0.88 =73 percent as an upper limit 
to be accounted for,'® the remaining 27 percent 
having been proved electrical. 

An energy distribution analysis based upon the 
latitude effect shows that these magnetically 
deflectable rays constitute a portion of what 
appears to be a continuous energy spectrum. 
Zanstra has pointed out!’ that, since there is a 
definite lower limit Vo imposed by the earth's 
magnetic field upon the energy of the rays 
reaching the earth vertically at each latitude, a 
measurement of the latitude effect for such 
vertical rays should give their energy spectrum. 
If ¢(V)dV is the ionization due to the vertical 
rays having energies between V and V+dI, 
Zanstra shows that 


o(V)dV = —(dy/dVo)dV, (1) 


where y is the observed ionization due to the 
vertical rays, and V» is the lower limit to the 
energy of the transmitted electrons, which, ex- 


15 E.g., (a) T. H. Johnson, Phys. Rev. 43, 834; 44, 856 
(1933); (b) L. Alvarez and A. H. Compton, ibid. 43, 835 
(1933); (c) B. Rossi and S. de Benedetti, ibid. 45, 214 
(1934). 

16 Using the data from their own latitude effect and 
Korff's directional effect measurements in America, Bowen, 
Millikan and Neher, reference 5, calculate this remainder 
on a similar basis as 0.85. 
17H. Zanstra, Naturwiss. 22, 171 (1934). 
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pressed in electron volts, is* 


Vo=1.92 & 10" cos‘ X, (2) 


\ being the geomagnetic latitude. 

Though latitude effect measurements on the 
vertical rays have been made with coincidence 
counting tubes,'® the observations have sta- 
tistical errors which are too great to be suitable 
for such an analysis. Measurements of the 
distribution of cosmic rays with varying zenith 
angle show, however, that their angular distri- 
bution is within experimental error the same 
at different latitudes.'® Accordingly, ionization 
measurements which average the intensity from 
all zenith angles should give very approximately 


* Note added in proof: The value of V» is given more 
exactly by the recent calculations of Lemaitre and Vallarta 
(Phys. Rev. 50, 493 (1936), Fig. (10). The difference from 
Stérmers formula (2) is however not sufficient to modify 
the conclusions here reached. 

18a. P. Auger and L. Leprince-Ringuet, Nature 133, 138 
(1934); b. J. Clay, Physica 2, 299 (1935). 

‘° Among the best published results are those of Johnson, 
reference 20, in Peru (equator), and of Kolhérster and 
Janossy, reference 21, in Germany. They find, respectively, 
for the ratio of the intensity at zenith angle 45° to 0° the 
values 0.475 and 0.495, each with a probable error of about 
1 percent. Auger and Leprince-Ringuet, comparing this 
ratio at latitudes 0° and 38° with somewhat lower precision, 
find a slight difference in the opposite sense. 

*° T. H. Johnson, Phys. Rev. 45, 584 (1934). 

*1W. Kolhdérster and L. Janossy, Zeits. f. Physik 93, 
111 (1934). 
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Fic. 2. Curve I represents data of Fig. 1 plotted against 
the geomagnetic energy minimum Vo. Curve II is the 
slope of curve I, and represents the energy spectrum of 
the rays (assumed to be electrons) transmitted by the 
atmosphere. 


the relative values for the vertical component of 
the rays. We may thus use in place of Zanstra’s 
Eq. (1), 


(3) 


—dI/dVo, 


o(V)= 
where J is the cosmic-ray intensity as observed 


with an ionization chamber. 


H. 
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Thanks to the of the Canadian 
Australasian Steamship Company, Captain Hill 
of the R.M.S. Aorangi, and especially the very 
skillful help of First Officer R. N. Turner, we 
have been able to obtain a series of repeated 
measurements the the 
Vancouver, 


courtesy 


of latitude effect over 


route from Canada, to Sydney, 
Australia, using one of the Carnegie Institution's 
model C meters” shielded with 12 cm of lead. The 
results of 4 trips are shown in Fig. 1, and repre- 
sent the most suitable data for our modified 
Zanstra analysis which are now available. Re- 
calling that each geomagnetic latitude corre- 
sponds to a given energy minimum V9, these data 
may be averaged and plotted against Vo, as in 
curve I of Fig. 2. For energies lower than 0.2 the 
curve is drawn flat, in accord with Clay’s result 
that north of \=54° no change as great as 0.6 
percent occurs.” When the slopes of this curve 
are taken, we obtain from Eq. (3) the values of 
¢o(V) shown in curve II. 

For values of Vo below 0.9 10!° the effect of 
the earth’s atmosphere becomes evident in 
limiting the rays that are received to a higher 
than minimum admitted by the 


energy the 


2 Cf. A. H. Compton, E. O. Wollan and R. D. Bennett, 
Rev. Sci. Inst. 5, 415 (1934). 

237. Clay, C. G. t’Hooft and P. H. Clay, Physica 2, 
1033 (1935). 
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Fic. 3. Data of Fig. 1 averaged for both hemispheres and plotted logarithmically 
against Vo, showing that the magnetic latitude effect is limited only by the finite magni- 
tude of the limiting energy V» corresponding to the earth’s magnetic field. 
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magnetic barrier. The form of the curve for the 


lower energies strongly suggests two components, 


one with a limiting energy imposed by the 


atmosphere of about 8 X 10° ev (+electrons?) and 
the other of about 2.5 10° ev (protons?). There 
are several lines of evidence supporting this 
interpretation. It will, however, be preferable to 
postpone discussion of these interesting details 
until the further data now being collected are 
available. 

For our present purpose the important point is 
that for energies over 0.910! ev, where the 
atmospheric absorption is not a determining 
factor, there appears a continuous energy distri- 
bution as far as 1.92 ev, beyond which the 
method is no longer applicable.** The trend of 
the curve indicates, however, that electrical 
particles of much higher energies are present in 
considerable quantity. 

In order to extrapolate to higher energies, a 
more suitable representation of the data is to 
plot log I against log Vo, as in Fig. 3, where the 
points represent the average measured intensi- 
ties. For energies greater than 0.9 < 10" ev it will 
be seen that the points fall within experimental 
error upon a straight line. If this line were to 
continue indefinitely, it would mean that for 
infinite limiting energy Vy) the percentage of 
transmitted cosmic rays J would be zero, i.e., all 
the cosmic rays would be proved electrical. In 
other words, this quantitative analysis indicates 
that the magnitude of the latitude effect is 
limited only by the strength of the earth's 
magnetic field, and if that field were very great, 
the cosmic rays might well be completely eX- 
cluded near the equator. 

The latitude effect at high altitudes has not 
been studied in sufficient detail to carry through 
a similar analysis. Its rapid increase with altitude 
is however well known. Fig. 4 shows typical 
results of recent high altitude measurements. 
Millikan and Neher have recently reported® that 
their airplane flights in the Philippines give 
results closely similar to those here shown for 
Peru, whereas Clay’s values in Java are presented 
with reservations. These data indicate a ratio 
near the top of the atmosphere of at least 4 to 1 

*4* This conclusion is similar to that reached by Clay, 


reference 18 (b), from his latitude effect experiments with 
counter tubes. 
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between 52° and the equator. There seems also 
little question but that near the poles the 
intensity at great altitudes should be con- 
siderably greater than at 52°.*° Moreover, as in 
the case of the latitude effect at sea level, 
directional experiments show that much of the 
rays reaching the equator must be electrical. We 
may accordingly conclude that of the cosmic-ray 
ionization at very high altitudes, not more than 
a few percent can be due to neutral primaries. 


2. EVIDENCE FROM COINCIDENCE EXPERIMENTS 


Perhaps the most striking proof that the 
coincidences observed with multiple counters in 
line are due to single high energy particles 
traversing them all is supplied by the recent 
experiment performed independently by Auger 
and Ehrenfest®*® and by Street, Woodward and 
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Fic. 4. High altitude measurements of cosmic rays at 
different latitudes. 52°, Compton-Stephenson- Millikan ; 
42°, Millikan-Bowen; Peru, Millikan-Neher; Java, Clay. 

25 Cf., however, M. Cosyns, Nature 137, 616 (1936), 
whose balloon flights suggest that the latitude effect at 
high altitudes may end sharply at geomagnetic latitude 
49° in Europe. 

2° P. Auger and Ehrenfest, Comptes rendus 199, 1609 
(1934). 
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Stevenson,”’ and illustrated in Figs. 5 and 6. 
Cloud photographs are taken when 3-fold coinci- 
dences are produced by rays traversing 45 cm of 
lead. With this arrangement, 90 percent of the 
photographs record single straight tracks. This 
is clear that 
produced by highly penetrating particles. Ex- 


evidence such coincidences are 
periments by Rossi** have shown that about half 
of these particles will penetrate a lead screen 
equivalent in mass to the atmosphere, which is 
about the same rate of absorption as is observed 
with the cosmic rays as measured by an ioniza- 
tion chamber. This strongly suggests that the 
penetrating particles are the primary cosmic 
rays, a confirmation of Bothe and Kolhdérster’s 
classical argument.’ 

Rossi?” and Hsiung?’ hay e performed in different 
forms an experiment which shows that these 
penetrating particles are not secondaries pro- 
far 


ducible at sea level, but must come from 


above the apparatus. Thus in Fig. 7, if pene- 
trating secondaries were excited in 20 cm of lead, 
the counts with arrangement C should be greater 
than in case B. After making allowance for the 





es ¢ 
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Fic. 5. Arrangement of coincidence counters, C;, Co, C. 
and expansion chamber Ch for photographing trails of 
penetrating cosmic rays (Street-Woodward-Stevenson). 


*7 J. C. Street, R. H. Woodward and E. C. 
Phys. Rev. 47, 891 (1935). 

25 B. Rossi, Naturwiss. 20, 65 
82, 151 (1933). 


Stevenson, 
1932); Zeits. f. Physik 


Rev. 46, 653 (1934). 


2? D. S. Hsiung, Phys 


R H. 
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Fic. 6. Trails of two cosmic-ray particles after traversing 





40 cm of lead (Street-\Woodward-Stevenson). 


slight effect of the well-known soft secondary 
rays, both Rossi and Hsiung find the coincidence 
rate in cases B and C identical within experi- 
mental error. Thus all the coincidence producing 
particles originate above the apparatus. 

It might be supposed that these penetrating 
particles are, however, secondaries excited high 
in the atmosphere by readily absorbable photons. 


The 


coincidence producing particles have been found 


This suggestion is ruled out by two facts. (1 


to show the latitude effect in substantially the 
same way as does the total cosmic-ray ioniza 
tion,'* and also to exhibit directional asymmetry 
in the earth’s magnetic field at the equator. 
They are thus at least in large measure due to 
high energy primary electrical particles, and our 
knowledge of their high penetrating power can 
leave little doubt but that they are themselves 
these primary particles. (2) We shall see below 
that very high energy photons and electrons of 
the same energy should be absorbed at approxi- 
mately the same rate. Thus if the high energy 
coincidence producing electrons are due to pri- 
mary photons, some of these photons should 
penetrate the atmosphere, in intensity compa- 
rable with the secondary electrons, and excite 
more high energy electrons in the lead block of 
the Rossi-Hsiung apparatus. That this does not 
occur may be taken to mean that the fraction of 
the coincidence producing particles which are 
secondaries of high energy photons is zero, with a 
probable error of about 3 percent.*° 

In addition to these coincidence producing 
particles, a ‘shower producing radiation” is like- 

39 If, as Bowen, Millikan and Neher assume, the photons 
are much more penetrating than the electrons of the same 
energy, this argument becomes the more cogent; for the 
conditions for producing high energy electrons in the lead 
block are then ideal 
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wise recognized.** This may be studied either 
with counting tubes out of line, or with a cloud 
chamber. The evidence strongly favors the view 
that this shower producing radiation consists of 
high energy photons. Johnson** has shown, 
however, that this radiation also is subject to 
latitude effect, which means that the photons 
must themselves be produced by primary elec- 
trical particles.” It would appear that nearly all 
of the cosmic-ray tracks observed in cloud 
chambers are of either the high energy coinci- 
dence producing type or the shower type. Since 
both of these are ascribable to primary electrical 
particles, very little room is left for any ionization 
that may result from primary photons. 

It should be added that cosmic rays in deep 
mines are measurable with multiple coincidence 
counters, which, as Kolhérster and Clay have 
pointed out,’ means that the rays there observed 
are corpuscular in the same sense as are those at 
sea level. 


3. EVIDENCE FROM THE TRANSLATIONAL EFFECT 


There is a third type of evidence regarding the 
nature of the primary cosmic rays which comes 
from the magnitude of the effect on the intensity 
of cosmic rays due to the earth’s rapid motion 
with the rotation of the Milky Way. The 
experimental data are now adequate to state 
with considerable confidence that this effect 
exists, and is of the magnitude predicted upon the 
assumption of electrically charged rays. The 
observed effect is, however, much too small to fit 
the prediction based upon photons as the primary 
rays. 

It will be recalled that according to the 
astronomers’ recent findings the rotation of the 
galaxy carries the earth in the direction of the 
constellation Cygnus, about 45° north, at a speed 
of about 0.001 that of light. For rays such as 
cosmic rays going with approximately the speed 
of light, this produces a Doppler effect difference 
in frequency of 0.2 percent between the rays 
striking the front and the back of the earth, if the 


origin of these rays is assumed to be at rest 


31 (a) Cf. B. Rossi, reference 2 (b); (b) C. D. Anderson, 
R. A. Millikan, S. Neddermeyer and W. Pickering, Phys. 
Rev. 45, 352 (1934); (c) T. H. Johnson, Phys. Rev. 47, 
318 (1935). 

32 Rossi (reference 2 (b)) reaches the same conclusion on 
the basis of absorption measurements. 


wn 


outside of the galaxy. The corresponding Doppler 
effect in intensity at the surface of the atmosphere 
would be 4 times as great, or 0.8 percent. When 
account is taken of the absorption by the 
atmosphere, and the fact that the measurements, 
being made with an ionization chamber, include 
rays from all angles, Dr. Getting and I calculate 
about 1.0 percent difference between the front 
and the back of the earth.* This, of course, 
assumes that if the earth were not moving the 
intensity would be the same from all directions. 

If the primary cosmic rays consist of photons 
or other neutral particles, this full effect should 
be observed. If they are electrically charged, 
according to Stérmer’s theory the particles will 
be so strongly deflected by the earth’s magnetic 
field that the effect will be much reduced. A 
comparison of the predicted and the observed 
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Fic. 7. Hsiung's coincidence experiment. 


33 A. H. Compton and I. A. Getting, Phys. Rev. 47, 
817 (1935). The value there predicted is 1.2 percent 
difference between the front and the back of the earth. 
The slightly lower value here given involves a minor 
correction to our published calculation. 
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TABLE I. Test of the translational effect. 
Percent variations. 
Theor 
Positron 
Photon Negatron 
Primaries Primaries | Experiment 
Latitude 47°N 
Amplitude of diurnal 
Variation 0.25 0.02-0.10 0.03 


Phase of diurnal varia- 


tion 20h 50m 


20h 40m | 20h 40m 


Latitude 34°S 
Amplitude of diurnal 
Variation 0.28 0.03-0.12 0.04 
Phase of diurnal varia- 
tion 20h 40m | 20h 40m | 21h 45m 
N —S Difference 0.5 0.5 0.67 
+ Triennial mean, Illing. 
X Triennial mean, Schonland, Delatizky and Gaskell 
= Provisional value, Compton and Turner 
values is shown in Table I for both the as- 


sumption of photon primaries and of positron- 
negatron primaries. This table shows the data 
that have been reported within the last few 
months. It may be noted that whereas the effect 
calculated for electrical particles is rather uncer- 
tain because of the unknown magnitude of the 
geomagnetic effect, if the cosmic rays are photons 
the only escape from an effect of about the 
predicted magnitude would seem to be the 
difficult assumption that there occurs a natural 
anisotropy in the source of the rays just sufficient 
to balance the translational effect. 

A suggestive confirmation of the prediction was 
immediately possible using data obtained by 
Hess and Steinmaurer during the year 1932. 
Since that time IIling** has published, in June of 
this year, a reduction of Hess’ measurements for 
the three years 1932, 1933 and 1934, corrected 
for barometer changes and plotted against 
sidereal time. The result, as indicated in the first 
row of Table I, shows a variation with precisely 
the predicted phase, and with an amplitude 
which lies within the range estimated for electron 
cosmic rays, though only 12 percent of that 
calculated for photon rays. 

Illing left the existence of the effect open to 
some experimental doubt, because of possible 
seasonal variations. Just now, however, a similar 
summary of three years’ measurements in the 


34 W. Illing, Terr. Mag. and Elec. June (1936) 
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southern hemisphere has appeared, by Schonland 
and his collaborators at Capetown,®® which 
shows an effect of the same order of magnitude 
and with almost exactly the same phase. This 
disposes of any possibility of a seasonal variation, 
since this should be opposite in phase in the two 
hemispheres. We may thus consider the galactic 
rotation effect to be real, and its magnitude to 
fit the rough estimates based upon the electron 
hypothesis of cosmic rays. 

Our studies of the intensity of cosmic rays on 
the Pacific Ocean should eventually give us a 
good experimental value of the difference in 
intensity between the northern and southern 
hemispheres. As yet our results are only pro- 
visional, because the southern end of the journey 
is not far enough south to reach the polar plateau 
of cosmic-ray intensity. Basing our estimate, 
however, upon the observed intensity at the 
first knee in the northern and southern curves, 
we find a greater intensity in the north by 
about 0.6 percent—very close to the theoretical 
difference. 

This tentative verification that 
tional effect is of the full magnitude given by the 
theoretical calculation makes it difficult 
indeed to reconcile the small diurnal variation 
with the assumption that photons constitute 
more than a few percent of the primary cosmic 


the transla- 


very 


rays 
4. ARGUMENTS FOR PRIMARY PHOTONS 


The most effective defense of the hypothesis 
of primary photons which has come to my 
attention is that of Bowen, Millikan and Neher,® 
presented before the London Congress of Nuclear 
Physics in 1934. They admit the electrical 
character of the primary rays responsible for 
about 15 percent of the ionization observed at 
sea level, as shown by latitude effect and direc 
tional experiments. They however consider prac- 
tically the whole of the remaining ‘“‘nonfield 
sensitive” portion, estimated as responsible for 
85 percent of the sea-level ionization, to be due 
to primary photons. More recent statements by 
Professor Millikan® indicate that this view has 
remained essentially unaltered. 


Nature 138, 325 (1936). 
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Fic. 8 AND Fic. 9. Cloud photographs of showers (Anderson, Millikan, Nedder- 
meyer, Pickering) that have been used as evidence of incoming primary photons. 
rhe showers are interpreted in the text as due to secondary photons. 


1) The one positive argument put forward for 
primary cosmic-ray photons is the citation® of 
two of Anderson’s cloud photographs, Figs. 8 and 
9, which show showers of positrons and negatrons 
that must have been produced by photons 
having in some cases total energies exceeding 
2.5 10° ev.** “There seems to be no escape,” 
they comment, “from regarding these showers as 


resulting from incoming photons themselves.”’ 


° Cf. reference 31 (b), Figs. 9 and 11. 


We would suggest the alternative interpreta- 
tion that 


showers 


the photons which produce these 


are themselves secondaries from an 
event which occurred not far above the cloud 
chamber. In Fig. 8 we note that at least two of 
the showers, the high energy pair from the top 
of the chamber and the large shower from the 
central lead plate, have total energies of the same 
order of magnitude, >5 x 108 ev. If it is assumed 


that these are caused by the same photon in 
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Fic. 10. Cosmic-ray intensity vs. latitude curves for, A, 
high magnetic latitudes, B, magnetic equator, C, difference 
between A and B, showing that equatorial rays B are 
absorbed in the same manner as are the magnetically 
deflectable rays C (Compton and Stephenson). 


successive ionizing events, it indicates a rate of 
absorption much too great for a ray which can 
penetrate the atmosphere. If on the other hand 
they are due to different photons, their identity 
in time and direction and their small lateral 
separation shows that these photons must have 
originated in the same event not far above the 
chamber. In either case their secondary origin 
appears to be necessary. Similarly an examina- 
tion of Fig. 9 shows at least 3 points of origin 
of the showers, close together in the central 
plate. Two of these showers are of the same 
energy. Their 


large order of magnitude of 
direction shows that no one can be considered as 
the parent of the others; it is clear that they are 
all due to a close bundle of nearly parallel, non- 
ionizing rays, presumably photons. Their close 
proximity indicates that their origin must have 
been only a short distance above the expansion 
chamber. Here also a secondary origin of the 
photons is therefore clearly indicated. The evi- 
dence presented for the occurrence of primary 
photons is thus unsatisfactory 

2) All the remaining arguments are attempts 
to show that incoming electrical particles are 
incapable of producing the observed cosmic-ray 
effects. Thus Bowen, Millikan and Neher state 
that the absorption of the high energy electrons 
“cannot possibly yield an altitude ionization 
resembling that found experi- 


curve remotely 
mentally from sea level up to 4.5 m.”’ This they 
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emphasize as ‘the most general and the most 
cogent argument against the attempt to explain 
the nonfield sensitive portion of the cosmic-ray 
high 
that 


ionization as due primarily to incoming 


energy electrons.” They suppose rather 


absorption of such high energy electrons should 


be of the “range’’ type instead of the ‘‘expo- 
nential’” type. Such arguments are ineffective 
when it is recognized that the main part of the 
energy loss by high energy electrons is through 
the production of photons (as x-rays are produced 
by cathode rays), and that this type of energy 


loss leads necessarily to an appreximately ex- 


ponential diminution of ionization along the 
path of the electrons. 
That high energy electrical rays are in fact 


absorbed exponentially was shewn some years 
ago by our high altitude measurements of cosmic 
rays at latitudes.*” Thus in Fig. 10, 


curve A represents the absorption in the atmos- 


different 


phere of the rays received at high magnetic 
latitudes, and B that of the rays at the equator 
A and B, and 


thus shows the absorption of electrically charged 


Curve. C is the difference between 


rays of the type so affected by the earth’s mag- 
netic field that they cannot reach the earth at 
the equator. The close similarity of curve C, 
which includes only electrically charged rays, to 
curve B, which represents the rays transmitted 
by the earth’s field at the equator, shows that 
electrically charged particles are absorbed in the 
manner characteristic of cosmic rays. As is well 
known, a curve of type C, being concave up 
wards, is expressible as a spectrum of expo- 
nentially absorbed components. In order to fit 
the absorption co- 
from 


values of 
to extend at 


this the 
efficients would 
0.035 to 0.12 per cm of mercury, corresponding 


to the electrical particles, perhaps of different 


curve 


need least 


natures, included in the energy range of from 
210° to 2X10" ev. This would seem to be an 
adequate answer to the argument against pri 
mary electrons quoted above. 

According to Bethe and Heitler’s theory of 
absorption,** for cosmic-ray energies there is a 
striking similarity in the dissipation of energy by 
photons and electrons. This is shown graphically 

7 A. H. Compton and R. J. Stephenson, Phys. Rev. 45, 
448 (1934). 

’H. Bethe and W. Heitler, Proc. Roy. Soc. London 
A146, 83 (1934). 
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in Fig. 11, adapted from their paper. We may 
note that for electrons it is the electron collisions 
that result in ‘‘range’’ absorption. For energies 
greater than about 200 mc? or 10° ev the radiation 
process is the predominant factor, and the ab- 
sorption should correspondingly be of the expo- 
nential type. This includes all of the important 
cosmic-ray energy range. Similarly, for energies 
than 30 mc* pair production 
becomes the predominant method of 
dissipation for photons. It is especially to be 
noted as shown by the solid line, that while in 
the low energy range electrons are stopped much 
more quickly than photons, for energies greater 
than 1000 mc, or 5108 ev, the predicted ab- 
sorption rate for electrons is but little greater 


greater about 


energy 


than for photons. 

These major results are common also to the 
extension of Bethe and Heitler’s theory proposed 
by Oppenheimer,*® and to the alternative theory 
of Swann.* The rate of energy loss by electrons 
in lead as measured by Anderson has substanti- 
ated the predictions up to 800 mc?, and shows at 
least no major departures up to 2000 mc*. The 
theory is likewise consistent with what little is 
experimentally known regarding the absorption 
of high energy photons. 

With this theory in mind, the remaining ob- 
jections® to the primary electron hypothesis for 
cosmic rays are immediately answerable: (3) 
The photons which occur with energies up to 
2.5 10° ev, as shown by photographs of cosmic- 
ray showers, may readily be interpreted as 
secondaries when we remember that primary 
electrons are found with energies more than 10 
times as great, as proved by equatorial direc- 
tional experiments. (4) The statement that the 
hypothesis of primary electrons would leave very 
few observable electrons at sea level in the inter- 
mediate energy range 3x 108 
310° ev leaves out of account the tertiary 
electrons excited by the secondary photons, 
which would be expected to have energies within 
this range. (5) The claim that photons are 100 


between and 


times as penetrating as electrons of the same 
energy is inconsistent with our present know!- 
edge. (6) ‘“Exponential’’ absorption for rays of 


39 J. R. Oppenheimer, Phys. Rev. 47, 44 (1935). 
40 W. F. G. Swann, Phys. Rev. 46, 828 (1934); 47, 250 
(1935). 


E 





1129 


CTRICAL PARTICLES 


cosmic-ray energy is not to be considered a 
unique characteristic of photons. 

(7) The increase in penetrating power with 
depth at the equator, where the low energy 
primary electrons are absent, has been con- 
sidered® to favor the primary photon hypothesis. 
It may be noted that Bethe and Heitler’s theory, 
resulting as it does in a practically uniform ab- 
sorption coefficient in air or water for all photons 
of energy greater than 10’ ev, provides no means 
whereby photon radiation can be hardened as it 
traverses the atmosphere. Swann has suggested *® 
a modification of the theory which would result 
in increased penetrating power of electrons on 
passing An in- 


terpretation is the assumption that the elec- 


through matter. alternative 
trical particles are of two or more types, such as 
and should 


very little radiation, and should be much more 


electrons protons. Protons excite 
penetrating. They should thus remain after the 
electrons are absorbed. There are some difh- 
culties with this suggestion, though there are 
also several lines of supporting evidence. In any 
case the difficulty of accounting for the harden- 
ing of the rays on traversing matter is not re- 
moved by substituting photons for electrons as 
primary rays. On the other hand, direct coinci- 
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Fic. 11. Effective cross section (units Z*r?/137) of water 
molecules for the production of radiative and electron 
collisions by high energy electrons (broken lines), and pair 
production and Compton effect collisions by photons 


(dotted lines). The solid line gives the ratio of the total 
absorptional cross-sectional areas for electrons and photons 
(adapted from Bethe and Heitler). 
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dence counter experiments by Rossi** and others, 
which record only electrical particles, show a 
continual hardening as the thickness of the ab- 
sorption screen is increased to 100 cm. This 
constitutes direct evidence that electrical rays 
are hardened by filtering through matter. 

(8) The observation by Anderson“ of nearly 
equal numbers of positive and negative particles, 
at least for energies up to 10° ev, is taken® to 
imply that the electrons are secondaries result- 
ing from pair production. This is undoubtedly 
the correct interpretation for most of the lower 
energy particles. Johnson*'* has however given 
apparently conclusive evidence that the 
called ‘“‘shower producing” cosmic rays are due 
to primaries consisting of approximately equal 
numbers of positive and negative rays. This 
evidence consists in showing that these rays are 
subject to the latitude effect and are hence 
electrical, but show very little directional asym- 
metry and must thus include equal numbers of 
positive and negative particles. It is thus unsafe 
to conclude from an equal distribution of posi- 
tive and negative electrons that they are of local 


so- 


origin. 

(9) In earlier discussions*® ‘ the flattening of 
the altitude vs. ionization curve at very high 
altitudes, as observed by Regener, was taken to 
indicate the photon character of the primaries, 
which should not ionize until they had produced 
secondary electrons, and should thus show feeble 
ionization near the top of the atmosphere. 
When, however, the Century of Progress bal- 
loon flight, made at a higher magnetic latitude, 
showed no such flattening, it became clear!* 
that this was a geomagnetic phenomenon, and 
indicated that the corresponding portion of the 
radiation is electrical in character. 


SUMMARY 


These arguments accordingly fail to give any 
support to the photon hypothesis of the primary 
cosmic rays. They receive, however, a straight- 


“1C. D. Anderson and S. H. Neddermeyer, Proc. Int. 
Conf. on Physics (London, 1935), Vol. I, p. 171. 
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forward interpretation on the hypothesis of 
electrical primaries if we assume the correctness 
of the current theories which make photon 
production the chief mechanism by which an 
electron loses its energy. 

On the other hand, near the top of the at- 
mosphere latitude and directional experiments 
show that definitely more than 75 percent and 
probably more than 95 percent of the ionization 
near the poles is due to charged primaries. 

At sea level, the latitude effect establishes the 
electrical origin of 17 percent of the ionization, 
and when analyzed in terms of energy distribu- 
tion indicates a similar origin of a much larger 
part, very possibly the whole, of the ionization. 
Directional experiments also show that at least 
12 percent of the remainder is electrical. The ob- 
served small magnitude of the diurnal changes 
of cosmic rays caused by the earth’s motion with 
the rotation of the galaxy, seems irreconcilable 
with more than a few percent of primary photons. 
Coincidence counter experiments confirm this 
result, showing that all the coincidence produc- 
ing rays, within a probable error of about 3 per- 
cent, are due to primary electrical particles. 

At great depths, the fact that coincidence 
counting tubes give absorption values in good 
agreement with those from icnization measure- 
ments is evidence that the most penetrating 
cosmic rays also are electrically charged particles 

It is rather surprising to note that there has 
not appeared positive evidence of a reliable 
character for any photons whatever among the 
primary cosmic rays It seems certain that they 
are not responsible for any of the major features 
of the rays as we now know them. 

In order to describe adequately the properties 
to the cosmic rays, it seems necessary to assume 
the presence of both positive and negative 
electrons among the primaries. There is also some 
suggestion of the existence of protons. We can- 
not now discuss these interesting details. It 
seems however fair to consider the electrically 
charged nature of the primary cosmic rays as 
one of their important characteristics estab- 
lished by the last decade of cosmic-ray research. 
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Significant improvements of the apparatus for the mul- 
tiple acceleration of light ions, now called the ‘‘cyclotron,” 
have recently been made. The apparatus produces many 
microamperes of deuterons having energies up to 6.3 Mev 
and a few tenths of a microampere of 11 Mev doubly 
charged helium ions. The means of deflecting the high 
speed ions from the accelerating system has also been 


INCE the first experiments' six years ago, 

work in our laboratory on the development 
of methods of accelerating charged particles has 
gone forward steadily with the endeavor of 
evolving methods and apparatus for nuclear 
physics that would find wide usefulness. Almost 
from the beginning, the method of multiple 
acceleration of light ions, wherein the ions cir- 
culate in a magnetic field in resonance with an 
oscillating electrical field, gave most promise, and 
accordingly has been given major attention in 
our laboratory. Four years ago, details of an 
apparatus that generates 1.2 Mev protons were 
described,? and two years ago an account*® was 
given of a much larger apparatus of this type, 
capable of generating a fraction of a micro- 
ampere of 5 Mev deuterons. In the meanwhile, 
significant improvements have been made not 
only in increasing the output of ions to many 
microamperes at higher voltages, and in the 
general convenience and reliability of the appa- 
ratus, but also in the means of withdrawing the 
ions from their spiral paths. It is the purpose of 
the present paper to describe these recent 
advances in design and also to discuss some of 
the questions involved, for some of the changes 
which we have made suggest the possibility of 
further improvements. 


APPARATUS 


A general view of the apparatus, which is now 
called‘ the magnetic resonance accelerator, or, 


1E,. O. Lawrence and N. E. Edlefson, Science 72, 
376-377 (1930). 

2 E. O. Lawrence and M. S. Livingston, Phys. Rev. 40, 
19-35 (1932). 

3 E. O. Lawrence and M. S. Livingston, Phys. Rev. 45, 
608-612 (1934). 

*E. O. Lawrence, E. McMillan and R. L. Thornton, 
Phys. Rev. 48, reference p. 495 (1935). 


improved, with the result that the ions can now be with- 
drawn completely from the apparatus into the air through 
a thin window in the chamber wall. The present paper is 
devoted to a detailed description of the cyclotron, including 
its adjustment and operation, with a discussion of the 
questions involved. 


more frequently, the ‘‘cyclotron,” is shown in 
Fig. 1. The magnet which was described in the 
publication two years ago has not been changed 
in any respect. For some time we have intended 
to enlarge the pole faces beyond the original 273” 
diameter, but this change still remains for the 
future. All of our attention has been given to the 
improvement of the cyclotron at its present size. 
The drawing of Fig. 2 is a general view of the 
chamber with cover removed, showing the ion 
source, the accelerating electrodes, which are 
called duants, or dees, the electrode for deflecting 
the ions from their spiral paths, the rotatable 
beryllium and tungsten target, and the window 
in the chamber wall through which the ions 
emerge into the air. These component parts are 
discussed in the following paragraphs. 


Ion source 


As in the original arrangement, the ions are 
produced in the central region by a beam of 
electrons passing along the magnetic lines of 
force between the dees and ionizing the gas 
filling the chamber at a suitable pressure. For- 
merly, spiral filaments above and below the 
central region were mounted on glass tubes in 
copper pipes which extended from the chamber 
wall to the center, along a radius between the 
dees. These tubular shields extended to the 
center on one side only, thus introducing an 
undesirable asymmetry in the field between the 
dees. This has been corrected® in the present 
apparatus by bringing the filaments to the center 
across one dee, as shown in Fig. 3. Here may be 


*In the Cornell cyclotron this correction has been 
achieved by the symmetrical placing of four filaments 
(M. S. Livingston, Rev. Sci. Inst. 7, 55-68 (1936). 
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Fic. 1. General view of the cyclotron. 


usually 53} turns of 30 mil tungsten with an 
over-all diameter of about }’’, are mounted by 
clamping to water-cooled copper blocks. The 
filament support is enclosed in a copper shield 


seen two filament supports, one above and one 
below one of the dees (which in this photograph 
has been removed). The flat spiral filaments, 





‘ Fic. 3. Photograph of a cyclotron chamber with cover 
plate and one dee removed. In this apparatus the ions are 


Fic. 2. Perspective drawing of the cyclotron chamber with not withdrawn completely but bombard targets within the 


iron cover plate removed. chamber wall. 








eo 





ne TE, rer. 2» 











APPARATUS FOR MULTIPLE ACCELERATION 1133 


(because of the radio frequency potentials on 
the dees) and the whole assembly is made as thin 
as possible for maximum clearance between the 
dee and the filament structure. In practice, it has 
been our experience that }” vertical or }” 
horizontal clearance is enough to prevent trouble 
from breakdown discharges. 

The filament support is mounted on a sylphon 
arrangement which permits moving the filament 
to various places in the central region. Theoreti- 
cally, the ion source should be near the center of 
the cyclotron, inasmuch as the ions start their 
spirals on paths of small radii of curvature in 
the oscillating electric field between the acceler- 
ating electrodes. In practice, however, it is 
usually found that the optimum position of the 
filament is not at the center but off the center 
one way or another. This arises from the circum- 
stance that in general the centers of the suc- 
cessive circular paths migrate from the geo- 
metrical center of the apparatus as a result of the 
lack of complete circular symmetry of the mag- 
netic field. This migration of the paths of the ions 
can be changed by correcting the asymmetry of 
the field. This can be done by inserting suitable 
pieces of iron (which we call ‘‘shims’’) between 
the cover plates of the vacuum chamber and 
the magnetic pole faces, but oftentimes it is 
much simpler to make the correction by shifting 
the point at which the ions begin their spiral 
paths. 

The whole filament assembly is bolted and 
waxed to a plate on the side of the chamber wall, 
and is made conveniently demountable for 
filament replacement. Water-cooling pipes pass 
around the filament shield and also are provided 
for the wax joints. The filament electrodes are 
insulated from the surrounding shield by glass 
tubes and emerge from the sylphon assembly 
through a rubber plug, the arrangement being 
made vacuum-tight by wax. The insulation is 
adequate for several thousand volts, but in 
general the filament is maintained at a negative 
potential not greater than 1000 volts with 
respect to the surrounding shield and the walls 
of the vacuum chamber. 

Although there are two filaments (one above 
and one below) only one at a time is used, the 


other being moved to one side, to be available 
when the former burns out. The electron emission 


from the ion source passing down between the 
dees, strikes the iron (lid or base plate) and, 
inasmuch as no water cooling is provided at the 
point, large electron currents can give rise to 
excessive heating and evolution of gas. In the 
future when larger currents are used, it will be 
necessary to provide water cooling where the 
electrons strike, but for the currents used until 
now (less than 500 ma at 1000 volts) this has 
not been necessary. 

From the magnitude of the current of high 
speed ions arriving at the periphery in com- 
parison to the estimated amount of ionization at 
the ion source, it is clear that the focusing 
action of the electric and magnetic fields is such 
that a very considerable fraction of the ions 
starting at the center reach the target. Because 
space charge effects are negligible for such small 
currents at such high voltages, it is to be expected 
that the output of high speed ions would be pro- 
portional to the production of ions at the center, 
and therefore, it is of importance to consider 
means of increasing the source of ions. 

The ion current will be proportional to the 
electron current and therefore it is necessary 
only to increase the filament emission in order to 
increase the ion output. Unfortunately, large 
electron currents are limited by space charge, 
and it is necessary to increase the accelerating 
electrical fields in order to increase emissions. 
In the present arrangement, it is not efficacious 
to increase the electron emission by increasing 
the electron accelerating voltage, because the 
probability of ionization of the gas—for example, 
the production of protons in hydrogen—falls off 
with increasing energy of the electrons, as has 
been shown by Bleakney.® The optimum energy 
of the electrons for the formation of protons or 
deuterons is in the neighborhood of 200 volts, 
and accordingly it is desirable to design the 
source of electrons to yield large currents at low 
voltage. 

One way to accomplish this is to increase the 
size of the filament, since in practice it has been 
found that there is a relatively large central 
region that is suitable for the ion source. Over 
distances of as much as two inches the output 
of high speed ions in relation to the production 
at the ion source does not vary more than 50 


® W. Bleakney, Phys. Rev. 35, 1180-1186 (1930). 
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Fic. 4. Detail of accelerator electrodes and mountings. 


percent. Space charge effects can be reduced also 
by the use of an accelerating electrode close to 
the filament. One suitable arrangement would be 
that of a heavy (perhaps 50 mil) straight tungsten 
wire several inches long, extending parallel to the 
diametrical edges of the dees, midway between 
them. With such a straight wire filament held 
accurately in position, it would be feasible to 
place, very close by, a water-cooled accelerating 
electrode having a narrow rectangular slot 
which would produce much higher fields at the 
filament surface than are obtained in the present 
arrangement. 

Thus it appears that there should be no great 
difficulty in producing 5 or 10 times as much 
effective electron current as at present used, and 
accordingly that much greater output of high 
speed ions. At the present time, the maximum 
output current of high speed ions is about 25 
microamperes and the usual steady-running 
currents are about 20 microamperes. With the 
ion source increased, the output of high speed 
ions should be more than one hundred micro- 
amperes. 


The accelerator electrodes 


The semi-circular hollow accelerator electrodes 
or dees, are clearly shown in Fig. 2 and i# more 
detail in Figs. 3 and 4. They are made entirely of 
copper, copper sheeting being soldered to a 
heavier copper frame, with water-cooling pipes to 
take care of the heating due to the high frequency 


charging currents and the bombardment by the 
circulating ions. The dees have an over-all 
diameter of 24}’’, while the exit slit, as may be 
seen in the drawing of Fig. 4, is such that the 
diameter of the final circular path of the ions is 
23’’. The dees have a width at the center of 13” 
which uniformly decreases to 1}” at the periph- 
ery. The widening of the dees at the center 
favors the production of ions there without 
unduly increasing the capacity of the system. 
The dees are supported on the ends of 1 
copper rods 22” long, which are in turn sup- 
ported by means of Pyrex pipe insulators, the 
details of which are shown in the drawings of 
Figs. 2 and 4 and the photograph of Fig. 3. The 
insulators are standard }”’ wall Pyrex glass flange 
pipes of 3” inside diameter which have been 
blown out in the central region to 4” diameter 
(the insulators shown in Fig. 3 are not blown 
out to 4”), the glass being connected to the 
metal parts by fuse wire gaskets, according to 
the suggestions of D. L. Webster and P. A. Ross. 
From a mechanical standpoint this arrange- 
ment for supporting the dees is very satisfactory, 
but there have been troubles with cathode-ray 
punctures of the insulators and high frequency 
hysteresis heating. These troublesome matters, 
however, have been practically eliminated by 
providing the cathode-ray shields shown in ‘the 
drawings and also by providing air cooling of 
the glass. This geometry of the electrode arrage- 
ment not only practically eliminates cathode rays 


7 
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but also distributes the high frequency electric 
fields in such a way as to minimize hysteresis 
heating effects; but nevertheless there remains 
enough heating of the glass to make it necessary 
to provide jets of air on the glass surfaces to 
keep them cool. With air flowing over the glass 
insulators, we have had no trouble with insulator 
breakdown for several months, and presumably 
the insulators will last indefinitely. 

Of course, the endeavor always will be to 
produce larger numbers of high speed ions, at 
ever higher energies, and this will mean an ever- 
increasing applied high frequency voltage, so 
that doubtless the point will be reached where it 
will be difficult to prevent breakdown of the 
insulators. It will then be desirable to eliminate 
them practically altogether. This can be accom- 
plished by mounting the dees on heavy and rigid 
inductance coils in vacuum. The high frequency 
power for such a circuit could be introduced 
through insulators at relatively low voltage and 
thereby it would be possible to increase the 
power input many fold without insulator dif- 
ficulties. This arrangement, which makes use of 
vacuum insulation in the manner of the Sloan’ 
radio frequency high voltage generator, will 
doubtless be used in our laboratory sometime in 
the future. 

The question of the effect of dee width (i.e., 
height of opening) on the output of high speed 
ions is one of much interest and importance. The 
fact that many microamperes of high speed ions 
are obtainable indicates that the focusing action 
of the curved electrical fields between the dees is 
very good, but it is impossible to predict exactly 
the effect on the focusing action of the geometry 
of the accelerating electrodes. Of course, it can 
safely be argued that the number of ions pro- 
duced at the center will be at least a linear 
function of the dee width at the center, but 
whether the peripheral width of the dees affects 
the ion output very much can only be determined 
by experiment. 

To examine this question we have recently 
tried dees having a uniform internal width of 2” 
and have found that steady output currents of 
50uA are obtainable under approximately the 
same conditions that 20uA are obtainable from 
the narrower dees. Thus it seems that the output 


7D. H. Sloan, Phys. Rev. 47, 62-71 (1934). 


of ions increases somewhat faster with dee 
width than proportionally, and therefore a 
cyclotron designed for very large current output 
should not only have an ion source capable of 
giving large electron emissions, but also should 
have very wide dees. 


Withdrawing the ions 


The arrangement for withdrawing the ions is 
shown in Figs. 2 and 4. The ions emerging from 
one of the dees are deflected by a transverse 
electrostatic field and travel between the de- 
flector plate and the dee on an approximately 
cricular path of radius 30 percent greater than 
the final circular path within the dee (of radius 
113’). Since the magnetic force acting on the 
circulating ions is equal to the centrifugal force, 
which is equal to twice the energy of the ions 
divided by the radius of curvature, it follows 
that the magnetic force on 5 Mev ions traveling 
on a radius of 30 centimeters would correspond 
to an electrical field of about 330,000 volts per 
centimeter. Thus, an applied electrostatic field 
of the order of magnitude of 100,000 volts per 
centimeter is used in the present arrangement. 
Since the pole face diameter is 27}’’, the ions 
emerging from the interior of the dee are at a 
distance of 2}”’ from the wall of the vacuum 
chamber and with the deflecting fields used in the 
present experiments they are brought to the 
tank wall in a distance less than one-fourth of 
the circumference. As may be seen in Figs. 2 and 
4, a water-cooled target mounted on a ground 
joint may be rotated in and out of the beam. 
Usually two targets make up the assembly, one 
of tungsten and one of berryllium, the former 
being used mainly for test purposes (to measure 
the beam intensity with a minimum of neutron 
background) while the latter is used as a source 
of neutrons. The arrangement of the inset in the 
chamber wall adjacent to the target makes it 
possible to place objects near to and directly in 
front of the source of the neutrons. Since the 
neutrons partake of the momenta of the bom- 
barding deuterons, in general there is an enhance- 
ment of the intensity of the neutrons in the 
forward direction so that there is real advantage 
for investigations with neutrons to be able to 
place objects near and directly in front of the 
beryllium target. 
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It is frequently desirable to expose objects 
directly to the deuteron beam and for this pur- 
pose it is convenient to bring the beam out into 
the air. A thin platinum window is provided for 
this purpose. The platinum foil 0.0001” thick 
(stopping power, 1.5 centimeters of air) is 
soldered to a brass grid which is water cooled 
and able to dissipate the heat liberated by the 
deuteron beam striking and passing through the 
window into the air. The holes in the grid are 
such as to allow a little more than half of the 
beam to pass through. A beam of several micro- 
amperes of deuterons emerging into the air 
produces a luminosity of lavender color of con- 
siderable intensity and a photograph of a 5 
microampere beam of 5.8 Mev deuterons passing 
out into the air from the window to a distance of 
25 centimeters is shown in Fig. 5. 

Although the arrangement here shown is quite 
convenient and useful in many investigations, 
for some purposes it would be desirable to with- 
draw a beam of ions through a vacuum tube to a 
point which is at a considerable distance from 
the accelerator system, in order to remove 
background radiations, produced by ions cir- 
culating between the dees that strike parts of 
the apparatus other than the intended target. 
This can be readily accomplished by attaching 
a vacuum tube of any desired length in place of 
the platinum window, and putting the platinum 
window on the other end of the tube. If no 
attention is given to refocusing, the ion beam 
will progressively diverge as it passes through 
such an extended vacuum tube, partly on 
account of the inhomogeneity of the magnetic 
fields external to the vacuum chamber. However, 
it is not a difficult matter to correct this inhomo- 
geneity; in fact, pieces of iron can be placed 
around the tube in such a way as to refocus the 
beam of ions at the extended exit window. This 
has not actually been done as yet in our labor- 
atory, but we intend to do so very soon for the 
purposes of certain experimental investigations 
in progress. 


Oscillators 


The source of high frequency power is a self- 
exciting (tuned grid—tuned plate circuit) oscil- 
lator which has the advantage of simplicity 
involving a minimum of equipment. Two oscil- 
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lator tubes of our construction which have been 
recently described,* are used in a push-pull 
circuit. These tubes are rated at 30 kilowatts, 
although in practice much smaller powers are 
used. The anodes are supplied with three phase 
full wave rectified voltage up to a voltage of 
about 10,000 volts. The rectifier consists of 
three 220 volt to 6600 volt transformers con- 
nected for three phase full wave rectification with 
872A rectifier tubes. In order to draw as large 
currents as might be desired, 12 of these rectifier 
tubes are actually used, parallel groups of two 
being connected with suitable stabilizing resist- 
ances in series. When the cyclotron is running, 
delivering 20 microamperes of 5.5 Mev deuterons, 
the anode voltage on the oscillators is about 
8000 with an anode current of about three 
amperes. Of the total input power of about 25 
kilowatts, it is found that about 12 kilowatts of 
high frequency power is delivered to the ‘‘tank 
circuit.”’ 

Although the simplicity, of the self-exciting 
oscillator commends itself for many practical 
purposes, it has the disadvantage of inconstancy 
of frequency and a certain instability in opera- 
tion; for variable frequency means variable 
output of high speed ions and for many purposes 
it is desired to have as great constancy as possible 
in the energy and intensity of the beam. We are 
at the present time engaged in the construction 
of a new oscillator system, using a master oscil- 
lator and successive stages of power amplifiers, 
an arrangement which should deliver to the 
accelerating system radio frequency power of 
constant frequency and great stability. With the 
present arrangement, it is possible to keep the 
beam steady within 10 percent by adjusting the 
magnetic field from time to time. The new 
oscillator should make it possible for the cyclo- 
tron to deliver a steady and homogeneous beam 
of high speed ions for considerable periods of 
time, perhaps hours, without requiring attention. 


ADJUSTMENTS 


The dees 


The sylphon arrangement in the dee supports 
makes it possible to change the position of the ' 








8D. H. Sloan, R. L. Thornton and F. A. Jenkins, 
Rev. Sci. Inst. 6, 675-82 (1935). 
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Fic. 5. Photograph of 5uA of 5.8 Mev deuterons (range 
25 cm) emerging from platinum window of cyclotron. 


dees during operation. As would be expected, it 
is often found that the optimum position is that 
in which the dees are accurately parallel to each 
other and symmetrically placed with respect to 
the magnetic field. Deviations from this sym- 
metrical adjustment, of one millimeter for 
example, are sometimes detectable by as much as 
a 10 percent reduction in the output of ion 
current. This fact emphasizes the importance of 
constructing the dees so that they do not tend 
to wobble or sag out of position, and also the 
desirability of the sylphon arrangement for 
making adjustments. In this latter regard, it is 
found that a slight lack of parallelism of the 
dees can be partly compensated by adjustments 
of their relative position. The sylphons make it 
possible also to vary the separation of the dees, 
which when parallel are about ?” apart. No 
great variations are produced by varying the 
separation, but it does seem in general that the 
further the dees are apart the less the output of 
high speed ions. It is not known whether this is 
due to change in the production of ions at the 
center or to the focusing action of the fields out 
toward the periphery. 

In the light of the above remarks, it would 
seem that with dees and supports of sufficient 
rigidity, the sylphon on the supports would not 
be needed, as the dees could be placed accurately 
in the optimum symmetrical positions at the 
time of assembly and would keep these positions 
during the course of operation. It is difficult to 
say whether in practice the simplification of the 


dee supports obtained by eliminating the 
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svlphons would more than compensate the loss 
of the possibility of making dee adjustments. 


The magnetic field 


Neglecting relativity corrections, which are 
unimportant, the simple theory of operation of 
the cyclotron would suggest that the ideal mag- 
netic field is uniform over the whole accelerating 
system. In practice, this ideal is not achieved, 
partly because of the imperfections of the 
geometry and the inhomogeneity of the iron of 
the magnetic circuit, and partly because the 
saturation of the iron pole pieces results in a 
decrease of the magnetic field toward the periph- 
ery of the dees; this latter, surprising as it 
may seem, is an inhomogeneity that is really 
desirable. The actual variation of the magnetic 
field from the center outward in the present 
cyclotron is shown in Fig. 6, where it is seen that 
from the center to the edge of the magnetic pole 
faces, the field falls off 25 percent, while even 
from the center to the point where the ions are 
withdrawn from the accelerating electrodes, 
radius 11}’’, the field actually diminishes by 7 
percent. At first sight, one would think that such 
a marked radial inhomogeneity of the magnetic 
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Fic. 6. The radial distribution of the magnetic field 
Che shaded areas indicate the magnet pole faces and the 
iron plates of the cyclotron chamber. 
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Fic. 7. Photographs of recoil protons in Wilson cloud chamber filled with hydrogen at pressure of 100 cm Hg produced 


by neutrons from Be target bombarded with 7uA of 5.5 Mev deuterons. Distance of cloud chamber from Be target, 
; 5s M : 


6 ft., (b) 40 ft. 


field would make it impossible to accelerate the 
ions to the periphery, but on second thought it 
is realized that if the high frequency accelerating 
voltage is of the order of magnitude of 100,000 
volts in the production of 5 Mev deuterons, the 
ions circulate only about 25 complete revolutions 
and under these circumstances an inhomogeneity 
toward the periphery of only 7 percent would not 
seriously affect the synchronization of the cir- 
culating ions with the oscillating field. The high 
frequency voltage used in the present cylotron 
has not been measured reliably, but it is esti- 
mated to be somewhere between 50,000 and 
100,000 volts. 

This radial gradient of the magnetic field, 
which becomes quite large near the periphery, 
has the effect of opening up the spiral paths. 
Instead ofrtraveling on;cifclés within the dees, 
the decreasing field toward the edge causes the 
ions that do not have their path centers exactly 
at the center of the field> to spiral around on 
continuously increasing radii of curvature. The 
extent of this effect of the field gradient for the 
present cyclotron is illustrated by a simple 
approximate calculation of the path of an ion, 


a 


on the assumption that the center of the spiral 
for the revolution, just before the ion is with- 
drawn from the dee, is }” off center, and that the 
high frequency accelerating voltage is 50,000. 
The path of the ion, as approximately deter- 
mined in this way, makes its exit from the dee at 
a distance of more than }” from the previous 
spiral path. Thus it is seen that this outward 
field gradient produces a very considerable 
opening up of the spirals. This is obviously a 
desirable characteristic, inasmuch as it facilitates 
the withdrawal of the ions from the accelerating 
system. 

But the fact that a shift of the center of the cir- 
cuit of thespirals by }” produces such a large effect 
on the spiral paths imposes a rigid requirement 
that the magnetic field be highly homogeneous in 
azimuth. If the average value of the magnetic field 
in one quadrant were appreciably different from 
that in the others, the effect would be to produce 
a migration of the centers of the spiral paths, and 
in turn, serious distortion of the paths, greatly 
magnified by the spiraling action of the radial 
field gradient. In adjusting the magnetic field of 
the cyclotron, therefore, it is necessary to insert 
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pieces of iron between the pole faces and the 
bottom and top of the cyclotron chamber (as 
indicated in Fig. 6, }’’ gaps are provided for this 
purpose) in a purely empirical manner to correct 
the azimuthal inhomogeneity. In general, it may 
be said that increasing the magnetic field, let us 
say on the southern portion of the pole faces, 
has the result of causing the centers of the spiral 
paths to migrate toward the east, for clockwise 
rotation as seen from above. If it is found, 
therefore, that the ion beam emerging from the 
exit slit of the dee has an energy less than that 
to be expected, the presumption is that the 
centers of the spiral paths have migrated toward 
the exit slits so that the final radii of curvature 
are less than the distance from the exit slit to 
the center of the apparatus. Under such circum- 
stances, the energy of the particles could be 
increased by inserting shims to increase the 
magnetic field in a particular region, which 
would cause the spiral paths to migrate back 
toward the center, and hence increase the radii 
of the final paths. Of course, it should be em- 
phasized that in practice the adjustment of the 
magnetic field by the insertion of shims is largely 
an empirical matter. Although the effects of the 
shims cannot always be understood, the shim 
adjustment is fortunately not difficult to make; 
inasmuch as without any shims at all, a beam 
of ions is obtainable. Thus, the homogeneity, 
intensity and energy of the ion beam can be 
observed continuously as the magnetic field 
adjustment is carried out. 
PERFORMANCE 

Hydrogen ions 

The present cyclotron delivers several micro- 
amperes of 5 or 6 million volt deuterons without 
any adjustments after the vacuum chamber has 
been placed between the poles of the magnet. In 
the beginning of operation, the circulating ion 
beam liberates gas from the various metal parts 
of the apparatus, so that during the first few 
hours of operation, the output beam of ions 
steadily increases as the gas liberation decreases. 
It is found that for satisfastory operation, the 
pressure of foreign gases must be kept below 10-5 
mm Hg, and that the largest currents are ob- 
tained when the hydrogen pressure is of the 
order of magnitude of 10-* mm Hg. By making 


the various adjustments (i.e., the deflecting 
potential and the magnetic field distribution with 
shims and the position of the dees, the deflecting 
plate, and the ion source filament, which usually 
is at the center) the ion currents are usually 
built up to optimum values of 20 to 25 micro- 
amperes of 5.8 Mev deuterons, when using dees 
having an internal peripheral width of 1}”, or 
about 50 microamperes of 4.3 Mev ions when us- 
ing 2’’ width dees of smaller radius. It is found 
that the obtainable current is not appreciably 
dependent upon the output energy of the ions, as 
we have, when using the narrower dees, obtained 
just as large currents above 5 Mev as below 5 
Mev. The maximum energy of deuterons that 
has been produced by the present cyclotron is 
6.3 Mev. 

The obtainable beams of high speed protons 
are usually equal in current to the deuterons, but 
less in energy on account of the required higher 
frequency of the electric oscillations. 


Helium ions 


The cyclotron, when adjusted for deuterons, is 
also in adjustment for doubly charged helium 
ions. The helium ions having twice the charge 
and mass relative to deuterons, receive twice the 
increments in energy each time they pass between 
the dees and arrive at the periphery with twice 
as much energy. Thus, apart from difficulties of 
producing doubly charged helium ions at the 
center, the cyclotron produces 10 Mev alpha- 
particles as readily as 5 Mev deuterons. 

On one occasion, we replaced the deuterium 
with helium in the cyclotron chamber and 
obtained 1/10uA of 11 Mev alpha-particles. 
(Several times larger currents doubtless could 
have been obtained by completing the adjust- 
ments.) That the beam of high speed ions was 
really made up of such high speed helium ions 
rather than deuterons of lower energy was 
definitely established in several ways. Thus, for 
example, it was observed that the alpha-particle 
beam attained its maximum intensity at a mag- 
netic field slightly lower than that for the 
deuteron beam, had a range one-half that of the 
deuteron beam under similar circumstances, and 
produced a luminosity in the air in intensity 
equivalent to that of 4 times as large a deuteron 
current. 
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As yet we have not carried out experiments 
with these energetic alpha-particles, excepting 
to observe that they do not give rise to as large 
a neutron emission from a bombarded beryllium 
target as deuterons of half their energy. 


Neutrons 


Among the many elements that emit neutrons 
under deuteron bombardment, beryllium is par- 
ticularly noteworthy. The neutron emission from 
a beryllium target bombarded by 7uA of 5 Mev 
deuterons is indicated by the cloud chamber 
photographs of Fig. 7. With the hydrogen-filled 
cloud chamber at a distance of 6 feet from the 
target, the number of recoil protons is too great 
for counting purposes, and a more suitable 
number is obtained at 40 feet. Although it is 
difficult to compare vastly different numbers of 
neutrons, estimates made in several ways indicate 
that the emission of neutrons from the beryllium 
target of the cyclotron is more than 10° times 
that from a mixture of one curie of radon and 
beryllium. 

This great neutron intensity facilitates many 
nuclear studies, but on the other hand, hampers 
some investigations. Among the practical dif- 
ficulties, there are the biological effects produced 
by the neutron rays, which make it necessary 
to provide considerable protection for the in- 
vestigator. Some studies*""' have 
already been carried out which have shown that, 
in biological action, the usual neutron emission 


biological 


® J. H. Lawrence and E. O. Lawrence, Proc. Nat. Acad. 
22 (2), 124-1933 (1936). 

10R, E. Zirkle and P. A. Aebersold, Nat. Acad. Sci. 
22 (2), 1934-138 (1936). 

1 J. H. Lawrence, P. A. Aebersold and E. O. Lawrence, 
Proc. Nat. Acad. 22, 543-557 (1936). 
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from the beryllium target of the present appa- 
ratus is equivalent to the gamma-radiation from 
100 grams of radium. Accordingly, at the present 
time, for the protection of the operator, the 
cyclotron is controlled at a distance of 40 feet 
from the target with suitable intervening ab- 
sorbing material. As the intensity of the neutron 
source is further increased, it will be necessary to 
provide additional absorbing screens. 


Artificial radioactive substances 


With 5 Mev deuterons, it has been found 
possible” to produce radioactive isotopes of 
many of the elements throughout the periodic 
table. In many cases, the yields of the radio- 
active substances are quite large; as for example, 
a day’s bombardment of sodium metal with 20 
microamperes of 5 Mev deuterons produces 
more than 200 milligrams—equivalent of radio 
sodium," i.e., an amount of radio sodium having 
a gamma-ray activity equivalent to that of 200 
milligrams of radium. That such large amounts 
of radioactive forms of many of the elements can 
be manufactured in the laboratory is of much 
importance in opening up new avenues of research 
both in the physical and the biological sciences. 


We are greatly indebted to our laboratory 
colleagues, who have contributed so much to 
the work here reported. We acknowledge with 
thanks also the continued support of the Re- 
search Corporation, the Chemical Foundation, 
and the Josiah Macy, Jr., Foundation. 

127. M. Cork and E. O. Lawrence, Phys. Rev. 49, 788 


(1936). J. J. Livingood, Phys. Rev. 50, 425 (1936). 
13 E, O. Lawrence, Phys. Rev. 47, 17. (1935). 
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Radioactivity Induced in Sulphur 


R. SAGANE, Radiation Laboratory, Department of Physics, University of California, Berkeley, California 
(Received October 5, 1936) 


Under bombardment of 4-6 Mev deuterons, sulphur becomes radioactive with the emission 
of both positive and negative electrons with 3 periods 3+0.1 m, 3341 mand 14+0.3 d. These 
decay periods were found to have their origin in P*®, Cl* (both positron emitters) and P® 
(negative electron emitter). The energy distribution of positrons from Cl was studied by a 
cloud chamber. When bombarded by neutrons two radioactive substances were produced 
which emit negative electrons. The half-value periods were found to be 2.6+0.2 h and 14 d, 


which are characteristic of Si*! and P®. 


INTRODUCTION 


INCE the discovery by Curie-Joliot of arti- 

ficial radioactivity, many radioactive iso- 
topes have been found produced by bombard- 
ment with protons, deuterons, a-particles and 
neutrons. Many light elements have been investi- 
gated in this laboratory by the use of the deuteron 
beam from the cyclotron. The work described in 
this paper is an investigation of radioactivity 
induced in sulphur by both deuteron and neutron 
bombardment. 


APPARATUS 


The source of deuterons was the large cyclo- 
tron' in the Radiation Laboratory. At the time of 
the bombardment it was adjusted to furnish 
deuterons of energy of 4-6 Mev. 

The targets were, in all cases except the first 
two exposures, sulphur fused on Pb plates and 
covered with two sheets of Pt foils, 4 mm in air 
equivalence, to avoid contaminations. In the first 
two cases, sulphur was fused on Al plates and 
covered with Al foils. The bombardments were 
made in vacuum on a special target holder with a 
cooling device to prevent the sample from 
evaporating. 

All the measurements of activity were carried 
out on a Lauritsen type quartz fiber electroscope 
with an Al window and absorption measurements 
were made by placing Al sheets between the 
source and the window. A Wilson chamber? was 
used to investigate the nature of the ionizing 
particles from the irradiated samples and also for 


1E. O. Lawrence and M. S. Livingston, Phys. Rev. 45, 
608 (1934). 

21 am indebted to Dr. Kurie for the use of the cloud 
chamber. 


the measurement of the energy distribution of the 
positrons from radiochloriné. 


I. ACTIVATION BY DEUTERONS 
Decay periods 


A thick sulphur target on an Al plate was bom- 
barded in vacuum at first for about 2 hours by an 
ion beam of 4 Mev in energy. At the end of the 
exposure the target was placed at the window of 
the electroscope. The measured activity as a 
function of time was plotted on a logarithmic 
scale and analysis of this curve shows the existence 
of two radioactive substances. The half-periods 
were observed to be 3341 m and 14+0.3 d, and 
also there were traces of 2.6 # and 14 h ac- 
tivity. The ratio of the initial intensities of 
these substances was estimated roughly to be 
360 : 1.6: 0.3 : 0.2, respectively. Repeating the 
same exposure, the same type of decay curve was 
reproduced. 

To find out whether there is a shorter period 
activity, a sulphur target on Pb covered with 
thin Pt foil, 0.0001” in thickness, was bombarded 
for 10 minutes. After due allowance for the 33 m 
activity, a half-life of 340.1 m was found. A 
short exposure of 5 m was made afterward and 
the same result was reproduced. 

To investigate the nature of the ionizing 
particles, they were observed visually in a Wilson 
chamber which was placed in a magnetic field of a 
few hundred gauss. Both positrons and negative 
electrons were emitted from the target. 


The active substances 
Ordinary sulphur contains three isotopes, S*, 


S* and S*, the relative abundance of which is 
reported to be 96:1: 3. If sulphur combines 
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Fig. 1A. Logarithmic plot of the activity of the chemi- 
cally separated P*°, Dots show the measured curve. x’s 
show the corrected curve for a small amount of contami- 
nation of Cl of about 1 percent in intensity. 


with deuterons, Cl*, Cl, Cl8* are expected to 
result with the capture of deuterons (this is 
supposed to be very improbable); P**, P*!, P% 
with the emission of a-particles; Cl, Cl§*, Cl 
with the emission of neutrons; and S*, S*, S* 
with the emission of protons. 

The identity of P*° and P® is apparent since the 
periods characteristic of these substances are 
known with good accuracy. As for Cl*, Cl** half- 
value periods of 40+5 m* is given for Cl* and 50‘ 
to 35 minutes for Cl** (probably 37 m®). The 
accuracy of these measurements is not good 
enough to identify the 33 minutes activity. 

Chemical tests were also made to check their 
identity. The sulphur sample was dissolved in 
CSs, and mixed with dilute HNQOs, in which a 
small quantity of red phosphorus and a drop of 
dilute HCl was added. Then the solution was 
heated for a short time until the phosphorus was 
completely dissolved. After filtration of the water 
solution, Cl ions were precipitated with AgNO; 


30. R. Frisch, Nature 22, 434 (1934). 

‘A. T. Alichanow, A. J. Alichanian and B. S. Dzelepow, 
Nature 135, 393 (1935); E. Fermi and others, Proc. Roy. 
Soc. A146, 483 (1934); F. N. D. Kurie, J. R. Richardson 
and H. C. Paxton, Phys. Rev. 47, 796 (1935). 

5A. R. Olson, W. F. Libby, F. A. Long, R. S. Halford, 


J. Am. Chem. Soc. 58, 1313 (1936). 
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Fic. 1B. Logarithmic plot of the activity of the 
chemically separated P®. 


then PO, ions were precipitated with magnesia 
mixture. 

The logarithmic plot of the decay curves of 
these chemically separated samples evidenced 
excellent straight lines as shown in Fig. 1; con- 
sequently half-lives of P®° (by short exposure), 
P® and Cl (by long exposure) were determined 
with sufficient accuracy. It is therefore obvious 
that the 3 m activity is due to P**, 14 d to P® and 
33 m to Cl. 

Examining the chemically separated Cl in the 
Wilson chamber it was found that the sample 
emitted chiefly positrons and a small number of 
negative electrons. The latter can be explained as 
Compton recoil electrons of annihilation radia- 
tion of positrons and also a small amount of 
contamination of P*. A study of the periodic 
table reveals that 33 m activity is probably due 
to Cl* or perhaps Cl*, because the intensity of 
this Cl seems to be too strong to explain by the 
following process, if the fact is taken into ac- 
count that the relative abundance of S* is only 1 
percent : S¥+ D?=Cl*+-2n. 

Thus the reactions are supposed to be 


S®+ DP? = P*°+ Het P30_,S}39 4+ et 

S84 DPD? = Cl*+ (S$? + D? = Cl+ 7! 
S? + ] )? — C]*4) 

C]/#—S*4 +et 


S44 >.>? = Pp? + Het P®_—,S%4 + ¢ : 














RADIOACTIVITY INDUCED IN SULPHUR 1143 


ACTIVITY 





! 2 3 4 
HOURS 


Fic. 1C. Logarithmic plot of the activity of the 
chemically separated radiochlorine. 


P® can also be expected to result from a neu- 
tron bombardment on sulphur and there is a 
probability that P® might be created by neu- 
trons from D+D. An experiment to test this 
point was made. Two similar samples of sulphur 
were put one upon another on the target holder, 
and the upper target was bombarded by deu- 
terons. As these targets were rather thin, in- 
tensities on both targets of the neutrons created 
on the upper target by D+D were supposed to 
be the same order of magnitude. After bombard- 
ment the intensities of the activity of the lower 
target was found to be less than 1/200 of the 
upper target. This definitely demonstrated that 
P® on the upper target was chiefly created by 
deuteron impacts. 


The radioactive radiations 


The beta- and gamma-rays from P*® and P® 
have been studied in detail by many others and 
there is no need to reiterate it here. However, 


Frisch is the only person who has succeeded in 
making Cl**. He created it by bombardment of 
P3! by a-particles, consequently the intensity 
must be very weak and the accuracy of the 
measurement is rather poor. In the present ex- 
periment radioactive chlorine was made by 


deuterons with strong intensity, the initial in- 
tensity of 2 hours’ exposure being 300 divisions 
per second, that is 3X10‘ times the natural leak. 
Some detailed investigation of the radiations 
from radioactive chlorine was made. 

Gamma-rays. Rather weak gamma-rays were 
found, the absorption coefficient of which was de- 
termined with lead sheets. A value of 5 mm of 
lead for the half-value thickness was obtained. 
Although the geometrical conditions were not 
satisfactory, this value probably corresponds to 
an energy of about 0.5 Mev. This is in agreement 
with the energy of the gamma-rays which is 
expected when positrons from Cl are annihilated. 

Beta-rays. The energy distribution of 8-rays 
from radioactive chlorine was studied by the use 
of the cloud chamber with a magnetic field of 
about 400 gauss. The source was chemically 
separated from the irradiated sulphur in half an 
hour after exposure and then mounted inside of 
the chamber, covered with a thin Al foil to pro- 
tect the chamber from contamination. This 
source was rather thick owing to the addition of 
too much chlorine in the course of chemical 
separation and one could not get rid of some 
irregularity in absorption; still the result shows a 
typical positron distribution and the Konopinski- 
Uhlenbeck® plot gives a fairly good straight line 
as shown in Fig. 3. 

The upper limit is estimated to be 3 Mev. The 
point corresponding to 3 Mev and 33 m placed on 
Sargent’s’ curves of beta-ray periods and energies 
lies right on the second group line. 


II. AcTIVATION BY NEUTRONS 


A block of sulphur (20153 mm) was put 
about 5 cm from the Be target of the cyclotron 
and irradiated for about 2 hours by neutrons. 
After irradiation the activity was measured as 
before. The decay curve was analyzed in two 
exponential curves showing the existence of two 
radioactive substances, the decay periods of 
which are 2.6 hours and 14 days. 

These periods are characteristic of Si*' and P® 
which are expected by the following reactions. 


®E. J. Konopinski and G. E, Uhlenbeck, Phys. Rev. 
48, 107 (1935); F. N. D. Kurie, J. R. Richardson and 
H. C. Paxton, Phys. Rev. 49, 368 (1936). 

7B. W. Sargent, Proc. Roy. Soc. 139, 659 (1933). 
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Fic. 2. Distribution histogram for the positrons emitted 
from chemically separated radiochlorine. 


S*+n =Si*!+ He‘, 
S24 =P2®+H), 


Rough chemical tests were made by dissolving 
activated samples with CS». A small quantity of 
Si powder and red phosphorus was added, as a 
mixture in dilute HNO;. Then it was heated 
slightly till the powder of phosphorus was dis- 
solved completely by HNO. After filtration of 
HNO; solution, phosphorus was precipitated by 
the magnesia mixture. The CS, solution was also 
filtered afterward and the Si powder was obtained 
as a residue. 

By examining these samples it was found that 
the 14-day activity was due to a phosphorus 
isotope and the 2.6 hour activity to an Si isotope. 

As for the CS: solution, it was dried to powder 
of sulphur. This was then oxidized by sodium 
chlorate and HNO; to H2SO,. The latter was 
precipitated with BaCl, as BaSO, in HCI solu- 
tion. But this sample showed practically no 
activity. 

A search for S* 

According to the consideration of reactions 
with S, a radioactive sulphur S*, which is said to 
have a long half-life of about 80 days,® is expected 
to result from both deuteron and neutron 
bombardment. 

A search for this isotope was made by making a 


SE. B. Andersen, Zeits. f. physik. Chemie B32, 237 


(1936). 


careful chemical separation of sulphur as before. 
The first two or three chemically separated 
samples showed some small activity, but by 
repeating the process more carefully the activity 
was found to be much smaller than before and by 
measurement of the half-life, this activity was 
found to be due to the contamination of P®. 

The result shows that the deuteron bombard- 
ment for 3 hours with 13 microamperes beam on 
the target, cannot make S* of intensity of more 
than 1/1000 of P®. 

As for neutron bombardment, the longest ex- 
posure was made on one sample which was put 
about 10 cm from the Be target of the cyclotron 
for about 20 days completely covered with 
paraffin. During these days the cyclotron ran 
very steadily with about 10-15 microamperes of 
5.5 Mev deuterons and nearly every day neutron 
bombardment on mice and wheat seeds, etc., was 
continued. Yet the activity of the chemically 
separated sulphur from this sample was about the 
same order of magnitude of the natural leak and 
even this activity was also found afterward to be 
due chiefly to P® contamination. 

In case of this thick target, it is quite difficult 
to estimate the lower limit, but it seems the cross 

















Fic. 3. K—U plot for radiochlorine. This plot extrapolates 
to an upper limit at E+/=6.9, about 3 Mev. 
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section for neutron capture by sulphur must be ing the cyclotron and for his helpful discussion 
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A Positive Energy Wave Packet Solution of the Dirac Equation 


CarL Eckart AnD B. L. GrirFinG, Ryerson Laboratory, University of Chicago, Chicago, Illinois 
(Received October 3, 1936) 


A wave packet solution of the Dirac equation is obtained by superposition of positive energy 
solutions only. The solution involves an arbitrary parameter, and it is shown that this can be 
chosen so that all three dimensions of the packet are less than h/mc. The possibility of con- 
structing such a packet had been doubted. It is also shown that in the limit Ax=0, AxAp, 
= 1.49h. It is not known whether this is the smallest value for the uncertainty product, subject 
to the condition that only positive energy solutions are used in the construction of the packet, 
or whether some other solution might yield a still smaller value. 


A QUESTION of some technical interest is 
the following: can a wave packet be con- 
structed by superposing only positive energy 
solutions of the Dirac equation for the free elec- 
tron, such that all its dimensions are smaller than 
h/mc? Apparently no conclusive answer to this 
question has been published.' In this note it will 
be shown that the omission of the negative energy 
solutions does not impose a lower limit on the 
size of the wave packet. This will be done by 
constructing a particular wave packet whose 
dimension at ¢=0 is arbitrary, and using only 
positive energy solutions in its construction. 
Because of the relativistic invariance of the 
theory, it is sufficient to consider only packets 
whose centroids are at rest in the system of co- 
ordinates chosen. It is known that in general, 
that packet for which AxAp, is smallest is ob- 
tained by using a Maxwellian distribution of 
momenta. In a preliminary investigation, it was 
found very difficult to discover the analog of 
this theorem when the positive energy condition 


1 It is known that several investigators have arrived at 
conflicting conclusions regarding this question. The only 
definite expression of opinion that has been published 
seems to be that of L. de Broglie, L’ Electron Magnétique, 
p. 286. 


is imposed. It was therefore decided to use a 
distribution such that the Maxwellian is ap- 
proached for large packets. A further considera- 
tion was that all necessary integrals should be 
calculable without dubious approximations. 

In the following, #/mc will be used as unit of 
length, and h/mc’ as unit of time. Then the Dirac 
equation can be transformed to the Schrédinger- 
Gordon equation by the substitution 


¥=1((0F /dt)+a-VF]—anF. (1) 
In order that y satisfy the Dirac equation, the 
“spinor potential’ / must satisfy 
(PF /d?)=V?F—F. (2) 
Wave packet solutions of this equation are easily 
written down in the form 
F,= SG,(p) exp [i(p-x+wt) |dp, (3) 
where ¢@ is the spin variable, and 
w=(1+")!>0. (4) 


If the negative sign of w is chosen in Eq. (3), 
then y will involve only positive energy solutions. 
For the present purpose, it is sufficient to take 


Gi=A exp (— fw), G.=G;=G,=0, (5) 
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where @ is an arbitrary constant determining the 
size of the packet, and A is the normalizing 
factor. Substituting into Eq. (3), transforming 
to polar coordinates in both configuration and 
momentum space, and using an obvious device 
to replace sin by exp: 

F,=(29A ir) | exp [ipr—w(8+it) |pdp. (6) 

“—o 

In order to obtain y from Eq. (1), the following 
functions are also needed : 


F,=(24rA/r) Sf exp [--- |pwdp, (6.1) 


F,=(2rA)S exp[--- ]p'dp. (6.2) 
In terms of these, the probability density is 
D= Ze Vo ° 
=|F,—F,|\?+|F,-F,/?/r (7) 


and is a function of r only. 

The solution thus obtained will be discussed 
only for =0, at which time the packet has its 
minimum extension. The integrals of Eqs. (6) 
can then be evaluated in terms of Hankel func- 
tions; setting 


gé=(r°+")}, sin 6=r/§, 


p=sinh a, 


it follows on substitution that 
arc 

F,=(2A ir) | sinh (2a) 
@ 


Xexp [—£ cosh (a—70) \da. (8) 


The path of integration can be displaced from 
the real axis to the line Imag(a) =@; then setting 
A\=a—10, 

Te 


F,=(sA, ir)| sinh (2+2i6) 


Xexp [—£ cosh dA ]dd 
=(7rA/r) sin (26) f cosh (2)) 


Xexp [—£ cosh \ ]dd. 


Since? 


? Jahnke-Emde, Tables of Functions, second edition, p. 
218. The formula given there is for the Hankel function of 
the second kind, but it is readily transformed to the one 
above when » is an integer and £ real. 
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f cosh (md) exp [ —& cosh \ ]dd 


x 


= (ri) 17, (it), (9) 
it follows that 
F, = (2A /r)(2i*)H2 (7) sin 28. (10) 
Similarly 


F,= (2A /2r)} (wit) Fs (2E) sin 36 
+ (rt?)H,"(7é) sin 6}, (10.1) 


F,= (9A /2) | (at*)H3 (2é) cos 30 
— (r?7)H, (2&) cos 6}. (10.2) 


These formulae do not furnish a very clear 
picture of the solution, but on plotting the prob- 
ability density as a function of r by the aid of 
tables,’ its graph is found to be very similar to a 
Gauss error curve. The half-width is a mono- 
tonically increasing function of 8 and is zero for 
8=0. 

As this is the principal point at issue, it was 
investigated in greater detail, using the known 
approximation formulae for the Hankel func- 
tions. When ¢ is small (i.e., 71 and B<1),* 


(wt"*") HT, (2£) = (m—1)1(2/E)"-++, (n 0), (11) 
which leads to F;=8xAB £&---, (12) 
F,=89A(362—r?)/#-++, (12.1) 


F,,=82AB(B?—3r*) /f--+, (12.2) 
whence 


D =647n°A*[ (1+6)*2+8(1—8)62]/E"---. (13) 


It is easily seen that this function has a maxi- 
mum for r=0; which becomes increasingly sharp 
as 8 approaches zero. The function diminishes 
monotonically with increasing 7 and has a single 
point of inflection, finally approaching zero as 
r-*. While this cannot be the true asymptotic 
behavior of D, yet the error introduced by using 
Eq. (13) for all values of r will vanish with 8. 
This can be seen in various ways: the correct 
asymptotic behavior of D can be obtained from' 


(wi"*") HH, (¢&) = (29 /é)' exp (—£), (ED 1). (14) 


3 Reference 2, p. 286. 
* Reference 2, p. 202. 
5 Reference 2, p. 204. 





POSITIVE ENERGY 
Thus the true curve will approach zero more 
rapidly than the approximate curve; because of 
the large negative exponent in Eq. (13), the 
resulting error will not be great, however. More 
precisely, when r&6<1, then D=D,)=647°A? 
(9 88); if on the other hand, B<r<1, then 
D=641°A?/r5<D,>. Thus most of the area under 
the D curve comes within the range of validity 
of Eq. (13). The validity of this argument can 
be checked, since the normalizing integral can be 
exactly evaluated by transforming to momentum 
space, and is given by a long, though finite, series 
of Hankel functions with the argument 2278. 
Using Eq. (11) to simplify this expression, we 
found that the normalizing integral is approxi- 
mately 
82r!A*[68->§—38-4--- ]. 


These first two terms of the series are identical 
with those of the series obtained by integrating 
Eq. (13). This result thus fixes the normalizing 
factor as A?=65/487*---. 

The evaluation of the mean-square radius Ar 
was also attempted by transforming to momen- 
tum space, but the integrals encountered were 
not very tractable. It was therefore decided to 
use Eq. (13) in obtaining an approximate value 
of the integral 


Ar? =8r | r*Dadr. 


“oOo 


It will be clear from the above discussion that 
the error involved in this procedure cannot be 
great provided that 8 is much less than unity. 
The result of an elementary integration is that 


SOLUTION OF 
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DIRAC EQUATION 
This verifies the result previously obtained by 
graphical methods, that Ar approaches zero with 
8. Since Ar’=3Ax* it follows that Ax=28/3 in 
the present units, or 


Ax = (2/3)Bh/mce, B<1. (15) 


While this would seem to be a conclusive 
answer to the question proposed above, it is of 
interest to note the value of the uncertainty 
product for this solution. The value of Ap, is 
exactly calculable by transformation to momen- 
tum space, and is again a finite series of Hankel 
functions. For small values of 8 it is approxi- 
mately given by 


Ap. = (5)'mc/B. (16) 


Hence AxAp,=1.49h in this limit; as has been 
remarked above, it is not known whether this 
is the smallest value that the uncertainty product 
can have when the packet is constructed only of 
positive energy solutions and Ax<<h/mc. It seems 
very likely that this minimum is a complicated 
function of Ax. 

It remains to discuss the case B>1. This has 
not been as carefully studied as the other limit, 
but the following considerations are not apt to 
be seriously in error. It is easily seen that the 
momentum distribution function is 


(27)*A[(w—1)?+p* ] exp (—28w). 


For large values of 8, this may effectively be 
approximated by the Maxwellian distribution 


A’ exp (—8p*). 


From this it may be concluded that when 6>1, 
AxAp,=h and Ax=8'mce. 
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Pressure Broadening of the Potassium Resonance Lines by Argon and Nitrogen* 


Gorpon F. Hutt, JR., Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
(Received October 3, 1936) 


Broadening of the potassium resonance lines in absorp- 
tion by argon and nitrogen pressures from 1 to 40 atmos- 
pheres has been observed. Half-widths and shifts of the 
lines, when plotted against “relative density,” yield curves 
that are straight lines up to about relative density 14, 
indicating that in this pressure region Lorentz type of 
broadening predominates. At about relative density 14, 
the curves depart from linearity. The results are in accord 


with theory, and partly in agreement with previous work. 
Theoretical contours for pressure broadened spectral lines 
have been fitted to those of experiment. Agreement is quite 
satisfactory, considering the approximations made in evalu- 
ating the theoretical curves. Three independent methods 
of calculation yield for the constant 6 (measure of van der 
Waals forces) the values 24(10)-** and 18(10)-** cm® sec.~! 
for argon and nitrogen, respectively. 





INTRODUCTION 


RESSURE effects of nitrogen on the po- 

tassium resonance lines in absorption have 
already been investigated by Watson and Mar- 
genau.' However, they were able to obtain only 
three spectrograms between relative density 12 
and 20. It is in this higher pressure region where 
departure of the linear relation between relative 
density and broadening occurs. In order to es- 
tablish more definitely this point of departure, 
and to observe the broadening at higher relative 
densities, the present work was undertaken. Since 
the publication of the paper by Watson and 
Margenau, the latter? has extended further the 
theory of pressure broadening and given a theo- 
retical contour for broadened spectral lines. Con- 
sequently in the present work it was possible to 
compare theoretical with experimental contours. 


EXPERIMENTAL PROCEDURE 


A detailed description of apparatus is given 
elsewhere.* It consisted of a steel absorption tube 
87 cm long, fitted with 1 cm thick Pyrex glass 
windows. A three-section electric furnace served 
to vaporize the potassium in the absorption tube. 
Argon and nitrogen were used as perturbing gases 
up to relative density 28 (40 atmospheres). Above 
this the potassium doublet coalesces making 
measurements of half-width inaccurate. Pressures 
were reduced to relative densities by the law for 

* Part of a dissertation to be presented to the Faculty 
of the Graduate School of Yale University in candidacy 
for the degree of Doctor of Philosophy. 

1'W. W. Watson and H. Margenau, Phys. Rev. 44, 748 
(1933). 

*H. Margenau, Phys. Rev. 48, 755 (1935). 

3H. Margenau and W. W. Watson, Phys. Rev. 44, 92 
(1933). 


perfect gases. The more exact van der Waals law 
was not used, since for even the highest pressures 
the correction term was found to be negligible. 
Correction was not made for thermal transpira- 
tion, as was done by Watson and Margenau.' 
Both experiment and theory have shown‘ that no 
correction for this phenomenon is necessary if the 
diameter of the outlet to the pressure gauge is 
many times greater than the molecular mean free 
path. In the apparatus used, the diameter of the 
outlet was several mm and the mean free paths of 
the order of (10)-* mm or less. 

Spectrograms of the resonance lines (7665A 
and 7699A) were taken on Eastman type 144-N 
plates in the first order of a 21-foot concave 
grating in a stigmatic mounting. A Pointolite 
lamp was used as continuous light source, and 
second order Fe lines served for the comparison 
spectrum. The characteristic log J vs. D curve of 
the spectrograms was obtained from the blacken- 
ing marks produced by a step filter. This con- 
sisted of a quartz disk on which were deposited 
five strips of platinum of different densities and 
had been previously calibrated for 47650. The 
filter was placed directly in front of the plate and 
just to one side of the resonance lines. Both 
blackening marks and potassium lines were re- 
corded simultaneously on the photographic plate, 
thus insuring the same continuous spectrum as 
well as the same exposure time for each. This 
method of obtaining the characteristic curve has 
the advantage of not relying on the reciprocity 
law. Microphotometer traces of the spectrograms 
were made on a Koch-Goos self-recording instru- 
ment, and those line contours which did not lie on 


4L. B. Loeb, Kinetic Theory of Gases (1927), p. 283. 
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Fic. 1. Wave-lengths of the maxima of the potassium reso- 
nance lines vs. relative density of argon and nitrogen. 


the straight line part of the characteristic curve 
were corrected to give true absorption contours. 
The wave-length of the maximum of each absorp- 
tion curve as well as its half-width was calcu- 
lated from the dispersion on the microphotometer 
trace as determined by interpolation between 
Fe lines. 
RESULTs 

The results are shown graphically in Figs. 1 
and 2. In Fig. 1 the wave-lengths of the maxima 
of the potassium resonance lines are plotted 
against relative densities of argon and nitrogen, 
and in Fig. 2 the averaged half-widths of the 
doublet are plotted. (No significant difference 
between the half-widths of the doublet com- 
ponents was detected.) Curves have been drawn 
through the ‘“‘mean”’ of the plotted points. The 
points in Fig. 1 lie more closely on the curves 
than do those in Fig. 2, which is to be expected, 
since wave-lengths of the maxima can be de- 
termined with greater accuracy than _half- 
widths. The latter depend upon corrections for 
nonlinear plate characteristics, which do not 
affect the line maximum. Extrapolating the 
curves back to the axis of ordinates, the values 
7664.84A and 7698.96A are obtained for the 


maxima of the potassium doublet at zero relative 
density, and 0.25 cm for the half-widths. It is 
seen that argon produces a greater half-width and 
shift of the maximum than nitrogen, and the 
shift of 47699 is slightly greater than \7665. The 
curves up to relative density 14 are straight lines, 
indicating that in this pressure region ‘‘velocity"’ 
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Fic. 2. Half-widths of the potassium resonance lines vs. 


relative density of argon and nitrogen. Points are averaged 
values for the two components of the doublet. 


broadening® predominates. At about relative 
density 14, the curves depart from linearity, 
showing that ‘“‘statistical’’ broadening® is be- 
coming effective. This is in accordance with 
Margenau’s theory.? and also with the previous 
work done by Watson and Margenau' using 
nitrogen as the perturbing gas. However, shifts 
and half-widths obtained for nitrogen in the 
present work are nearly twice as great as their 
values. Half this discrepancy can be accounted 
for by the fact that they made correction for 
thermal transpiration in reducing pressures to 
relative densities. In Table I are listed shifts and 
TABLE I. Summary of shifts and half-widths of the potasstum 


resonance doublet. 


Shift Half-Width 


Av Ay 
d d 
7665A 0.360 cm™ ) 
Argon > 1.01 cm™ 
7699 0.420 
7665 0.300 
Nitrogen 0.82 
7699 0.360 


’ H. Margenau and W. W. Watson, Rev. Mod. Phys. 8, 
22 (1936). 
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Fic. 3. Experimental (solid) and corresponding theo- 
retical (dashed) contours of \7699 perturbed by nitrogen 
(curves A), and by argon (curves B), relative densities 3.5 
and 5.9, respectively. 
half-widths per unit relative density as de- 
termined by the slope of the straight line part of 
the curves. According to theory, the half-width 
should be approximately twice as great as the 
shift when velocity broadening predominates. 
From Table I, we obtain Av; Av=2.6 for argon 
and 2.5 for nitrogen. 

To compare experimental with theoretical line 
contours, two microphotometer curves in the 
lower pressure range were chosen which did not 
need correction for the characteristic of the 
photographic plate. These were \7699 contours 
for argon and nitrogen, relative density 5.9 and 
3.5, respectively. In this pressure range, Mar- 
genau gives for an approximate theoretical 
contour 2 


1 Tr? , rh? i 
wr$o 2) of 2))} 
A v’—A/2 v’+A/2 


o(x) =(2 rf e~ “dy. 
0 


where 


Here, v’ is the frequency measured from the line 
origin vo, A the half-width, and \= 327 | b| '» where 
nis the number of perturbers per unit volume and 
ba constant, which is a measure of van der Waals 
forces. Eq. (1) was evaluated by numerical in- 
tegration for various values of A and rd?. The 
best fit with the experimental contour for argon 
was obtained with 7\?=0.5A where A=4.28 cm7 
and for nitrogen 7A? =0.4A, A=2.50 cm~. In Fig. 
3 experimental (solid) and theoretical (dashed) 
contours are shown fitted together; curves A are 


for nitrogen, B for argon. Agreement is good on 
the long wave-length side of the maximum, but 
on the short wave-length side the theoretical 
curves fall off too rapidly. This is to be expected, 


GORDON F 
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TABLE II. Values of the van der Waals constant 6 calculated 
by three different methods. 


van der Waals Constant } 





Calculation by Eq. Argon Nitrogen 
(1) 24(10)-82 18(10)~*? cm® sec.~! 
(2) 24 17 


(3) 23 18 


however, because of the approximations made in 
deriving Eq. (1). 

The relation between rd? and A for the theo- 
retical contours which best fit those of experi- 
ment, affords a means of calculating the van der 
Waals constant b. The values of 6 obtained in this 
manner are given in Table II. As a check on these 
values, two other methods of calculating 6 are 
available. The first is to apply Lorentz’ formula 
as modified by Weisskopf :* 


Av; =2.2b?/5( V)3/5n, (2) 


where V is the molecular root mean 
velocity, and where we shall take for Av; the 
half-widths of the theoretical curves. Secondly 
Margenau’s theory states? that the curves of half- 
width and shift versus relative density depart 
from linearity when the velocity broadening half- 
width becomes approximately equal to the sta- 
tistical half-width; i.e., when 


square 


Av; =r’. (3) 


Experimentally this point of departure is at 
about relative density 14, and the half-widths of 
the potassium doublet perturbed by argon and 
nitrogen at this relative density are 12.34 cm“ 
and 11.29 cm“, respectively. Values of 6 calcu- 
lated from (2) and (3) are also given in Table IT. 
The results in Table IT are in excellent agreement 
among themselves, but } for nitrogen does not 
agree with the value [7(10)-* cm® sec.~'] calcu- 
lated by Margenau? from the previous work on 
the broadening of the potassium resonance lines. ' 
The reason for this discrepancy has already been 
noted; the half-widths obtained by Watson and 
Margenau are about half as large as obtained in 
the present work. 

The author wishes to thank Professor W. W. 
Watson for suggesting this problem and for his 
helpful criticism during the course of the work, 
and also Professor H. Margenau for his advice in 
the theoretical interpretation of the results. 


®° \V. Weisskopf, Zeits. f. Physik 75, 287 (1932). 
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The Infrared Absorption of Alcohol-Acetone Mixtures 


WaLTER Gorpy,* Department of Physics, University of North Carolina, Chapel Hill, North Caroltna 


Received October 12, 1936) 


Several mixtures of methyl alcohol and acetone have been studied in the region 2.5 to 8.5. 
The acetone was found to shift the OH vibrational band of alcohol to the shorter wave-lengths 
and to increase appreciably its intensity; the effects of the alcohol on acetone were to increase the 
intensity of the vibrational band of the CO group and to shift it to the longer wave-lengths. 
No other bands were measurably affected. The results indicate association of the alcohol and 
acetone through interaction of the OH and CO groups. The possibility of the formation of 


proton bonds is suggested. 





NFRARED absorption studies' of aqueous 

and alcoholic solutions of inorganic acids, 
hydroxides and salts have revealed bands which 
indicate molecular association between the sol- 
vent and the solute molecules. More recently 
infrared studies? have been made of mixtures of 
water and some organic liquids. The present 
work on alcohol-acetone mixtures was under- 
taken to ascertain the influence of each liquid on 
the spectrum of the other in various mixtures of 
alcohol and acetone. 

A Hilger infrared spectrometer with a fluorite 
prism was used as a resolving instrument 
throughout the investigation. Cell windows were 
of fluorite, and the absorbing layer was main- 
tained at a constant thickness of 0.002 cm by 
mica washers. The chemicals Baker's 
analyzed. The methyl alcohol was dried for 
several hours over CaO, and the acetone over 


were 


CaCl, as an additional precaution against the 
presence of water in the samples. 

A study was made of the 2.94 alcohol band and 
of the acetone bands appearing at 5.73u and at 
8.174 for several methyl alcohol-acetone mix- 
tures. Fig. 1 shows the transmission curves for 
methyl alcohol and acetone and for four methyl 
alcohol-acetone mixtures in the region 2.54 to 
3.154. The broad band appearing in the region of 
2.91 in the methyl alcohol curve is charac- 
teristic of all alcohols, and arises from changes in 


* Now at Mary Hardin-Baylor College. 

1E. K. Plyler and W. Gordy, J. Chem. Phys. 2, 470 
(1934); E. K. Plyler and E. S. Barr, J. Chem. Phys. 2, 
306 (1934); E. K. Plyler and F. D. Williams, J. Chem. 
Phys. 2, 565 (1934). 

* D. Williams and E. K. Plyler, J. Chem. Phys. 4, 154 
(1936); D. Williams, R. D. Weatherford and E. K. 
Plyler, J. Opt. Soc. Am. 26, 149 (1936). 


the vibrational energy of the OH group. For 
alcohol concentration of 12 percent and below, 
this band appears at 2.75u, having been shifted 
to the shorter wave-lengths by 0.15. In addition 
to the shift to higher frequencies, the intensity of 
the OH band is increased appreciably. That is, 
the intensity of the band does not decrease in 
proportion to the decrease in the number of 
absorbers. As the alcohol concentration is in- 
creased to 25 percent and to 50 percent, the band 
approaches that for pure alcohol. It may be that 
the bands for these concentrations are made up 
of two components with maximum of absorption 
at 2.754 and 2.94. However, the bands are very 
broad, and it was not found possible to resolve 
them into two components. Fig. 2 shows the 
effects of alcohol on the vibrational band of the 
CO group of acetone. In contrast to the shift to 
the shorter wave-lengths of the OH alcohol 
band, the effect of the alcohol on the acetone is 
to shift the CO band to the longer wave-lengths. 
This shift is accompanied by an increase in 
intensity, as was the case of the OH alcohol band. 

To show that the above mentioned variations 
of the OH and CO bands are not merely apparent 
effects due to the superposition of alcohol and 
acetone absorption, theoretical curves have been 
computed for mechanical mixtures in which no 
interaction occurs. In Figs. 3 and 4 these curves 
are shown in comparison with the experimental 
curves for the same mixtures. We have assumed 
that the transmission of the mixture for a given 
frequency is 


oT” sam Simi—Som 
T=e 2, 


where f, and fz are the fractional parts of alcohol 
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Fic. 1. The percent transmission of alcohol and acetone 
and alcohol-acetone mixtures in the region 2.54 to 
Su. 


and acetone, respectively, and where m, and mz 
are defined by 


Ti=e™, T2=e-™, 


T, and T, being the transmission of pure alcohol 
and pure acetone, respectively, for the given 
wave-lengths. It should be remembered that the 
cell thickness was the same in all cases. 

From Fig. 5, which gives the absorption of 
several alcohol-acetone mixtures in the region 
7.85u to 8.45y, it will be seen that the strong 
acetone band appearing at about 8.174, bottom 
curve, apparently is not shifted in position for 
any of the mixtures, nor is its intensity appreci- 
ably increased as was the case of the CO acetone 
band and the OH alcohol band. A curve was 
obtained for 6 percent acetone in alcohol which 
is not included in the figure. The intensity of the 
band was so low that the position of maximum 
absorption could not be ascertained. 


GORDY 


An attempt was made to study other bands in 
the region 2.54 to 9u, but because of the over- 
lapping of alcohol and acetone bands no definite 
results could be obtained. A study was made of 
the effects of benzene on the CO acetone band 
for concentrations of acetone in benzene varying 
from 50 percent to 6 percent. There appeared to 
be no shift in the CO band for any of the concen- 
trations. A similar study was made of the effects 
of benzene on the OH band of alcohol. Likewise, 
no measurable shift was observed in the OH 
band, although there appeared for low alcohol 
concentrations a slight depression on the short 
wave-length side of the band, which may indi- 
cate the presence of vapor-like molecules of 
alcohol. Evidence for vapor-like molecules of 
alcohol in carbon tetrachloride solutions has 
been reported by Kinsey and Ellis.’ 

The pronounced changes in the electric mo- 
ments of the CO and OH groups of acetone and 
alcohol indicate dipole interaction. It may be 
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3 E. L. Kinsey and J. W. Ellis, Phys. Rev. 49, 105 (1936) 
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that we have here the formation of proton bonds, 
the hydrogen of the OH group being shared 
between the oxygen of the two groups. This 
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would account for the extreme changes in the 
polarization of these groups. Similar results‘ 
indicating dipole interaction in dioxane-water, 
methyl cyanide water, acetone-water mixtures 
and in mixtures of aniline with several other 
organic liquids may be explained by the proton 
bond theory. This theory has received other 
emphasis recently.6 Abnormal results obtained 
in dielectric constant measurements of liquid 
mixtures are likely due to similar changes in the 
electric moments of certain groups which may be 
observable by infrared absorption measurements. 
The work is being continued. 

The writer wishes to thank Dr. E. K. Plyler 
for the use of his instrument and Dr. J. A. 
Wheeler for suggesting a convenient method of 
making the theoretical computations. 





* Unpublished data obtained by the author. 
°>L. Pauling, J. Amer. Chem. Soc. 57, 2680 (1935); 58, 
94 (1936); L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 
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The thermionic emission from round filaments and flat strips of platinum was measured 
giving a value for the work function of 5.32 volts and 32 amperes/cm? deg? for the constant A. 
Highest possible vacuum conditions and heat treatment near 2000°K were found necessary 
for obtaining reliable results. The lack of such extreme heat treatment probably accounts for 
the abnormally high values of the constant A formerly found. Special theoretical explanations 


for platinum appear no longer necessary. 





ECAUSE of the wide variation found experi- 
mentally for the constants in the thermionic 
equation, the study of the thermionic emission of 
platinum has been of special interest. Experi- 
mental values for the work function have ranged 
from 2.18 to 6.62 volts, and the values for the 
constant A from 10.7X10- to 10.710’. 
Dushman and DuBridge have listed the results 
of various experimenters.!:? 

Since simple theory predicts A equal to 120 
amperes/cm? deg? for all pure metals, many 
attempts have been made to explain differences 
actually found. Simple theory assumes a clean 
smooth metal surface, a negligible reflection 
coefficient for electrons, and a constant number 
density of electrons within the metal. If any of 
these conditions are not satisfied, one may expect 
variations in the values of A. Of course, most 
experimenters have attempted to have clean 
metal surfaces by using only the purest platinum 
obtainable, and subjecting it to long heat treat- 
ment under the best possible vacuum conditions. 

If the surface of the metal is not perfectly 
smooth, the increased area due to the surface 
roughness would lead to an A larger than 
normal. Fowler* considered such a possibility as 
this in an effort to account for DuBridge’s value 
of A, which was 1.7X10*. This value would 
require a surface roughness factor of about 120, 
and a surface under such conditions would have 
to have a sort of sponge-like structure. Fowler 
pointed out that even if this were so, the acceler- 
ating field could have practically no effect for 
most of the assumed additional surface. At most, 
surface roughness could account for an increase 
in A of only about 50 percent. 


1S, Dushman, Rev. Mod. Phys. 2, 381 (1930). 
?L. A. DuBridge, Phys. Rev. 31, 236 (1928). 
’>R. H. Fowler, Proc. Roy. Soc. 122, 36 (1929). 


A value of A smaller than 120 could be caused 
by a reflection coefficient for electrons greater 
than zero, but there is reason to believe that 
with a clean metal, such a coefficient could have 
but little effect, and probably could not be 
detected experimentally. 

Bridgman‘ developed a theory to explain large 
departures of A from the universal value by 
connecting this departure with a temperature 
derivative of the photoelectric frequency and a 
nonvanishing of (S,—S,,). This theory assumes 
that the surface of the metal has reached an 
equilibrium condition. DuBridge® points out 
that on the basis of this theory the thermionic 
equation can be written 


I=A,T* exp —(b—aT)/T. 


This equation is consistent with the observed 
linearity found for different values of log A and 
the work function as the specimen goes through 
succeeding stages of heat treatment. 

Reimann® has also suggested a method of 
explaining departures of A from the universal 
value by allowing for the change in electron 
number density with temperature. On the basis 
of this theory the form of the thermionic equation 
remains the same, but an additional factor is 
introduced which, multiplied by the theoretical 
A of 120, gives a new effective A. For platinum 
the effective A turns out to be about 27. It 
should be added that this method is open to 
question from the theoretical point of view, 
since the expansion of the metal with tempera- 
ture was allowed for in the final thermionic 
equation rather than earlier in the derivation. 


4P. W. Bridgman, Phys. Rev. 31, 862 (1928). 
°>L. A. DuBridge, Proc. Nat. Acad. Sci. 14, 788 (1928). 
6 A. L. Reimann, Nature 133, 833 (1934). 
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DuBRIDGE’s AND VAN VELZER'S RESULTS 


Most of the older thermionic work can be dis- 
missed entirely because of poor vacuum tech- 
nique. However, the results obtained recently by 
DuBridge and Van Velzer are of interest. In 
DuBridge’s’ latest work he reports a value of 
6.27+0.07 volts for ¢ and 1.7X10* for A. His 
tube was first baked for a short period, then the 
nickle cylinder heated with an induction furnace 
until no more gas was given off, then the tube 
rebaked for 24 hours, and finally the cylinder 
reheated with the induction furnace. During this 
time his platinum strip was at about 1573°K, 
occasionally being raised to 1773°K. His readings 
covered from about 1400°K to 1750°K. 

Van Velzer® heat treated his specimens for 
longer periods and went to higher temperatures. 
He found abnormally high values of A and ¢ 
which persisted until he raised his filament to 
1990°K. After a short heat treatment at this 
temperature, he obtained ¢ equal to 5.4 volts 
and A equal to 170+20. 


APPARATUS FOR PRESENT WORK 


For measuring thermionic emission, single 
platinum filaments were mounted in the center 
of a three-section guard ring type molybdenum 
cylinder. The diameter was 25 mm and the 
central section 25 mm long. Thermionic currents 
were measured with a shunted galvanometer 
having a sensitivity of 4.4 10~-'° amperes/mm. 
Storage batteries were used for the plate poten- 
tial and for the filament current. 

Temperatures were measured with an optical 
pyrometer and readings were taken through a 
flat Pyrex window sealed to the experimental 
tube. To obtain a temperature scale, the emis- 
sivity of platinum was measured using rolled 
cylinders as described by the writer. The tem- 
perature scale found agreed with that given by 
Worthing.'° 


RESULTS 


The thermionic emission from five specimens 
was measured. The first specimen was a round 
12 mil wire, the purest platinum obtainable from 

7L. A. DuBridge, Phys. Rev. 32, 961 (1928). 

8H. L. Van Velzer, Phys. Rev. 44, 831 (1933). 


*L. V. Whitney, Phys. Rev. 48, 458 (1935). 
10 A. G. Worthing, Phys. Rev. 28, 174 (1926). 
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PLATINUM 


the American Platinum Company. The experi- 
mental tube was baked for over 1000 hours at 
450°C, and pressures while baking were finally 
better than 510-7? mm. Heat treatment of the 
filament 1600 hours. 
After the first 500 hours steady heat treatment 
was maintained at 1640°K; after 1000 
hours the temperature of heat treatment was 


was continued for over 
about 


gradually increased to about 1730°K, where it 
was maintained for the last 400 hours. No 
readings from this specimen were included until 
after more than 1000 hours of heat treatment, 
and when readings were taken pressures as read 
by an ion gauge were better than 310-5 mm. 
The values of the thermionic constants showed 
no change during the last 700 hours of heat 
treatment. With readings taken up to 2000°K, 
A was found to be 35 and ¢ equal to 5.32 volts. 

Air was then let into the system, and a side 
tube with an aluminum getter sealed on. Both 
DuBridge and Van Velzer had sealed off their 
experimental tubes with getters, and this pro- 
cedure was used in the present work as a further 
check. Within experimental error, the results 
agreed with those found previously, ¢ and A 
being 5.31 volts and 32, respectively. 

The filament was then removed, rolled flat to 
a thickness of about 0.05 mm, and remounted. 
Baking was continued for 200 hours and heat 
treatment for 350 hours. Readings were taken, 
but were given little weight since a temperature 
gradient, due to unevenness of the filament after 
rolling, was in evidence along the length of the 
filament. A simple theory was developed to take 
into account this gradient. Readings were suf- 
ficiently good to show that the thermionic 
constants were similar in magnitude to those 
found before. A measured 70 and ¢ 5.4 volts. 

The second specimen was also flat, but it was 
rolled from an non-heat-treated round filament. 
This filament appeared perfectly uniform and 
there was no temperature gradient evident. Heat 
treatment was continued for only 270 hours, but 
the filament was flashed to high temperatures for 
short periods for forced aging. A was 33 and ¢ 
5.28 volts. 

The unrolled specimens, and specimens rolled 
before and after heat treatment, were measured 
to discover any differences that might result from 
preferred” crystal orientation. The thermionic 
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Fic. 1. Plot of log. i—2 log. T against 1/T. 


constants from one face of a crystal might differ 
from those of another face. Rolling might tend to 
orient the crystals and result in certain preferred 
crystal faces being exposed to the outside. 
However, no differences in thermionic constants 
could be detected. 

The third specimen was a strip from the supply 
of platinum originally used by DuBridge at 
Wisconsin.» ® Spectroscopic tests of this plati- 
num showed traces of palladium which were not 
present in the other specimens used in this work. 
Due to the thinness of this third specimen, it did 
not stand the rigorous heat treatment given the 
others. The tube was baked for 120 hours and 
heat treatment of the filament continued for 
nearly 300 hours. Readings were made with 
pressures better than 3X10-* mm. High values 
of A were found which dropped with heat treat- 
ment from more than 10° to 310%. The final 
value of ¢ found for this specimen was 5.43 
volts. These values are probably not character- 
istic of the clean metal as the specimen burned 
out before sufficiently high temperatures could 
be reached for steady and repeatable readings. 

The fourth sample was a round filament from 


tL. A. DuBridge, Proc. Nat. Acad. Sci. 12,162 (1926). 
2L. A. DuBridge, Phys. Rev. 29, 451 (1927). 


the same supply as the first and second, and was 
measured to find the effect of a shorter period of 
heat treatment but at higher temperatures. The 
temperature was raised to 1890°K and readings 
taken after the first 100 hours of heat treatment. 
Readings at this time were comparatively stable, 
¢g measuring 5.8 volts and A 1.510%. These 
results, which are consistent with those obtained 
by Van Velzer for similar heat treatment, are 
shown on curve A of the Richardson plot, Fig. 1. 
After this heat treatment, thermionic currents 
showed an almost imperceptible but definite 
decrease with time. However, after final heat 
treatment at 2000°K, thermionic currents ap- 
peared stable and values of A and ¢ dropped to 
31 and 5.31 volts, respectively. Heat treatment 
was continued for over 100 hours after these 
lower values were found, but there was no notice- 
able further change. 

The fifth specimen was a flat strip specially 
purified and kindly furnished by Dr. J. A. Becker 
of the Bell Laboratories. This burned out twice 
before readings above 1900°K could be obtained, 
but fortunately the failures were near the ends 
both times and the specimen cou!d be remounted. 
After heat treatment at 1900°K, A measured 10* 
and ¢ 5.75 volts; heat treatment at 1930°K, 
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A became 300 and ¢ 5.6 volts. Still lower values 
probably would have been found, but the filament 
failed when raised to higher temperatures. 

Curve B, Fig. 1, of the Richardson plot is 
based on the data from the three specimens, 
No. 1, No. 2, and No. 4, which withstood the 
severe heat treatment. The average from these 
specimens give a value of 5.32 volts for g and 32 
for A uncorrected for any surface roughness. 
Results from the other specimens indicate that 
similar values were being approached. The linear 
relationship between log A and ¢, which has been 
discussed at length by DuBridge,® is seen to 
hold as shown by Fig. 2. 


The value of A found is seen to be well within 
the region predicted by theory. Special explana- 
tions to justify the exceptional behavior of 
platinum now seem unnecessary. The present 
work indicates the necessity of the best possible 
vacuum conditions, and verifies Van Velzer’s 
conclusions as to the importance of high tem- 
peratures in the heat treating process. Appar- 
ently, heat treatment at nearly 2000°K for short 
periods is necessary for dependable results. The 
writer wishes to thank Professor H. B. Wahlin, 
under whose direction this work was done, for 


his suggestions and assistance. 
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The Variation of the Internal Friction and Elastic Constants with Magnetization 
in Iron. Part I 


W. T. Cooke, Columbia University, New York, New York 
(Received October 5, 1936) 


A new and precise method is described for measuring 
the elastic constants and coefficients of internal friction of 
solid substances. The method is applied to investigate the 
variation of Young’s modulus and the longitudinal coeffi- 
cient of internal friction of annealed and unannealed 
Armco iron with magnetization at room temperature. Data 
are given which describe the behavior of the iron in three 
different magnetic states, vz., at different points on the 


normal induction curve and on a hysteresis cycle, and in 
the anhysteretic state obtained when an alternating field, 
the amplitude of which is gradually reduced to zero, is 
superposed on a steady magnetizing field. The observed 
coefficients of internal friction of the unannealed specimen 
require correction, and this, together with the torsional 
data and theoretical interpretation of the results, are given 
in Part II of this paper. 





HE amplitude of vibration at all points of 

a free body which is vibrating in one of its 
normal modes diminishes with the time in a 
manner described by the formula 


u=tUy exp (—Aft), (1) 


in which u denotes the amplitude, f the fre- 
quency of vibration, ¢ the time, and A a quantity 
whose value depends, in the first instance, upon 
the substance and upon the nature of the stress 
associated with the particular mode of vibration. 
The name “internal friction’’ has been given to 
the property of solid substances of which this 
attenuation in amplitude is a measure. An alter- 
native view of the phenomenon follows from a 
consideration of the rate of transfer of organized 
vibrational energy into heat.' Thus, if W denotes 
the energy dissipated per cm* per cycle of vibra- 
tion, then it can be shown that, in accordance 
with Eq. (1), 

W=éS,?, (2) 
where So is the stress amplitude, £=A/P, and P 
is the appropriate elastic modulus. 

The coefficients of internal friction, the é’s, 
and the logarithmic decrements, the A’s, of a 
substance resemble the elastic moduli in that 
their definition requires the specification of a 
stress. The types of stress here dealt with are 
the stretch and twist associated, respectively, 
with the longitudinal and torsional vibration of 
a slender bar. The corresponding decrements are 
hereafter denoted by the symbols A” and A®%, 
and, conformably, the Young’s and rigidity 
moduli by the symbols E and G. 

The energy dissipated in a vibrating ferro- 
~ 1 Kimball and Lovell, Phys. Rev. 30, 948 (1927). 


magnetic body receives contributions from three 
distinct sources: first, internal friction of the 
sort present in nonferromagnetic materials; sec- 
ond, eddy currents induced in the medium by 
the change in magnetization produced by the 
varying stress; and third, internal friction of 
purely ferromagnetic origin. No theory of inter- 
nal friction in nonferromagnetic solids has yet 
been written, nor has any attempt been made to 
lay an experimental basis for one. While it is 
hoped that the experimental method here de- 
scribed may encourage such an effort, the present 
research is concerned primarily with the two 
contributions to the internal friction which are 
associated with the magnetization. For, as will 
be shown in Part II of this paper, recent develop- 
ments in ferromagnetic theory permit these to 
be correlated with other magneto-elastic proper- 
ties of the material, and the success of such 
correlation lends added support to the theory. 


EXPERIMENTAL METHOD 


The experimental method utilizes the proper- 
ties of a separately excited composite piezoelec- 
tric oscillator constructed by cementing a suit- 
ably cut right circular cylinder of crystalline 
quartz to one end of a cylinder of specimen 
material of identical cross section. The cement 
is hard shellac, softened by heating, and the join 
is made under pressure. The system is freely 
suspended, and a sinusoidally varying potential 
difference is applied between electrodes of gold 


leaf affixed in proper position to the quartz.’ 


? Details of the method for preparing the quartz crystal 
are given in papers by Balamuth, Phys. Rev. 45, 71 
(1934), and Rose, Phys. Rev. 49, 50 (1936). 
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ELASTICITY AND 
In consequence of the piezoelectric stress which 
accompanies the electric field in the quartz, a 
stationary state of forced longitudinal or tor- 
sional vibration is established in the composite 
oscillator. A system of this sort offers two meth- 
ods for measuring the elastic moduli and co- 
efficients of internal friction of the specimen 
material.* 

The first is that adopted by Zacharias,‘ who 
measured the variation of E and A” with tem- 
perature in unmagnetized single and polycrystal- 
line nickel over the temperature range 30°C to 
400°C, and Siegel,> who measured the variation 
of E and A®¥ with magnetization and temperature 
in polycrystalline nickel over the same tempera- 
ture range. In this method the amplitude of the 
applied potential difference is held constant, and 
means are provided for varying its frequency 
and for measuring the amplitude of the current 
which flows to the oscillator. The current ampli- 
tude varies critically with frequency in the neigh- 
borhood of certain ‘‘resonance frequencies’’ at 
which the amplitude of the forced vibration 
passes through a maximum. Values of the moduli 
and decrements are deduced from the observed 
nature of this variation. The method is generally 
applicable for the measurement of the moduli, 
but the precision of the decrement measurement 
is bad when the value of this quantity is less 
than about 10~-*. Accordingly it is not suitable 
for studying the internal friction of iron and 
most nonferromagnetic solids. 

In the present method the composite oscillator 
forms one arm of an a.c. bridge which is excited 
by an alternating voltage of constant amplitude 
and variable frequency. Values of the moduli 
and decrements are deduced from the observed 
variation of the electrical impedance of the oscil- 
lator with frequency. The precision of the decre- 
ment measurement fails when this quantity is 
large of the order 10-*, so that the two methods 
perfectly supplement each other. 

The electro-mechanical theory of the single 

3 As a matter of fact, the composite piezoelectric oscil- 
lator has been used to measure the internal friction of 
solids by four different methods, of which the two here 
described are regarded as superior. For descriptions of the 
—— see Quimby, Phys. Rev. 39, 345 (1932), and refer- 
ence /. 


4 Zacharias, Phys. Rev. 44, 116 (1933). 
> Siegel and Quimby, Phys. Rev. 49, 663 (1936). 
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FRICTION IN SOLIDS 
crystal oscillator,* and the mechanical theory of 
the two-part composite oscillator’ have been set 
forth in detail by previous writers. The theoreti- 
cal treatment of the present system follows so 
closely that of these others that it need here only 
briefly be outlined. The charge which must be 
placed on the gold leaf electrodes in order to 
establish a potential difference & between them 
can be separated into two parts: first, a part of 
the ordinary capacitative sort, equal to C’&, 
where C’ is the interelectrode capacity ; and sec- 
ond, a part required to neutralize the piezoelec- 
tric charge produced by the vibrational strain 
in the quartz, and proportional to the space 
average value of this strain, taken along the 
cylinder. Accordingly, the current J, which flows 
to the oscillator, may be written in the form 


Il=C'd&/di—K'edS,, /dt, (3) 


where K’ is a geometrical constant, ¢ the appro- 
priate piezoelectric coefficient, and S,, the aver- 
age strain. Now Sy, is proportional to &, so that 
the ratio Z=&/I is independent of & and is the 
electrical impedance of the oscillator. The fur- 
ther development of the theory is concerned with 
the evaluation of Syy. 

In order properly to estimate the effect of the 
layer of cement on the decrement measurements, 
the composite oscillator is initially regarded as 
composed of three coaxial cylinders of identical 
cross section and different lengths and materials. 
The actual piezoelectric stress in the quartz 
cylinder is represented by harmonically varying 
surface tractions over its end faces. The equation 
of motion of each medium is assumed to be of 
the form 


pu dP =P(A+T9d dt)?u dx’, (4) 


where u is the particle displacement, p the den- 
sity, P the appropriate elastic modulus, and T 
a dissipative coefficient which, in accordance 
with the theory of free vibrations in bars, is 
related to A by the formula 


A=2nfT. (5) 


The particle displacement in each medium is 
* Dye, Proc. Lond. Phys. Soc. 38, 399 and 457 (1926); 
Van Dyke, Proc. Inst. Radio Eng. 16, 742 (1928). 

7Quimby, Phys. Rev. 25, 558 (1925); Balamuth, 
reference 2. 
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given by an expression of the form 


u=|A exp (a+i8)x+B exp 
—(a+i8)x} expint, (6) 


where 


n=2rf, B=n(p/P)', a=3npT, 


and it is assumed that 7 is sufficiently small that 
n?T?<1. The six amplitude coefficients (A and B 
for each medium) are evaluated by means of the 
six simultaneous equations which express, re- 
spectively, the continuity of stress and displace- 
ment at the two interfaces and of stress at the 
end faces of the oscillator. S,y is calculated from 
the resulting expression for u in the quartz, and 
the electrical impedance of the oscillator follows 
at once from Eq. (3). 

The expression for Z can be arranged in the 
form 1/Z=inC’+1/Z,. Furthermore, in the 
neighborhood of a resonance frequency the form 
of the expression for Z,, is precisely the same 
as that for the electrical impedance of an in- 
ductance, capacity and resistance in series. It 
follows that, near resonance, a composite oscil- 
lator is electrically equivalent to a series resonant 
electrical circuit of impedance Z,, shunted by 
a capacity C’. If Z, be written in the form 
Zm=R+1X, then, for a three-part oscillator, 


R=(Kn/D){(M21/y1)(sec.* yi — p1) 
[1 —(Moys?/ Ms) pops |+ (Moz2/y2) 
»4 (sec.? V2 —p2)[1 ~“— (Myy;?, ‘Ms) pips | 
+ (Ms23/‘ys)(sec.? ys — ps) 
<[1—(M, Meys?/ M3?) pip) |}, (7) 
X =(Kn/D){ Mipit+ Mepe+ Mops; 
—(M,Moys?/Ms3)pipeps}, (8) 
where 
D=2(1—sec. ve)[1— (Miys?/ M3) pips ] 
—_ (pove" “‘M2)(Mipi a M3ps ), 
K =a constant,® 
M,;=mass of a cylinder, (¢=1, 2, 3) 
L;=length of a cylinder, 


si:=aiLi, 

yi=h/fi, 

fi=(1/2L,)(Pi/p;)', 

pi=(tan y,)/¥i, (9) 


and the subscripts, 1, 2, 3, refer, respectively, to 
the specimen material, the quartz, and the 


* For longitudinal vibrations, K =1/e, where 6 is the 
width of a gold leaf electrode. 


COOKE 


cement. It will be noted that the quantities f; 
are the fundamental frequencies of free vibration 
of the cylinders separately. The resonance fre- 
quencies of the oscillator are those at which 
X=0. 

In practice the lengths of the quartz and speci- 
men cylinders are so adjusted that f> is equal to 
an integral multiple of f; within +50 cycles at 
50 ke. Furthermore, the thickness, L;, of the 
layer of cement is small of the order 10-* cm. 
A scrutiny of Eqs. (7) and (8) reveals that, under 
these circumstances and near a resonance fre- 
quency of the system, the presence of the cement 
is without appreciable effect upon the values of 
Rand X. Accordingly, the resonance frequencies 
are the solutions for f of the equation 


Mipit+ Meope=0. (10) 


Again, the variation of X with frequency near 
resonance is given by the formula 


X =39K(M,+ M2)éf, (11) 


where 6f is the departure from the resonance fre- 
quency, and the resistance of the oscillator at 
resonance by the formula 


R={K(M,Ai+ MeA2) fo, (12) 


where fo is the resonance frequency. 

The equivalent electrical inductance L, and 
capacity C, of the two-part oscillator are related 
to X by the formulae L=1/n,°;C=X/4rdf, from 
which and Eq. (11), 


L=}K(M,+M,), (13) 
and 
C=8/Kn,?(M,+™.). (14) 


The experimental procedure for obtaining the 
elastic modulus and decrement of the specimen 
material is as follows: It will be noted that the 
behavior of the quartz crystal alone is described 
by Eqs. (10) to (14) with M, set equal to zero. 
Hence f2 is the observed resonance frequency of 
the crystal vibrating alone, and A; is related to 
the observed resistance of the crystal by the 
equation R=}KM2Acfe. K is evaluated with 
Eq. (11) from the observed variation of the react- 
ance of the composite oscillator with frequency. 
Finally, P; and A; are calculated with Eqs. 
(9), (10) and (12), from the observed resonance 
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CYCLES DEPARTURE From RESONANCE 


Fic. 1. Comparison of the observed and calculated values 
of the resistance and capacity required to balance an a.c. 
bridge in one arm of which is a composite piezoelectric 
oscillator. 


frequency and resistance of the oscillator at 
resonance. 

The a.c. bridge employed is similar to that 
described by Stratton.® Essential features are 
adequate shielding and the use of a Wagner 
ground system. The ratio arms are those of a 
General Radio Co. type 210 unit, with an auxili- 
ary variable condenser connected in parallel with 
the coil of lower resistance to compensate the 
distributed capacity of the other and the leakage 
resistance between the electrodes on the quartz. 
A condenser of about 30 cm capacity is con- 
nected in parallel with the oscillator in order that 
the reactance of this arm should be capacitative 
over the entire resonance range. The fourth arm 
of the bridge is a variable condenser connected 
in parallel with a Leeds and Northrup No. 4746 
decade resistance box. Reliable measurements 
can be made of resistances as high as six meg- 
ohms at 50 kc. 

Power is supplied the bridge by a very stable 
vacuum tube oscillator’ and amplifier. The fre- 
quency of the oscillator is controlled, in part, by 
a vernier condenser whose full scale represents a 
frequency change of 30 cycles at, roughly, 60 kc. 
Frequency differences can be measured with an 
accuracy of one-twentieth of a cycle. The vernier 
is frequency calibrated at intervals of 3 to 5 
cycles by beating high harmonics of the oscillator 
against those of a piezoelectric clock." 


® Stratton, J. Opt. Soc. Am. 13, 471 (1926). 
10 King, Beil Sys. Tech. J. 2, No. 4, 31 (1923), Fig. 66. 
1 Dye, reference 6; Quimby, reference 3. 


FRICTION IN SOLIDS 1161 


f 
v O} 
Fao ‘ 2445 
: rie 
xe 4 
be | e xt o 
Q x Ay c 
“4 “ 164 
ul . 8 
5 Fel 
S| 4 
=| F 
oO As 
> 2 a 
e . 
| - J 
| 10 20 30 40 





MASS w GRAMS 


Fic. 2. The resistance of a composite oscillator varies 
linearly with the mass of the specimen cylinder, and its 
inductance with the sum of the masses of the quartz and 
specimen cylinders. 


The apparatus functions remarkably well. The 
observations shown in Figs. 1 and 2 were made 
on a brass oscillator of square cross section 5 mm 
on edge. The brass cylinder was, initially, 18.64 
cm long, and the resonance frequency 56,251 
cycles. The points of Fig. 1 indicate the observed 
values of the resistance and capacity required to 
balance the bridge in the neighborhood of reso- 
nance. The curves, as drawn, represent the for- 
mulae of the theory. The agreement is note- 
worthy. The graphs of Fig. 2 show, in accordance 
with Eqs. (12) and (13), a linear variation of R 
with M, and L with (14,+ M2). These data were 
secured with a succession of six oscillators, each 
obtained from its predecessor by sawing a half- 
wave of vibration from the brass cylinder, the 
resonance frequency thus remaining constant. 

Fig. 3 is a cross section of the assembly for 
measuring the variation of the elastic moduli and 
internal friction of iron with magnetization. The 
oscillator is mounted in an evacuated glass tube 
on silk threads located at displacement nodes of 
vibration. The tube is surrounded, in turn, by a 
water jacket and by the magnetizing coils. The 
latter are designed to secure, as nearly as pos- 
sible, a uniform magnetization over the length 
of the specimen. The end coils, B of Fig. 3, are 
each 190 turns of 5 mil copper strip 5 cm wide. 
The middle coil is six layers of No. 16 D.C.C. 
wire 21 cm long. Typical magnetic flux distribu- 
tion curves are shown in Fig. 4 

Numerous experiments on iron oscillators of 
different lengths, oscillating at different har- 
monics of the fundamental mode, indicate that 
the experimental error in A; due to nonuniform 
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Fic. 3. Cross section of the assembly. 


magnetization near the ends of the cylinder is 
never greater than 5 percent.” An error of about 
the same magnitude may occur in the measure- 
ment of the percent change of the elastic modulus 
with magnetization. 


RESULTS 

Specimen materials 

The specimen materials upon which the obser- 
vations here reported were made are cold-rolled 
Armco iron (hereafter designated “‘unannealed’’), 
and the same substance after annealing at 930°C 
for 2 hours in hydrogen at atmospheric pressure 
and cooling to room temperature in 8 hours 
(hereafter designated ‘‘annealed’’). The modulus, 
decrement and Wiedemann effect measurements 
were made on cylinders of circular cross section, 
0.3 cm to 0.5 cm in diameter and about 30 cm 
long. The resonance frequencies at which the 
observations were taken lay in the neighborhood 
of 56 ke for the longitudinal vibrations and 39 kc 
for the torsional, so that there were seven half- 
waves of vibration in the specimen. The mag- 
netic and magnetostriction measurements were 
made on an ellipsoid of revolution whose semi- 
axes are 0.197 cm and 15.24 cm. 


Magnetic measurements 


Modulus and decrement measurements were 
made subsequent to three different magnetic pro- 
cedures: first, the iron was carried through many 
hysteresis cycles of magnetization at a chosen 








2 The same experiments reveal that, over the frequency 
range 36 kc to 72 kc, the coefficient of internal friction of 
unmagnetized, unannealed Armco iron is proportional to 
fi. This is consistent with the data of Wegel and Walther 
(Physics 6, 141 (1935)) who measured the frequency 


variation of the internal friction of many solids, and found 
values of the exponent of f ranging from +0.54 to —0.31. 





maximum value of the applied field, and left in 
the state corresponding to that maximum value 
(the data are hereafter designated ‘normal in- 
duction’’); second, the iron was carried through 
a hysteresis cycle of magnetization and left at 
a chosen point on the cycle (the data are here- 
after designated ‘“‘hysteresis’’); third, the iron 
was placed in a steady field of chosen value, 
and a 60-cycle alternating field was superposed 
thereon, the amplitude of which was gradually 
reduced to zero (the data are hereafter desig- 
nated ‘‘anhysteresis’’). 

The normal induction relations between mag- 
netization J, and magnetizing field //7, in the 
annealed and unannealed specimens are given in 
Table I. The hysteresis and anhysteresis data 
insofar as they differ from the normal induction, 
are given in Tables II and III. Values of J corre- 
sponding to values of H7 higher than those which 
appear in these two tables are in each case to be 
obtained from Table I. 
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Fic. 4. Distribution of magnetic flux in the specimen. 
The curve with the crosses shows the distribution when the 
specimen is placed in a long solenoid. 
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Fic, 5. The variation of Young’s modulus with magnetiza- 
tion in unannealed and annealed Armco iron. 


The coercive force, remanence, and initial sus- 
ceptibility of the unannealed material are, re- 
spectively, 2.3, 570 and 7.4; and of the annealed 
material 0.46, 1130 and 28.7. 


Young’s modulus and longitudinal decrement 


In order to test the validity of the theoretical 
expression for the contribution of the eddy cur- 
rents to the internal friction, observations were 
made on two unannealed specimen bars of the 
same length and different diameters, viz., 0.538 
cm and 0.397 cm. Of these only the larger was 
subsequently annealed, since it was found that 
the internal friction in the annealed material is 
so large that the eddy current contribution is 
almost negligible. 

The curves of Fig. 5 show the normal induction 
variation of Young’s modulus £, with magneti- 
zation in unannealed and annealed Armco iron. 
The results obtained with the cylinders of differ- 
ent diameter are indistinguishable, as they should 
be, and supply a check on the method. The an- 
hysteresis variation of E with J is the same as the 
normal induction. The value of Young’s modulus 
is 18.610" dynes/cm?® for the unmagnetized 
unannealed material and 19.910" dynes /cm?* 
for the unmagnetized annealed material. 

The points and circles of Fig. 6 show the ob- 
served variation of the longitudinal decrement 
A¥, with normal magnetization in the two un- 
annealed specimens. The circles of which the 
right or left halves are blackened indicate the 


differences between the observed decrements of 
the two specimens and the contributions arising 
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Fic. 6. The variation of the longitudinal internal friction 
with magnetization in unannealed Armco iron. The half- 
blackened circles indicate values of the decrement after the 
contribution arising from eddy currents has been subtracted. 


from eddy currents, when the latter are calcu- 
lated in the manner described in Part II of this 
paper. It is clear that all effect of the form of the 
specimen is removed by this calculation, and that 
the course of the corrected curve in fact repre- 
sents the variation with magnetization of that 
part of the internal friction which is of purely 
ferromagnetic origin. 

The anhysteresis longitudinal decrement in the 
unannealed material differs inappreciably from 
the normal induction. The variation of the hys- 
teresis decrement with magnetization is nearly 
the same at high magnetizations, but at low 
magnetizations the hysteresis values are all larger 
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Fic. 7. The variation of the longitudinal decrement with 
magnetization in annealed Armco iron. 


than the normal induction, and the variation 
does not show the marked asymmetry about the 
axis J=0 which characterizes the behavior of 
the torsional decrement (cf. Fig. 2 of Part II). 
The minimum value of the hysteresis decrement, 
which occurs at /=0, is about twice that of the 
normal induction. 

The circles and points of Fig. 7 show, respec- 
tively, the normal induction and hysteresis varia- 
tion of the longitudinal decrement in annealed 
Armco iron. It will be noted that here the effect 
of hysteresis is just the reverse of that in the 
unannealed material; the hysteresis values at 
low magnetizations are smaller than the normal 
induction. The eddy current contribution is re- 
sponsible for most of the change in the slope of 
the curve at the highest magnetizations. Cor- 
rected values of A¥ are given in Table I of Part II. 
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The experimental method described in Part I of this 
paper is applied to measure the variation of the rigidity 
modulus and torsional decrement with magnetization in 
unannealed and annealed Armco iron. Formulae are de- 
rived which evaluate the contributions to the moduli and 
decrements arising from eddy currents induced in the 


Ricipiry MopuLUS AND TORSIONAL 
DECREMENT 


N Part I of this paper' W. T. Cooke has 

described a precise experimental method for 
measuring the elastic moduli and coefficients of 
internal friction of solid substances, and has 
presented data which show the observed varia- 
tion of Young’s modulus £, and the longitudinal 
decrement A¥, with magnetization in specimens 
of unannealed and annealed Armco iron. The 
observed variation of the rigidity modulus G and 
the torsional decrement A® with magnetization 
in the same specimens is shown in Figs. 1 to 3. 

Both the longitudinal and torsional decrements 
are subject to correction, for, as was first pointed 
out by Kersten,” the energy dissipated in vibrat- 
ing magnetized substances receives a contribu- 
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Variation of the rigidity modulus with magnetiza- 
tion in unannealed and annealed Armco iron. 


Fic. 1. 


! Cooke, Phys. Rev., preceding paper. 
? Kersten, Zeits. f. tech. Physik 11, 463 (1934). 


vibrating specimen by the change in magnetization pro- 
duced by the varying stress. These formulae contain the 
magnetostriction and Wiedemann effect coefficients, which, 
accordingly, have been measured. All effects of the eddy 
currents are thus eliminated from the final data, which are 
discussed in the light of recent ferromagnetic theory 


tion arising from eddy currents induced by the 
change in magnetization which accompanies the 
vibrational stress in the medium. It is evident 
that the magnitude of this contribution must 
depend on the shape of the specimen. In order, 
therefore, to arrive at a quantity characteristic 
of the material alone the observed values of A’ 
and A®% must be corrected by amounts which, it 
will be shown, depend in part on the longitudinal 
magnetostriction and on the Wiedemann effect, 
respectively. Accordingly these properties of the 
specimen material must be measured. 


MAGNETOSTRICTION AND THE WIEDEMANN 
EFFECT 
Magnetostriction 
The 
made on the ellipsoid which yielded the magnetic 


magnetostriction measurements were 


data reported in Part I.* The “‘form effect" dis- 
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3 The experimental method is described in a forthcoming 
paper by D. Kirkham. 
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Fic. 3. Variation of the torsional decrement with mag- 
netization in annealed Armco iron. 


cussed by Becker‘ is negligible for such an ellip- 
soid. The results are shown in Fig. 4. The pre- 
cision of the measurement of 6/// is +0.1010~-°. 
The unannealed specimen is evidently not 
isotropic, probably in consequence of a non- 
random distribution of internal stresses, or of 
microcrystalline axes. The anisotropy is manifest 
in an abnormally low value of Poisson's ratio, 
computed on the assumption of isotropy; in the 
value of the magnetization at high fields, which 
is larger before than after annealing; and in the 
fact that the Villari reversal does not appear in 
the unannealed material. The last named phe- 
nomenon has been discussed by Bitter.® 


Wiedemann effect 


It has been shown’ that the Wiedemann effect 
is merely a complicated case of magnetostriction : 
the twist produced when a current is passed 
through a magnetized rod is a result of the 
rotation of the principal axes of strain in each 
volume element which occurs when a circular 
field, Hy, is superposed on a longitudinal field, 
H,. This twist may be regarded as due to a 
magnetostrictive shearing stress — n’H», where n’ 
is a function of H7, alone provided H7?«H . For 
a solid cylinder of radius R carrying a longitudinal 
current J e.m.u., the circular field at a distance 


* Becker, Zeits. f. Physik 87, 547 (1933). 

5 Bitter, Phys. Rev. 42, 705 (1932). 

° Fromy, J. de phys. 7, 13 (1926); Ann. d. physique 8, 626 
(1927). 
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r from the axis is 2/r,/ R?. If the cylinder is twisted 
through an angle ¢ radians per cm length, the 
shearing strain is rd, and the resulting stress of 
elastic origin is therefore Grd. The observed 
twist is determined by the condition — n’2/r/R® 
+Gro@=0, whence 7’=GR*¢/2J. It is this 
quantity 7’ that must be known in order to P 
evaluate the contribution of the eddy currents 
to the torsional decrement. It can be determined 

by measurement of the Wiedemann effect at \ 
currents sufficiently low that H?«<H,?; under 
this condition ¢ is proportional to J. 

An accurate measurement of »’ is possible only 
when the specimen is a long thin wire. Pre- 
liminary, calculation indicated, however, that the 
cylinders used in the decrement measurements 
would yield data adequate for the present pur- 
pose. The experimental method is as follows: 


The magnetizing coil assembly described in 
Part I is mounted on a rigid frame with its axis 
vertical. The specimen is clamped at the upper 
end, and to the lower end is attached a brass rod 
which terminates in mercury. The twist of the 
uniformly magnetized portion of the specimen is 
measured by observing the deflection of a light 
beam successively reflected from two mirrors, 
lying in approximately perpendicular planes 
through the axis of the rod, and at the same 
height, but mechanically attached to opposite 
ends of the portion of the rod under test. The 
deflection of the beam is thus unaffected by 
rotation of the system as a whole, and depends 
only on the difference of rotation of the two 
mirrors. To measure the small deflections (a few 
tenths of a mm at 1 m distance) a slit illuminated 
by monochromatic light is used as a source, of 
which a lens between the source and the mirrors 
forms an image about 1 m from the mirrors. The 
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Fic. 4. The longitudinal magnetostriction of unannealed 
and annealed Armco iron. 





MAGNETIZATION AND 
image is viewed through a microscope equipped 
with a 32 mm objective and a filar micrometer 
eyepiece. 

The deflection which accompanies a reversal 
of the current through the specimen is observed 
for different values of the current, the mag- 
netization being held fixed. The initial slopes of 
the current deflection curves thus obtained are 
found by graphical extrapolation. 

Because of the shortness and thickness of the 
cylinders used in the measurements, it is neces- 
sary to consider the effect of torques exerted on 
the ends of the specimen. It has been stated by 
Williams’ that the end of the specimen, being a 
pole in a circular field, is subject to a torque of 
amount 2mJ, where m is the pole strength. This 
torque, if present, would introduce an error not 
greater than 1.5 percent of the maximum ob- 
served value of n’. Williams’ formula, however, is 
obtained by attributing physical existence to 
the fictitious magnetic volume and surface dis- 
tributions p= —V- J and c=n- J. Actually these 
distributions are obtained by a transformation 
of a volume integral, and lead to correct ex- 
pressions for the total force and torque on a 
magnetized body, but not for the force and torque 
on a part of the body. Direct application of the 
formula dF=J-VHdr for the force dF on a 
volume element dr shows at once that the only 
component of force exerting a torque about the 
cylinder axis is 


dFy=(J 0H,/d2+J,0H¢/dr)dr. 


Calculation of the resulting torque on the 
assumption of the worst possible conditions of 
nonuniformity of J and H at the ends of the 
specimen indicates that the error is less than 1 
percent of the maximum value of 7’, and is 
therefore within the precision of the measure- 
ments, which is about 10 percent of this maxi- 
mum value. 

The variation of »’ with intensity of mag- 
netization in the annealed specimen is shown in 
Fig. 5. The relation between twist and current 
in the unannealed specimen is nonlinear, and 
shows much hysteresis. Hence it is possible only 
to set an upper limit of 0.5 10° on the value of 
n’ in this material. 


7 Williams, Phys. Rev. 32, 296 (1911). 
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Fic. 5. The Wiedemann effect in annealed Armco iron. 


THE DissiPATION DUE TO Eppy CURRENTS 


Kersten® has developed a theory which permits 
an evaluation of the contribution to A® arising 
from eddy currents in a medium for which the 
simplifications of Becker’s theory of ferromag- 
netism are valid, and which is in the remanent 
state of magnetization. For the present purpose 
a more complete theoretical analysis is necessary. 

The specimen will be assumed to be an in- 
finitely long, thin cylinder of radius R and 
resistivity p e.m.u., placed in a longitudinal 
magnetizing field /7. Associated with the field 
are an induction B, and a longitudinal mag- 
netostriction e. It will further be assumed that 
the changes resulting from a small variation of 
strain or of field are reversible. Under these 
assumptions, simple thermodynamic consider- 
ations® lead to relations between the alternating 
components of induction and stress, and those of 
field and strain. 

Alternating longitudinal elongations, dw dz, 
and associated lateral contractions (determined 
by the condition that the lateral tension shall 
vanish) produce an alternating induction B,; the 
circular currents induced set up an alternating 
longitudinal field /7,, with consequent alternat- 
ing longitudinal tension Z,. The relations are 


B,=pH,+4rndw/ dz, (1) 


and Z.= —nH + Edw/dz, (2) 


where u is the incremental permeability, dB d//, 


* Cf. Houstoun, Phil Mag. 21, 78 (1911); Becker, Zeits. 
f. Physik 87, 547 (1933); Voigt, Lehrbuch der Kristallphysik 
(1928), pp. 223, 563, 939. 
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at constant strain,? E is Young’s modulus at 
constant field, and n= E(de/dH). 

Alternating torsional strains, due, dz, produce 
a circular induction Bs; the longitudinal and 
radial currents induced set up an alternating 
circular field /7, with consequent torsional stress 
és." The relations are 


By=p'llg+4rn’duX4 02, (3) 


and 6z = — n' II, +Goug, dz, (4) 


where y’ is the permeability, B/ 77 (not dB/dH)," 
and G is the rigidity at constant field. 

If the specimen is undergoing longitudinal dis- 
placements described by the formula 


wn 


W = Wy exp 1(8z+n1), (5) 


the components of magnetic field, magnetic 
induction, electric field and current density are 
determined by the electromagnetic field equa- 
tions, in which the terms may be 
neglected; by the auxiliary relations between 
magnetic field and induction, and between elec- 
tric field and current density; and by the con- 
ditions of finiteness of the field vectors at r=0 
and r=, and of continuity of H7, and B, at 
r=R. The usual relation, B.=ul/., is to be 
replaced by Eq. (1). Solution of these equations 
subject to these boundary conditions gives" 


H ,=4ni(nBw/u) —1+ FIo(hr) J, 


radiation 


(6) 
where 


F=hK,(8R)/[uBl(AR)Ko(BR) +hl(hR)Ky(BR) }, 
h? = Bu we’ +4arnyi Pp, 


and the J’s and K’s are Bessel functions as 
defined in reference 12. At the frequencies of 
vibration here considered the real part of h is 
negligibly small, so that the last equation reduces 
to h=my/i, where m= (4rypn/p)}. 


® The uw obtained from the magnetization curve is the 
incremental permeability at zero stress, and exceeds the 
uw of Eq. (1) by 42n?/E. This difference, as well as the 
differences between the adiabatic and isothermal values of 
u, n, and E, is small enough to be neglected in this compu- 
tation. 

10 The stress and strain notation is that of Love, Treatise 
on the Mathematical Theory of Elasticity, 4th edition, pp. 
56, 90, 288. 

1 A small increment of H perpendicular to the original 
direction merely rotates the H and B vectors through the 
small angle By/B=He/H. 

12 See Gray, Mathews and Macrobert, Bessel Functions, 
2nd edition, Chap. III; Russell, Alternating Currents, 2nd 
edition, Chap. VII. 
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In accordance with Eq. (2) the field H, pro- 
duces a stress whose average value over the 
cross section of the rod is given by the expression 


— (7 -)ay =42 (n/n) 
X(1—(2F/AR)I\(hAR) jdw/dz. (7) 

Like the stress of elastic origin, this contains , 
dw, dz as a factor, but the coefficient of dw/dz is, 
in general, complex. The real part is equivalent 
to a change of Young’s modulus, and measures 
the nondissipative reaction of the currents upon 
the mechanical vibration. The imaginary part is 
equivalent to a viscous coefficient, and measures 
the dissipation of energy by the currents. The 
changes in E and A® thus produced by the 
currents are evaluated by comparing Eq. (2), in 
which the value of —7H, given by Eq. ( 
been substituted, with the usual 
relation in a cylinder with elastic and viscous 


( 
7) has 
stress-strain | 


stresses. Thus 
SE 4r¥' ; 2 O+ (ub a)(P?+0*) 
E Exl  a(1+ubQ/a)?+(ubP/a)*! 
4n’n? 2 P 


Ep a (1+pbQ/a)?+ (ubP, a)? 


and A# 


where 
P = (ber a ber’a+bei a bei’a) /(ber’a+bei’a), 
Q = (ber a bei’a — bei a ber’a) /(ber’a + bei’a), 
a=(f/f.)', b=BRKo(8R)/Ki(8R), 

and f.=p/8uR?°. 


At frequencies far below or above f, the fol- 
lowing simplifications are possible : 


AE 1 x*n* f 
=Q and A,f= “. 
E 2 Ep f. 


f 
1- 12/ 
f 


if 
2( 
f 
where 


c=yub/(av2) =u(f./2f)'BRKo(BR)/K,(8R). 


When f<f., 


When f>f., 
1+cv2 | 


AE 4arn? | ) 
(1te)?+e1" 


E Es| 
and 
4r°n* 


i 1 
) 9) 
(1+c)?+ec 





MAGNETIZATION AND 

Introduction of the values of » and uy at 
remanent magnetization, as given by Becker's 
theory," and of a new reference frequency f,, 
such that f.=(452/512)f,=0.2766f,, leads to 
formulae which agree with Kersten’s'™ provided 
c<K1. 

If the specimen is undergoing torsional dis- 
placements described by the formula 


Ug = 1d = rho Exp 1(8’s+n1), 


the quantities to be evaluated instead of /7, and 
(Z.)ay are 7, and the total torque 


R 
anf 6zr2dr 
0 


exerted across a surface normal to the cylinder 
axis. The analysis parallels that for longitudinal 
vibrations and leads to the following formulae: 
AG 4m” | 4 
= 1—-P’ 
G Gu’ | a |’ 


and 


where 
P’ = (ber a ber’a +bei a bei’a) /(ber”a+bei”a), 


Q’ = (ber a bei’a — bei a ber’a) /(ber”a+bei’a), 


a=(f/f.’)' and f.’=p/8x*u'R?®. 
AG ie f 
When I<. # =() and A& =. ; 
G 6 Gy py 
When f>f.’, 
AG 4 | f. | 
= 1 -22( ) ; (10) 
G Gy’ | fl | 
ry” e\} 
and 4.6 =— 2v2( ) (11) 
TM f 


APPLICATION OF THE THEORY 


In all the measurements under discussion the 
formulae for f>/f. are valid. Since the second 
terms in Eqs. (8) and (10) are small compared 


18 Kersten, Zeits. f. Physik 71, 553 (1931); 82, 723 (1933). 

4 Kersten, reference 2, Eqs. (3) and (4). 

The terms involving c result from the finite wave- 
length prescribed by the mechanical oscillation. Kersten’s 
analysis assumes this wave-length to be infinite. 
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ELASTICITY IN IRON 
with the first, the measured values of the elastic 
moduli are very nearly those at constant B 
rather than at constant /7. For it follows im- 
mediately, on setting B, and By equal to zero in 
Eqs. (1), (2) and (3), (4), that 


[Ee - [E Ju =49n?/p 


(G)s—(G)a=4rn"/x’, 


and 


and these are precisely the differences given by 
the first terms of Eqs. (8) and (10). As was 
pointed out by Kersten, this means that the 
alternations of magnetic flux occur only near the 
surface of the specimen. 

The effect of the eddy currents on the torsional 
decrement is small, being nowhere more than one 
percent of the measured value for either of the 
specimen materials. 

In the measurements of the longitudinal decre- 
ment for the unannealed material, two cylinders 
of different diameter were used (cf. Part I, Fig. 6). 
The agreement between the results obtained on 
these two specimens at low magnetization indi- 
cates that here the eddy current contribution is 
negligible, since, by Eq. (9), it should vary 
approximately in the 
diameter. The second rise in the curves at high 
magnetization, however, is entirely due to eddy 
currents, and, together with the difference in 
behavior of the two specimens, disappears when 
the contribution of the eddy currents to the 
decrement is subtracted. Similarly, the eddy 


inverse proportion to 


currents are responsible for the change of slope 
observable at the highest magnetizations in the 
curve of A® vs. J for the annealed specimen (Part 
I, Fig. 7). The corrected values of A*® for this 
substance are given below in Table I. 

It should be noticed that the formula, Eq. (9), 
involves two factors, 7 and yw, which must be 
evaluated by graphical differentiation of experi- 
mental curves, and which, moreover, are raised 
to powers higher than the first. It must be 
remembered also that the specimen is not actu- 
ally infinitely long and not uniformly magnetized. 
Consequently the correction for eddy currents 
cannot be computed with great precision; but 
the success of the theory in explaining the second 
rise of the curves in Fig. 6, and in removing the 


‘discrepancies between the observations on dif- 


ferent specimens at high magnetization, leaves 
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TABLE I. Variation of the decrements and elastic moduli at 
constant B with magnetization in annealed Armco 
iron. Corrected for the effect of eddy currents. 

Normal magnetization. 


Jgauss|| A” x10 (AE/E) X 102 Aa® x 108 AG G) X102 
0 25.91 0 25.75 0 
40 25.90 25.79 —0.0020 
80 26.13 25.88 —0.0023 
119 26.30 26.06 —0.0017 
159 26.28 —().0037 25.97 —0.0003 
199 26.45 25.84 +0.0003 
239 26.75 25.66 0.0005 
279 || 26.30 25.44 0.0013 
318 || 26.34 +0.0058 | 25.26 0.0103 
358 | 24.94 
398 25.88 24.76 0.0187 
477 25.83 0.0105 24.09 0.0244 
557 25.19 23.78 0.0307 
637 24.55 0.0289 23.55 0.0366 
716 23.98 23.51 0.0446 
796 23.19 0.0360 23.51 0.0481 
875 22.11 23.51 0.0546 
955 21.43 0.0556 23.46 0.0595 
1034 20.43 23.51 0.0650 
1113 19.13 0.0779 23.55 0.0700 
1193 || 17.09 0.1041 23.60 0.0776 
1233 15.48 23.64 
1269 || 14.27 0.1532 23.45 0.0938 
1342 11.60 0.1926 22.94 0.1167 
1376 || |} 22.34 | 0.1295 
1410 || = 8.81 0.2323 21.73 0.1412 
1444 || 20.98 0.1571 
1475 | 6.84 0.2625 | 20.01 0.1743 
1506 18.75 0.1914 
1537 || 5.30 0.2910 17.27 0.2106 
1565 15.66 0.2257 
1593 || 4.14 13.95 0.2421 
1619 |) 11.98 0.2582 
1623 || 3.65 0.3077 11.66 


little doubt that it is adequate for the present 
purpose, and that the corrected values of the 
decrement are characteristic of the material 
itself. 

Correction of the observed variations of E and 
G to obtain the values at constant B is easily 
accomplished, since the correction terms are the 
second in Eqs. (8) and (10) and these are small. 
The correction to constant #7 includes the first 
terms in these formulae and is less reliable. 
[AE/E]q for the unannealed specimen follows 
the observed curve (Part I, Fig. 5) to /= 1400; 
beyond this point the corrected curve flattens 
out, reaching a maximum of 0.12 percent at 
about 1490 gauss and decreasing to 0.11 percent 
at 1630 gauss. For the annealed specimen, the 
flattening out begins at /=1250, and the curve 
rises more slowly, reaching about 0.17 percent 
at 1630 gauss. The effect of eddy currents on the 
rigidity is negligible for the unannealed specimen ; 
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for the annealed specimen it is negligible except 
between J=1100 gauss and J=1500 gauss, 
where [ AG /G ]y is lower than the observed value, 
differing therefrom by perhaps 0.025 percent at 
J = 1300 gauss. 

The corrected values of the decrements and of 
the changes in the elastic moduli at constant B, 
for the annealed material in normal magnetiza- 
tion are given in Table I. 


THEORETICAL DISCUSSION 


According to the domain theory of ferro- 
magnetism,'® a ferromagnetic crystal consists of 
a large number of ‘‘domains,’’ each of which is 
permanently magnetized to saturation; the mag- 
netization process for the crystal as a whole 
consists of a reorientation of the magnetization 
vectors in the domains. At high fields all the 
domains are magnetized in the field direction ; in 
the demagnetized state their magnetization 
vectors are distributed at random among the 
various stable directions determined by crystal- 
line forces and internal stresses. The transition 
from the demagnetized condition to technical 
saturation consists of two stages. In the steep 
part of the magnetization curve the random 
distribution changes to one in which most of the 
domains are magnetized along the stable direc- 
tion nearest that of the field—this direction itself 
having meanwhile rotated toward the field. 
Above the knee of the magnetization curve the 
magnetization vectors are gradually rotated into 
the direction of the field. In the first stage, the 
redistribution occurs partly by a continuous and 
reversible growth of the more favorably oriented 
domains at the expense of the less favorably 
oriented, and partly by a sudden transition of 
the domain magnetization from a less stable to 
a more stable orientation, with dissipation of 
energy. 

Stress plays a role similar to that of magnetic 
field in determining the possible directions of 
magnetization and the relative stability of these 
directions; and conversely, any change in the 
direction of magnetization of a domain (except 
a reversal) is accompanied by magnetostrictive 
changes in its dimensions. The application of a 


16 See especially Bozorth and Dillinger, Phys. Rev. 41, 
353 (1932); Becker, Physik. Zeits. 33, 909 (1932). 
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tension or shearing stress to an unmagnetized 
ferromagnetic therefore produces redistributions 
and rotations similar to those produced by the 
field, and these give rise to a magnetostrictive 
component of strain in addition to the strain of 
purely elastic origin. At very high fields this is 
no longer true, because the direction of mag- 
netization is that of the field and is unaffected 
by small stresses. There are two observable 
results of this magneto-mechanical process, both 
evident in the data on the elastic constants and 
decrements of annealed iron. First, the ratio of 
strain to corresponding stress is greater in the 
demagnetized state, and, in fact, throughout the 
steep part of the magnetization curve, than at 
high fields ; and second, the energy dissipated per 
cycle of vibration is greater at low than at high 
magnetizations."” 

It remains to explain why the internal stresses 
present in the unannealed material, which in- 
crease the magnetic hysteresis, nevertheless 
decrease the internal friction at zero magnetiza- 
tion. This is not difficult. For there is this 
important difference between the effect of field 
and stress: the former brings about a difference 
of stability between orientations 180 degrees 
apart and the latter does not. In relatively stress- 
free iron the stable directions are the six (positive 
and negative) crystal axes, and stress and field 
alike can reorient through 90 degrees, partly 
reversibly and partly irreversibly. In unannealed 
material, on the other hand, the stable directions 
are determined by the internal stresses, and are 
two in number, differing by 180 degrees. Transi- 
tions produced by the field are possible and are 
attended by considerable loss of energy, so that 
magnetic hysteresis is large. But transitions 
produced by stress cannot occur, and the only 
magneto-mechanical process which contributes 
to the AE and AG effects is a reversible rotation. 
As the field is increased, however, opposite 
orientations begin to differ in energy sufficiently 
that a small stress, acting on a domain mag- 
netized in the less stable direction, can make this 
orientation unstable and produce an irreversible 
transition. Thus at fields sufficiently large to 
cause considerable difference in stability, yet not 


'7 Cf. also the experiments of Becker and Kornetzki, on 
torsional stress cycles in annealed iron, Zeits. f. Physik 88, 
634 (1934). 
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large enough to complete the transition process, 
a high mechanical decrement is observed. 

A more detailed explanation of the variation 
of ferromagnetic internal friction with magnet- 
ization must await a more complete analysis of 
the magnetization process. In connection with 
the AE and AG effects, however, certain calcu- 
lations are of interest. 

Kersten’s'* theory of the AE effect for material 
under high internal stress is not applicable to 
iron, but Akulov’s” theory for iron under low 
internal stress should be applicable to the 
annealed material. The theory contains formulae 
which describe the variation of E with magnet- 
ization at low magnetizations, and for the total 
change between demagnetization and saturation. 
The observed variation at low magnetizations is 
too small to yield any information, and even the 
highest fields used in this work are still consider- 
ably below that necessary for saturation. Never- 
theless a comparison between the greatest 
observed change in £ and G and the total change 
predicted by the theory should at least indicate 
whether the theory is correct as to order of 
magnitude. 

The total changes in the elastic moduli 
between demagnetization and saturation in 
quasi-isotropic polycrystalline iron are given by 
the formulae 

AE/E=?xar00°E/J,., (12) 
and AG 'G= (9/5) xoA100°G/J...”, (13) 
where /,, is the saturation magnetization, Ajo 
the longitudinal magnetostriction for a single 
crystal saturated in the (100) direction, and x» 
the observed initial susceptibility of the material. 
The first of these formulae appears in Akulov’s 
paper. The second can be derived therefrom in 
the following manner: It is assumed that the 
compressibility is independent of magnetization 
(since no magneto-mechanical process follows the 
application of a uniform pressure), and in both 
the demagnetized and saturated state the 
material is isotropic. If the compressibility of an 
isotropic substance is to remain constant then 

18 Kersten, Zeits. f. Physik 85, 708 (1933). 


19 Akulov and Kondorsky, Zeits. f. Physik 78, 801 
(1932); 85, 661 (1933). 
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variations in E and G must be related by the 
formula AG/G=(3G/E)AE/E. 

Eq. (12) has proved only moderately successful 
in accounting for the AE effect in polycrystalline 
nickel.”® In this connection it should be noted 
that it has never been possible to obtain entirely 
satisfactory agreement between the measured 
values of the elastic constants of a polycrystalline 
material and those calculated from the properties 
of the constituent crystals. This problem has 
been thoroughly discussed by Bruggeman.”' In 
Akulov’s theory, uniform stress is assumed 
throughout the aggregate of crystals and the 
strains are averaged. The alternative method, 
originally used by Voigt,” is to assume uniform 
strain and to average the stresses. If this method 
of averaging is used in Akulov’s theory, the factor 
1/5 in Eq. (12) is replaced by (1/20) (¢1: —¢i2)?/ G’, 
TABLE II. Comparison of theoretical and experimental values 


of the elastic moduli and of their variation 
with magnetization. 


AE ee | 9x10 
> “| > A » . 
EX10-"| GX10""' | kX 10 E G 
Uniform 
stress 19.45 7.42 17.3 0.309 0.393 
Theoretical 
Uniform 
strain 22.80 8.90 17.3 0.101 0.128 
Experimental 19.98 846 | 104 0.17 0.26 


20 Siegel and Quimby, Phys. Rev. 49, 663 (1936). 
*1 Bruggeman, Zeits. f. Physik 92, 561 (1934). 
% Voigt, Lehrbuch der Kristallphystk, p. 962. 





WILLIAM FULL 


E 


R BROWN, JR. 


where the c’s are the elastic constants of the 
crystals. 

The experimental values of Young’s modulus 
and the rigidity, and of their changes with mag- 
netization up to the highest magnetizations 
attained, are compared with the values predicted 
by the theory on the alternative assumptions of 
uniform stress and uniform strain, in Table II. 
The following values of the constants were used 
in the calculation: J,,=1710; Aioo = 16.2 K 107° 8 
Xo = 28.73% =, = 2.37 X10",  cye = 1.4110", 
C= 1.162X10".25 The bulk modulus & is in- 
cluded to show the unreliability of the methods 
of averaging even when the question of uniform 
stress or uniform strain is not involved. The 
agreement between theory and experiment 
affords as good a confirmation of the theory as 
can be expected from measurements on poly- 
crystalline materials. The desirability of meas- 
urements on single crystal specimens is obvious. 

In conclusion the writers wish to acknowledge 
their indebtedness to Dr. S. L. Quimby, who 
suggested the research and followed its progress 
with helpful counsel and encouragement, and to 
those who at various times assisted them in the 
experimental work. 


23 Average value from measurements of Honda and 
Mashiyama, Sci. Rep. Tohoku Imp. Univ. 15, 755 (1926). 

* Part I of this paper. 

* From measurements of Goens and Schmid, Naturwiss. 
19, 520 (1931). 
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REVIEW 


Ferromagnetic Anisotropy of Iron Crystals at Various Temperatures* 


Raymonp G. Piety, Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
(Received August 14, 1936) 


Magnetization curves in important crystallographic directions in small oblate spheroids 
cut from hydrogen-purified iron crystals show ferromagnetic anisotropy of the same sort as 
that previously reported for iron of less purity. At room temperature agreement with earlier 
work is good. At 296°C and at 500°C the principal (fourth-order) term in the anisotropy is 
larger, at 650°C and 696°C the same term is smaller than in earlier studies, anisotropy not being 
detected. The subsidiary (sixth-order) term is not accurately determined. The pendulum 
magnetometer method, which has certain advantages for such purposes over classical methods, 
was used and results for different spheroids agree well, especially at the higher temperatures. 





INTRODUCTION 


HE ferromagnetic anisotropy of iron crystals 

has been measured by Beck,' Webster, 
Honda and Kaya,’ Honda, Masumoto and 
Kaya.* The last of these studies covered the 
range in temperature from 5°C to 770°C, the 
others were confined to room temperature. Beck 
and Webster measured the torque on disks 
suspended in the air-gap of an electromagnet, 
Honda and Kaya measured the flux change 
through a search coil at the center of which an 
oblate spheroid was suddenly turned from one 
azimuth to another in the field of a magnetizing 
solenoid. Honda, Masumoto and Kaya used one 
of the same spheroids but changed the method of 
measurement by leaving the specimen in one 
azimuth during each cycle of heating and cooling, 
getting /—H curves at selected temperatures by 
reversing various currents in the solenoid. 

All the data agree in fixing the tetragonal axes, 
of form <100>, as the directions for easiest 
magnetization, the trigonal axes, of form <111>, 
as the directions for most difficult magnetization. 
The digonal axes, of form <110>, are inter- 
mediate between these extremes. The quantita- 
tive differences between results from different 
laboratories are larger than the errors in measure- 
ment can easily explain, so that they must be 
ascribed in part to differences in chemical purity, 








* Based on part of a dissertation presented to the 
Faculty of the Graduate School of Yale University in 
candidacy for the degree of Doctor of Philosophy. 

1K. Beck, Dissertation (Ziirich, 1918); Vierteljahrsschrift 
der naturforschenden Gesellschaft in Ziirich 63, 116-186 
(1918). 

2 W. L. Webster, Proc. Roy. Soc. A107, 496-509 (1925). 

3K. Honda, S. Kaya, Sci. Rep. Tohoku Imp. Univ. [1] 
15, 755-770 (1926). 

*K. Honda, H. Masumoto, S. Kaya, Sci. Rep. Tohoku 
Imp. Univ. [1] 17, 111-130 (1928). 


in mechanical and thermal treatment. Under 
these conditions it seemed worth while to re- 
measure the anisotropy in very pure and nearly 
unstrained iron by a method capable of fair pre- 
cision and relatively free from systematic errors. 
Such a method is that developed in this labora- 
tory® and first applied here by J. W. Shih*® to 
iron-cobalt alloys. It uses an accurately made 
oblate spheroid carried on a light pendulum at 
the common center of a magnetizing solenoid and 
a pair of coils producing a weak but spatially uni- 
form magnetic field gradient. These various 
windings have a common horizontal axis so that 
the resultant force on the specimen is also hori- 
zontal and can be measured by the deflection of 
the pendulum. The computed value of J then 
depends, except for a factor constant for each 
specimen, only upon the measured deflection of 
the pendulum and the measured current in the 
gradient coils. The absolute value of H within 
the specimen depends upon the measured current 
in the magnetizing solenoid and upon the value 
of J, but the difference between values of // along 
different diameters of the spheroid for the same 
value of J does not depend upon J, and it is this 
difference which is most important in the analysis 
of the results. Points on J—H curves near J: are 
particularly important in anisotropy analyses. 
These points are determined with greatest pre- 
cision by the new method, with least precision by 
the older methods. 

Several improvements in technique, recently 
described’ were devised and tested during the 
progress of the work reported herein. 

5L. W. McKeehan, Rev. Sci. Inst. 5, 265-268 (1934). 

*]. W. Shih, Phys. Rev. 46, 139-142 (1934). 


*L. W. McKeehan, R. G. Piety, J. D. Kleis, Rev. Sci. 
Inst. 7, 494-497 (1936). 
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SPECIMENS 


Comparison of properties in different directions 
is safest when as many as possible of the impor- 
tant directions are accessible in one specimen. 
This condition is best met, as is well known, in an 
oblate spheroid with its equatorial plane of the 
crystallographic form {110}. If we choose to call 
this the (011) plane we have lying in it the follow- 
ing important axes: [100], [011], [111], [111]. 
Principal effort was therefore expended on 
spheroids of this sort, (011) spheroids. As a check 
a (001) spheroid was also cut, containing the 
important axes: [100], [010], [110], [110]. 

The iron was very kindly furnished by the 
Bell Telephone Laboratories in three flat pieces 
about 3 mm thick and several centimeters in 
diameter. These were parts of boxes in which P. 
P. Cioffi had purified small iron specimens by 
washing them with pure hydrogen at about 
1450°C for 18 hours or more. His studies on this 
process,* which produces extreme magnetic and 
mechanical softness, lead to the expectation that 
the only nonmetallic impurities remaining in 
these boxes in measurable amounts are nitrogen, 
sulphur and phosphorus. The first three of these 
together should not exceed 0.01 weight percent. 
Silicon may be in excess of this limit, especially in 
those parts of the box which are in contact with 
the alundum support, but these parts were not 
used here. A provisional limit of 0.03 weight 
percent silicon is suggested. Manganese and 
copper may amount to about 0.02 weight percent 
each. It can be stated with confidence that this 
iron is more nearly pure Fe than any previously 
used in anisotropy studies. 

The crystals present in the pieces as furnished 
were of ample size for the pendulum magnetom- 
eter. The surface was first etched rather deeply 
with a mixture of water solutions of HNO; and 
HCl, restrained with C,H;OH. By observing re- 
flections from etch planes—here of form {100}— 
crystals of suitable size and orientation were 
selected. These were cut out by hand with a 
jeweler’s saw and mounted on an optical goni- 
ometer of special design® for more accurate ob- 
servation. The orientation of the <100> axes 





§P. P. Cioffi, Phys. Rev. 39, 363-367 (1932): 45, 742 
(1934). 

*L. W. McKeehan, H. J. Hoge, Zeits. f. Krist. 92, 476- 
478 (1935). 
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being thus determined with respect to the surface 
and edges of the fragment a rectangular iron bar 
was cut off at such an angle that when the frag- 
ment was cemented (with chaser’s cement) on 
the cut in the proper azimuth, the axis of the 
spheroid to be cut was parallel to the axis of the 
bar. The blank was now faced off by extremely 
light slow cuts on a Brown and Sharpe surface 
grinder with a diamond-dressed wheel running at 
high speed. It was remounted on the squared end 
of a new holder and recut in the same way until 
the thickness was less than 0.5 mm but still 
somewhat greater than that of the proposed 
spheroid. It was now ready for an x-ray inspec- 
tion of its orientation and freedom from distor- 
tion by the Laue method, using transmitted 
radiation so as to get a volume average. If the 
plane of the blank was found to be more than 
about 1° from the proposed crystallographic 
plane a new holder was cut and both sides were 
reground, otherwise only one side was reground. 
In either case the thickness was brought down to 
0.006 cm over the proposed minor axis of the 
spheroid and the blank was filed, on its holder, to 
a circle a little larger than the proposed equatorial 
section. From this point on, the procedure has 
been described in the previous paper.’ 

The spheroids were finally given a second x-ray 
inspection by the Laue method to make sure that 
no serious distortion had resulted in the last 
stages of shaping. The criterion was that the 
Laue spots should be symmetrical ellipses with- 
out more than traces of radial asterism. 

The applied field intensity, H+N/, is very 
much greater than the effective field intensity // 
in an oblate spheroid of practicable dimension- 
ratio (diameter/thickness). It therefore becomes 
very important to prove that the demagnetizing 
field, NJ, is reasonably uniform in intensity 
throughout the volume, as the method of meas- 
urement supposes. Partly to test this and partly 
for other more obvious reasons two spheroids 
having the same orientation but differing in 
eccentricity were tested over the whole available 
range in temperature. The equatorial plane in 
each of these lay within the accuracy of the x-ray 
analysis (about 1°) parallel to the (011) plane. 
Their characteristics are given in Table I. Under 
Novs is given the demagnetizing factor neces- 
sary to make the effective field sensibly zero up 
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TABLE I. Data on spheroids (011). 


Dimen Major | Minor 
Desig- sion Axis Axis | Volume; Mass 
nation| Ratio cm cm) |10-%cm*; (g) | Nobs | Neale 


P10 | 10/1 | 0.318 | 0.0318) 1.768 | 0.0139) 0.912 | 0.875 


P12 | 12/1 | 0.307 | 0.0256) 1.234 | 0.0097 0.710 | 0.743 


to half-saturation in the actual measurements 
of J and H. This differs appreciably from Neate, 
the demagnetizing factor for a perfect spheroid of 
the measured dimensions, but since N,,, was not 
dependent upon the direction along which H was 
applied it seems reasonable to suppose that the 
slight departures from ellipticity thus indicated 
are symmetrical with respect to the minor axis 
of each spheroid and that any resultant inhomo- 
geneity in // and J throughout the volume is not 
so great as to make the measured averages mis- 
leading. It should be noticed that N,., was 
greater than N..). in one case, less in the other. 
A third spheroid cut with its equatorial plane 
parallel to (001) was tested less completely, but 
showed directional independence of N,),, and fair 
agreement between N,», and Neate. 

Temperatures higher than 700°C could not, 
with available means, be kept constant enough 
to make /—// curves precise enough for anisotropy 
work. The difficulty is that temperature fluctua- 
tions change /,,, by relatively great amounts at 
temperatures near the Curie point so that it 
becomes very hard to take even one complete 
1-H curve under sufficiently constant conditions. 
So far as some unsatisfactory experiments indi- 
cated, however, there was no detectable ani- 
sotropy in magnetization above 700°C in this 
iron. 

At 296°C and at 500°C the differences between 
I-H curves for the two (011) spheroids were so 
small that single curves were drawn to fit the 
combined data. At room temperature, 18°C, the 
data were treated separately. 

Fig. 1 shows representative J—/J7 curves and 
Table II gives areas included between such 
curves and the two lines 7=/,,, and H=0, to- 
gether with derived quantities. All curves were 
carried to saturation, applied fields as high as 
1800 oersted being used if necessary. Some of the 
areas differ slightly from those published in my 
dissertation. This occurs because all the areas 


have been redetermined by a more precise method. 
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TABLE II. Work needed to magnetize. 


Speci- 104 erg. cm™ 

mens Win Wino Wio=Ko Ky Ky 
P10 18 16.10 13.16 2.48 +42.7 —17 
P12 18 16.59 12.38 1.33 +44.2 +14 
P10, 12 296 8.08 6.45 0.65 +23.2 — 8 
P10,12 500 2.56 0.21 +94 121 
P12 675 : 0.19 0 0 
P12 696 0.31 0 0 


ANALYSIS OF RESULTS 


In the manner proposed by Gans " we suppose 
that the work needed for isothermal magnetiza- 
tion to saturation, measured by the area defined 
above, is for any arbitrary direction in an uncon- 
strained cubic crystal given by 


W. = Kot+ K(a*a;? +a va; + a;"a9") 
+ Ko(a;*as*a;"). (1) 


In this equation a, a, a; are the direction 
cosines of the chosen direction referred to the 
three <100> axes, Ko, K,; and Ky are not direc- 
tion dependent and the only assumption made is 
that an equation of the sixth degree in the cosines 
is adequate."' The symmetry of the crystal annuls 
all other combinations of a, ae, a3. For the 
important directions of form <100>, <110> 
and <111> Eq. (1) reduces to 


Wioo=Ko; Wii=Kot+Ki 4; Win = Ko 
+Ki/3+K2/27. (2) 


These equations are not new. Some authors” use 
instead of K, a constant which in this notation 
is 2K. 

It is also possible, as Akulov™ has shown, to 
compute /—// curves for particular values of K,, 
assuming K, to be zero. As might be expected the 
computed curves agree best with experiment at 
room temperature and are not good approxima- 
tions at high temperatures. 

10 R. Gans, Physik. Zeits. 33, 924-928 (1932). 

1 The analysis which follows differs in one important 
respect from that presented in my dissertation. R. M. 
Bozorth has pointed out in private discussion that the 
sixth degree term in Eq. (1) is capable of explaining aniso- 
tropy hitherto regarded as completely anomalous. I have 
therefore used his method of computing the coefficient of 
this term, and have included it in Tables II and III. 
I am much indebted to him for the opportunity of using 
his method of analysis prior to its publication. 

2 Especially N. S. Akulov, Zeits. f. Physik 57, 249-256 
(1929); 67, 794-807; 69, 78-99 (1931). 
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Fic. 1. Magnetization curves for specimen P12 at three temperatures, 18°C, 500°C and 675°C. In the last 
case the J—(H-+ NIJ) curve is also shown (in part) to illustrate the importance of N. 





ANISOTROPY OF 


It will be noticed that in Table II the constant 
Kz varies widely in absolute value and even in 
sign. As far as these data go the necessity for the 
sixth degree correction is apparent only in the 
fact that Wi;; and Wi, are indistinguishable at 
and above 500°. (Precise equality requires 
K,= —9K,/4.) If the values of K, are plotted 
against temperature they fall on a straight line 
which intersects the line K,;,=0 at 630°C. That 
K, does not actually become negative, however, 
is suggested by the tabulated values for 675° and 
696°C. No permanent effect of heating and cool- 
ing the specimens during the measurement was 
detected at room temperature. 

The the 
Masumoto and Kaya‘ is very good at 


Honda, 


room 


agreement with work of 
temperature, not so good at the higher tempera- 
tures. For comparison, values of Ko, K, and Ke 
have been computed from their tables. In making 
this computation it has been assumed as in the 
present work that the correct demagnetizing 
factor to be used in getting // is that which makes 
the J-H curve rise vertically from the origin. 
This makes it necessary to decrease the areas 
derived from their tabulated HT and J by an 
amount not exceeding 0.5 < 104 erg. cm~* in any 
case, and usually much less. Where their max- 
imum value of / does not occur at their maximum 
value of Hit may be concluded that J has been 
over-corrected for flux through the air surround- 
ing the specimen and the final drop in J has then 
been ignored in computing the area. Finally, 
small differences in J... have been corrected for 
by multiplying the computed areas by appropri- 
ate factors nearly equal to unity. The result of all 
this is to give a series of areas from which Ko, K, 
and Ky may be computed as already explained. 
Table III presents the results. The values of K, 
and Ko, especially the latter, have been made to 


progress more smoothly up to 700°C by using 
faired curves for Wioo, Wiio and W,,,; to find the 
differences used in computation. Our values for 
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TABLE III. Work needed to magnetize. Derived from data of 
Honda, Masumoto and Kaya. 


10‘ erg. cm™ 

Cc Ko K, K, 

5 0.60 +41.3 +14 
100 0.27 +35.3 + 20 
200 0.26 +27.0 +22 
300 0.38 +18.9 +19 
400 0.015 +11.4 +11 
500 0.014 + 5.4 + 9 
600 0.065 + 2.6 — | 
700 0.14 + 1.0 — 2 
750 0.21 + 1.7 —12 
770 0.33 + 0.8 — 4 


K, at 296°C and 500°C are considerably higher 
than those at 300°C and 500°C in Table III. It 
will be noticed that Ko seems to rise a very little 
as the Curie point is approached in both sets of 
experiments. No one would be likely to suggest a 
linear drop in A, with increasing temperature 
from the values given in Table III. An exponen- 
tial drop like that nickel by 
Akulov™ in a very recent paper is a fair ap- 


suggested for 


proximation above 200°, but not below. 

In conclusion it can be stated that the purest 
iron we could procure is both qualitatively and 
quantitatively much like the iron used by Honda 
and his co-workers. It seems, however, to retain 
its anisotropy better at 300°C and 500°C and to 
lose it more suddenly in the neighborhood of 
600°C. The sixth-order term in the anisotropy 
equation is relatively unimportant except at 
500°C the fourth- 
order term falls to negligible values. 

The author wishes to thank Professor L. W. 
McKeehan for suggesting the problem, for advice 
and encouragement during the progress of the 
work, and for analyzing the tables of Honda, 
Masumoto and Kaya. He also thanks J. D. Kleis 
for valuable assistance in manipulations and 


temperatures above where 


measurements. 


18 N.S. Akulov, Zeits. f. Physik 100, 197-202 (1936). 
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Ferromagnetic Anisotropy of Nickel-Iron Crystals at Various Temperatures* 


J. D. Kvets, Sloane Physics Laboratory, Yale University, New Haven, Connecticut 


(Received August 14, 1936 


Magnetization curves in important crystallographic directions at temperatures up to near 
the Curie points in small oblate spheroids cut from annealed face-centered cubic nickel-iron 
alloys with from 10 to 50 percent iron show a change in ferromagnetic anisotropy at about 24 
percent iron. With less than this amount of iron the order of increasing difficulty for magnetiza- 
tion to saturation is, as in nickel at room temperature: <111>, <110>, <100>. With more 
than this amount of iron the order is <100>, <111>, <110>. A small decrease in the ap 
parent demagnetizing factor of the spheroids with rising temperature was noted. 


TUDIES of ferromagnetic anisotropy in solid 
solution alloys have already shown that the 
direction for easiest magnetization may change 
from one important crystallographic axis to 
another at room temperature as the composition 
changes. The first case of this sort to be reported 
occurs in the face-centered cubic alloys of nickel 
and iron. Lichtenberger! examined single crystal 
rods with their lengths parallel to a number of 
different directions in the crystal lattice. The 
number of specimens available in each composi- 
tion was not large and the cross sections were not 
strictly uniform, so that the data scatter rather 
badly. He that the 
easiest magnetization changed from the form 
<111> to the form <100> as the iron content 
increased through 29 percent. 
The second and third cases have been reported 
by Shih.? In body-centered iron-cobalt alloys, 
which he studied in this laboratory, he found the 


concluded directions for 


direction for easiest magnetization changed from 
the form <100> to the form <111> as the 
cobalt content increased through 42 percent. In 
face-centered nickel-cobalt alloys he found two 
transitions, from <111> to <100> at about 5 
percent cobalt, from <100> to <111> at about 
19.5 percent cobalt. 
Lichtenberger’s result has received partial 
confirmation from work by Burgers and Snoek,* 
who conclude from a change in the ferromagnetic 


* Based on part of a dissertation presented to the 
Faculty of the Graduate School of Yale University in 
candidacy for the degree of Doctor of Philosophy. 

1F, Lichtenberger, Ann. d. Physik 15, 45-71 (1932). 

2J. W. Shih, Phys. Rev. 46, 139-142 (1934); 50, 376- 
379 (1936). 

3 W. G. Burgers and J. L. Snoek, Zeits. f. Metallkunde 
27, 158-160 (1935); J. L. Snoek, Nature 137, 423 (1936). 


anisotropy of rolled polycrystalline sheets that 
the critical composition in the nickel-iron series 
occurs at 34 percent iron. In view of the uncer- 
tainties already noted in Lichtenberger’s work a 
more careful examination of a few alloys in the 
critical range was undertaken here before the 
notes by Burgers and Snoek appeared and was 
continued with even greater interest thereafter. 

The ferromagnetic anisotropy of nickel has 
been studied over a wide range in temperature by 
Honda, Masumoto and Shirakawa.* Although 
they did not comment upon the fact, their tables 
and curves show that the direction for easiest 
magnetization changes in this face-centered cubic 
metal from <111> at temperature to 
<100> at temperatures of 200°C and higher. 
The technique of measurement with the pen- 
dulum magnetometer, which has been described 
elsewhere,’ made it relatively easy to extend the 
range of investigation on nickel-iron alloys from 


room 


room temperature nearly to the Curie points, 
which in this series are not higher than 570°C. 


SPECIMENS 


Nickel-iron alloys of selected compositions 
were kindly furnished by the Bell Telephone 
Laboratories in the form of small rods. The 
original castings, rods about 2 cm in diameter, 
had been homogenized at high temperatures and 
rendered still more uniform in composition in the 
process of working. The pieces furnished were 
similar in quality to those used by Buckley, 


*K. Honda, H. Masumoto and Y. Shirakawa, Sci. Rep. 
Tohoku Imp. Univ. [1] 24, 391-410 (1935). 

5L. W. McKeehan, Rev. Sci. Inst. 5, 265-268 (1934); 
L. W. McKeehan, R. G. Piety, J. D. Kleis, Rev. Sci. Inst. 
7, 494-497 (1936). 
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TABLE I. Data on spheroids (011). a is angle between 
rotation axis of spheroid and nearest axis of form <100>; 
8 ts angle between rotation axis and the plane (100); N ts 
the demagnetizing factor necessary to make initial portion of 
I-H curves at room temperature coincide with the line H =0. 


| } | 
Volume 
cm! |a°® g° N 


Diam- | Thick- 
eter ness 
cm cm 


Desig-| Iron | Nickel 
nation | percent* | percent* 














K10| 10 90 | 0.3128) 0.0284/0.001462/44 2 | 0.818 
K30| 30 70 
K35| 35 65 


KS0 50 50 


| 0.3104) 0.0259|0.001304/40 0 | 0.84 
| 0.3090} 0.0256/0.001295/45 4 | 0.761 
0.3070) 0.0276'0.001350)43 2 | 0.82 





* By synthesis. Analysis for nickel and iron requires more material 
than could be spared. Castings made in the same way, however, always 
gave nickel and iron within one percent of the intended values 


McKeehan and Cioffi® in studies of magnetiza- 
tion and magnetostriction in the nickel-iron 
alloys in the permalloy range. Impurities, by 
analysis, are principally sulphur, silicon and co- 
balt, amounting to about 0.15 percent. 

Crystals a few millimeters in diameter and a 
few centimeters long were grown in magnesium 
oxide crucibles by the slow solidification process 
using a furnace similar to one described by 
Quimby.’ After breaking away the crucible as 
carefully as possible the castings were annealed in 
hydrogen at atmospheric pressure at about 
1100°C for one hour and slowly cooled. Laue 
photographs show but a small amount of residual 
lattice distortion. 

The method of selecting crystals for study was 
similar to that described by R. G. Piety® in con- 
nection with iron crystals. The only important 
difference lay in the occurrence of etch planes of 
forms {111} and {100} instead of {100}! alone. 
Etching is more difficult, also, agua regia being 
necessary for the more resistant compositions. 
After selection the process of forming an oblate 
spheroid in the desired orientation was that al- 
ready described.® It was not generally feasible to 
cut the spheroid exactly as proposed, with the 
(011) plane equatorial, and a method of de- 
scribing the actual orientation of the equatorial 
plane is therefore necessary. Let a be the angle 
between the rotation axis of the spheroid, which 
should be [011] and the nearest axis of form 
<100>, which must be either [010] or [001], 
and let 8 be the angle between the rotation axis 


*Q. E. Buckley and L. W. McKeehan, Phys. Rev. 26, 


261-273 (1925); L. W. McKeehan, P. P. Cioffi, Phys. Rev. 
28, 146-157 (1926). 

7S. L. Quimby, Phys. Rev. 39, 345-353 (1932). 

®R. G. Piety, Phys. Rev. 50, 1173-1177 (1936). 
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TABLE II. Work needed to magnetize. 
104 erg -«cm™=? 
Specimen ~~ Wire — Ww Win — Wie Ki K> 
K10 14 —0.18 —0.33 —0.72| — 2.4 
150 —0.08 —0.15 —0.32) — 1.2 
300 0 —0.03 0 — 08 
K30 14 0.17 0.17 0.68; — 1.5 
150 0.06 0.06 0.24) — 0.5 
300 0 0 0 0 
454 0 0 0 0 
K35 14 0.36 0.22 1.44) — 7.0 
150 0.32 0.25 1.28) — 4.8 
298 0.24 0.20 0.96| — 3.2 
490 0.12 0.12 0.48) — 1.1 
K50 14 0.83 0.43 3.32; —18.3 
150 0.69 0.52 2.76| —10.8 
300 0.46 0.36 1.84; — 6.8 


and the plane (100). In a perfectly cut spheroid 
a=45° and 8=0. It turns out that deviations of 
several degrees from the ideal values have a very 
small effect upon anisotropy in the equatorial 
plane because this passes through maxima and 
minima of the energy surface. The results can be 
corrected for such errors in cutting, as Shih? has 
pointed out in his latest paper, but this has not 
seemed necessary here. 

Table I gives information about the spheroids 
used in the present work. The demagnetizing 
factor N is that observed to be necessary to make 
the initial portion of J—/7 curves at room temper- 
ature coincide with the line 7=0. 


RESULTS 


The method of measurement and the reduction 
of observations to the form of J—H curves have 
already been described.* The results are presented 
in Figs. 1 to 4 and energy differences and energy 
coefficients are collected in Table II. For the 
method of computation and the meaning of 
symbols see the preceding paper by R. G. Piety.* 
The most unexpected result is that the direction 
for easiest magnetization changes not from 
<111> to <100> as iron increases beyond a 
critical amount, but from <111> to <110>. 
The transition occurs as nearly as these data can 
indicate at 24 percent iron, as is seen from Fig. 5. 
This is to be compared with Lichtenberger’s 29 
percent and Burger and Snoek’s 34 percent. (All 
compositions herein are in weight percent.) 

In the case of specimen K10, which is most 
nearly like nickel in its anisotropy, Akulov’s 
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method? of computing /—H curves, which neglects 
Kz, is reasonably successful. Assuming /—// along 
[111] as a basis, several points have been com- 
puted for J—// along [011] and [100]. These are 
plotted on Fig. 1 and fall near the curves actually 
found. 

One other peculiarity should be mentioned. As 
the temperature rises above room temperature 
the value of NV which is apparently best at room 
temperature does not make the /-// curve rise 
vertically from the origin. A smaller value of N is 
called for. The variation of N with temperature 
for a typical case is as follows: 


Specimen Temperature N 
35 14°C 0.761 
150 0.759 
298 0.754 
490 0.749 


This variation makes W, uncertain, so that its 
values have not been presented in Table II. No 
satisfactory explanation for this behavior pre- 
sents itself. It may be due to a variation in the 
thickness of a surface film magnetically very 
different from the interior. If so, the nature of 
this film and the factors upon which its thickness 


’N. S. Akulov, Zeits. f. Physik 67, 794-807; 69, 78-99 
1931). 


DECEMBER 15 1936 PHYSICAI 


Electrical Resistivity of Single Crystals 


Haroip E. Way, Physical Laboratory, 























Tene CRYSTALS 1181 
Or x10*erg «cm? 
BL 
6 
4- 
2 
0 
2e 
War Wie Wrem Wioo 
Ar 14°C é— aww 
150°C aio e-- 
i i L i i 
50 60 7T0%Ni 80 90 


Fic. 5. Variation of energy differences with composition. 


depends remain mysterious. A film of nonmag- 
netic material, everywhere of the same thickness, 
and thicker at higher temperatures, would change 
the apparent value of N in the observed manner. 

It is with great pleasure that the author takes 
this opportunity to thank Professor L. W. Me- 
Keehan, who suggested this problem, for his ad- 
vice and aid throughout the entire investigation. 
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of Some Dilute Solid Solutions in Zinc 


University of Iowa, Iowa City, Towa 


Received September 15, 1936) 


Sets of single crystals were grown of binary alloys con- 
sisting of dilute solid solutions of Cd, Cu, Ag, Au, Ni and 
Fe, each in zinc. Electrical resistivities at 20°C are ex- 
pressed in terms of principal resistivities, po and poo, these 
values coming from the usual cos? @ (@=orientation) plot 
of measured resistivities of a set of crystals. The resistance 
increases in a given series of alloys with increasing concen- 
tration, in a nonlinear fashion. The initial increase in 


INTRODUCTION 
HE resistivity and temperature coefficient 
of resistivity of dilute solid solutions of one 
metal in another have been studied frequently 
and the type of result to be expected is known at 
least qualitatively. Such studies in the past are 
apparently lacking for metals crystallizing in the 


nonregular system, probably due to the difficulty 


resistivity (above that of pure zinc) in micro ohm-cm pet 
atomic percent of solute is: Cd 0.94, Cu 0.4, Ag 1.1, Au 
2.3, Ni 51.0, Fe 300.0. This is correlated with the nearness 
of the solute metal to zinc in the periodic table. The ratio 
of principal resistivities is slightly higher than that for the 
zine crystal. Temperature coefficients decrease in such a 
way that Mathiessen’s rule is satisfied (with a maximum 


deviation of 4 percent) for all the alloys. 


of obtaining polycrystalline samples with grains 
oriented completely at random. It seemed worth 
while, therefore, to make a somewhat systematic 
study of single crystal samples, not only for 
the above reason, but because of the more funda- 
mental nature of the approach using single 
crystals. Zinc was chosen as the solvent material. 
Binary alloys were made, each of zinc and some 
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Designation of alloys and concentration of the 
solute material. 


TABLE I. 


Designation 





of alloys Percent by weight Atomic percent 

Cd 0.125 | 0.0721 

Cur, 2,3 125 0.250 0.500 | .1285 0.2571 0.5142 
Agi, 2, 3 .125 .250 500 | .07575 .1515 .3030 
Aut, 2,3 -125 .250 500 | .04144 .08289 -1657 
Fe1, 2 .005 .010 00585 01171 

Nit, 2 .005 .010 | .00557 1114 

FeNi .005 (Fe) +.005 (Ni) .00585 (Fe) +.00557 (Ni) 





one other metal in various concentrations. The 
resistivity and the temperature coefficient of 
resistivity from 0° to 40°C were measured for 
sets of single crystals of the unalloyed zinc and 
the various alloys. 


PREPARATION OF CRYSTALS 


The zinc used was “Bunker Hill,” a com- 
mercial zinc of exceedingly high purity, obtained 
from Platt Brothers and Company in the form 
of a single 50 lb. slab. By spectroscopic com- 
parison with a test sample! of known purity it 
was estimated that the B. H. zinc contained: 
Cd, appreciably less than 0.0008 percent; Pb, 
slightly less than 0.0047 percent; Cu, slightly 
less than 0.002 percent; Ni, Au and Ag, none 
detected; Fe, appreciably more than 0.0004 
percent (estimated by resistivity data, as shown 
later, to be about 0.002 percent), altogether 
certainly less than 0.01 percent impurity. 

Binary alloys were made of zinc and certain 
selected concentrations of each of the following: 
Copper, silver, gold, iron, nickel and cadmium. 
See Table I. The method of preparation was 
usually to add to a known weight of the molten 
zinc the appropriate weight of the second con- 
stituent. Occasionally the most concentrated 
alloy (of a series) was first prepared and the 
lesser concentrations made by dilutions from it. 
Where both methods were used to prepare a 
certain concentration a good check on their 
equality was obtained in terms of the resistivity 
measurements. A single ternary alloy (Fe Ni in 
Table I) was prepared by mixing in equal pro- 
portions by weight the 0.01 percent Fe and 
0.01 percent Ni alloys, thus giving an alloy of 
0.005 percent by weight of each. 

For the pure zinc and for each alloy a set of 
single crystals (from four to six) was grown with 





1A. W. Hanson, Phys. Rev. 45, 324 (1934), Table I. 





E. WAY 


orientations evenly distributed within the pos- 
sible range of orientation,? 0°-90°. The crystals 
were prepared as described by Cinnamon’ and 
in the same furnace. They were about 10 
centimeters long and with trapezoidal cross- 
sectional area of about 0.65 cm?. For the pure 
zinc Cinnamon’s ‘“‘growth conditions’’ were sub- 
stantiated, but for each alloy it was found 
necessary to determine appropriate growth con- 
ditions, which were toward higher gradients for 
about the same rates of growth used by Cin- 
namon. Keeping the ratio of gradient to rate of 
growth constant solely by decreasing the rate 
of growth did not produce single crystals. More- 
over, it was found that for the alloys of lowest 
concentration single crystals were grown more 
easily than for pure zinc, both because the width 
of the range of successful growth was increased 
and also because ‘‘optical mosaics,’"* were easily 
prevented. This was particularly noticeable with 
the smallest concentration (0.005 percent) of 
iron. With increasing concentrations growth 
conditions became more and more critical and 
it was not found possible to prepare a full set of 
crystals for greater concentrations than those 
given in Table I.° Also in most cases the highest 
concentrations were close to the probable limit 
of solid solubility. 


MEASUREMENT OF RESISTIVITY AND 
TEMPERATURE COEFFICIENT 

The resistance of a known length (about 8 cm) 

Th t fak length (about 8 
of the crystal was measured by comparison with a 
standard resistance using the customary poten- 
tiometer set-up. The crystal was placed in a 
wooden box and allowed to come to a constant 
temperature (room temperature, about 25°C) 
before resistivity measurements were made. The 
temperature was measured to 0.01 degree C by 
three thermocouples placed at each end and in 
the middle below the crystal. The distance 
between the potential contacts on the crystal was 
measured with a micrometer microscope and the 
cross-sectional area of the same part was meas- 

2 The angle between the length of the specimen and the 
vertical axis of the crystal. 

§C, A. Cinnamon, Rev. Sci. Inst. 5, 187 (1934). 

4H. K. Schilling, Physics 5, 1 (1934) and 6, 111 (1935); 
W. J. Poppy, Phys. Rev. 46, 815 (1934). 

5 With the exception of Cd, for which only one concen- 


tration was tried since a more thorough investigation would 
duplicate work known to be in progress elsewhere. 
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ured by the immersion method described by 
Hanson. Resistivity at 20°C was computed from 
these measurements and the temperature coef- 
ficient. 

The temperature coefficient was measured by 
placing the crystal in a large, sealed, brass tube 
which was immersed in water. The temperature 
was again measured by means of three spaced 
thermocouples, the brought out 
through a small brass tube joined on the larger 


leads being 


one and projecting above the surface of the 
bath. Resistivity measurements were made at 
three temperatures. For the first, ice was placed 
in the water surrounding the tube and stirred 
with a mechanical stirrer until the temperature 
of the crystal was constant as measured by the 
three thermocouples. Water at about 40°C was 
then substituted for the ice water and when 
equilibrium conditions were reached a second 
measurement was made. The third measurement 
was made after the crystal had come to, and 
remained for some time at, the temperature of 
the room. A resistance-temperature plot of these 
three points lay, within experimental error, along 
a straight line. The temperature coefficients of at 
least three crystals of each concentration were 
measured. The results for each concentration 
appeared constant within possible experimental 
error, there being indicated in particular no 
evident dependence on orientation. 


RESULTS AND DiscussION 
The resistivities of the Bunker Hill zinc and 
the zinc-copper alloys are shown in Fig. 1, each 


TABLE II. Principal resistivities in micro ohm-cm, their 


ratio and their difference, all at 20°C. 





Material po px po/ ps0 po — ps0 
B.H, Zn 6.253 5.902 1.059 0.351 

Cd; 6.345 5.970 1.063 375 
Cu 6.352 5.966 1.065 386 
Cur 6.406 6.023 1.064 383 
Cus; 6.466 6.086 1.062 380 
Ag: 6.348 5.986 1.060 362 
Ag: 6.425 6.062 1.060 | .363 

Ag; | 6.598 6.212 1.062 386 
Aw | 6.357 5.988 1.062 369 
Auz 6.436 6.072 1.060 .364 

Aus 6.608 6.219 1.063 389 
Fe, 6.450 6.070 1.063 380 
Fes | 6.590 6.203 | 1.062 387 

Ni; | 6.301 5.940 | 1.061 | 361 

Nie 6.336 5.954 1.064 382 
Fe Ni 6.490 6.111 1.062 379 


* A. W. Hanson, Rev. Sci. Inst. 7, 109 (1936). 
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Fic. 1. Resistivity at 20°C as function of square of 
cosine of orientation @ for B.H. zinc and series of alloys 
containing copper. Cu, =0.129, Cu,=0.257, Cu;=0.514 
atomic percent. 


plotted point representing one single-crystal 
specimen. Similar sets of curves were obtained 
for the other alloys. From these curves were 
read the values of the two principal resistivities, 
po and pg, parallel and perpendicular, respec- 
tively, to the vertical axis of the crystal. Table II 
includes these principals resistivities and other 
data for the B.H. zinc and the various alloys. 

It will be noted that the ratio po ‘py is appreci- 
ably increased (about 0.4 percent) by the smallest 
addition of any solute and is not changed very 
much by further additions. The behavior of the 
difference py—ps0 is not so regular. The first 
addition percent to 8 
percent) but there is a further increase in some 
cases as additional alloying material is added. 
The writer thinks it exceedingly probable that 
all the alloys are solid solutions with the solute 


causes an increase (3 


atoms replacing zinc atoms in an entirely random 
fashion. The initial changes in po/p99 and po— pso 
show that the crystal is slightly altered, sup- 
posedly in lattice dimensions, from the true zinc 
lattice by the addition of amounts of 
solute atoms. The disturbance due to further 


small 


additions is not so marked, as witness the con- 
stancy of po poo. The alloy designated Fe Ni was 
prepared to check the additiveness of the re- 
sistance change. Thus considering poo, the addi- 
tion of 0.005 percent Fe (Fe,) and 0.005 percent 
Ni (Ni,) cause an increase in resistance of 0.168 
and 0.038, respectively, summing to 0.206. For 
the Fe Ni the resistance change is 0.209 or the 
total effect of two additions is equal to the sum 
of the individual effects within experimental 


error. 
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Fic. 2. Principal resistivity poo at 20°C, as a function of 
atomic percent of solute. For dashed curves multiply 
abscissa scale by 210°. 


Although it seemed unlikely that any type of 
“ordered alloy’? would be produced with these 
small concentrations of solute, this possibility 
was tested. Since the attainment of an ordered 
condition would depend on the rate at which the 
alloys were cooled after solidification it might be 
that the cooling during growth of the crystals 
was too fast.’ A subsequent anneal and even 
slower cooling would test this point. Therefore, 
twenty-seven representative crystals had their 
resistivities remeasured several weeks after the 
first measurement.’ Twenty-four of these not 
only checked the first measurement but also 
showed no change in resistivity after two anneals, 
one at 200°C for 72 hours and the other at 
400°C for 36 hours. The cooling to room tem- 
perature after these anneals was extended over 


TABLE III. Relative effec ts of added elements. 


Initial Ap per Reciprocal* 


Material Group atomic percent effect 
Cd 2 0.94 
Zn 2 
Cu 1 4 0.29 
Ag l 1.1 a 
Au 1 2.3 9 
Ni &C 51.0 
Fe 8A 300.0 


* The figures in the last column are taken from Norbury's work 


and show the reciprocal effect, i.e., one atomic percent of zinc in copper 
produces an increase in resistivity of 0.29*10~-* ohm-cm et Nor- 
bury’s work was done with polycrystalline specimens 

7W. L. Bragg and E. J. Williams, Proc. Roy. Soc. 


A151, 540 (1935). 

5 Actually the crystals were left in the furnace for a 
period of fifteen hours while cooling from the melting 
point to room temperature. 

® First measurements were after the 


made directly 


growth of the crystals. The remeasurement was made, 
first as a check that no change had occurred due to a very 
slow relaxation at room temperature, and second to guard 
against possible changes due to accidental strains. 
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“1G. 3. Temperature coefficient of resistivity as function of 
atomic percent of solute. 


about 48 hours. To test the reverse point, i.e., 
whether an ordered condition had been attained 
during growth of the crystals, seven crystals 
were kept at 400°C for 24 hours and then 
quenched to 0°C by dropping them into iced 
water. No departure from the original resistivity 
occurred. It, therefore, seems likely that a 
simple random arrangement is stable both for 
slow and fast cooling. 

In Fig. 2 the resistivities are shown as a func- 
tion of atomic percent of the solute materials. 
Only one of the principal resistivities, pyo, is 
plotted but the other, po, would yield a similar 
plot. The relation between resistivity and con- 
centration is definitely not linear, the effect 
becoming less for equal additions as the con- 
centration increases. Nor do these curves fit a 
relation of the type proposed by Guertler :' 
R=R,+K(C—C?) in which Ro is the resistance 
of the pure solvent, K a constant, and C the 
concentration of the solute. The curves in fact 
bend the more 
rapidly than the above relation predicts. The 


towards concentration axis 
relative effects of the different added elements 
can be compared by computing (from the lowest 
the 


resistivity per atomic percent. This has been 


concentration in each case) increase in 
done and the results are arranged in Table III, 
the materials being listed in order of closeness to 
zinc in the periodic table, both horizontally and 
vertically. The data in the table could also be 
plotted to give a diagram similar to some given 
by Norbury" for various solutes in a single 
solvent metal. 

The relation between temperature coefficient 

10W. M. Guertler, J. Inst. Metals 6, 135 (1911). 


1! A. L. Norbury, Trans. Faraday Soc. 16, 570 
and J. Inst. Metals 33, 92 (1925). 


1921) 
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and atomic percent of solute is given in Fig. 3. 
The 
increased concentration and with the nature of 
the the 
resistivity previously discussed. The relation 


decrease in temperature coefficient with 


solute is very similar to increase of 
between these two effects is, however, brought 
out much better in Table IV, in which a test ts 
made of Mathiessen's rule. This rule states that 
QmPm = ap" in which a,, and p,, are temperature 
coefficient and resistivity, respectively, of an 
alloy, computed solely by “‘law of mixtures’’ and 
a and p are the actual measured quantities for 
the same alloy. For all the above described work, 
@» and p» are practically exactly the values for 
pure zinc, because of the use of such small con- 
centrations. Therefore, for Mathiessen’s rule to 
be obeyed the product ap for each alloy should 
be constant and equal to the value for the sol- 
vent material. The constancy of the last column 
in Table IV shows excellent agreement with this 
rule. The maximum variation from the value for 
the solvent (2206) is about 4 percent. 

The resistivity of the B.H. zinc was appreci- 
ably higher than the resistivity determined by 
Poppy‘ for two samples of zinc of approximately 
the same grade of purity. Moreover, Poppy's 


TABLE IV. Principal resistivity, pw, at O°C, temperature 
coe ficient and thetr product. 

Material py at 0 is a aps 

B.H. Zn 5.461 0.00404 2206 
Cd, 5.533 .00395 2185 
Cu 5.514 .00397 2189 
Cue 5.601 .00390 2184 
Cus 5.648 .00386 2180 
Agi 5.575 .00392 2185 
Age 5.674 .00385 2184 
Ags 5.768 .00381 2198 
Au 5.568 .00394 2194 
Aue 5.673 .00387 2195 
Aus 5.777 .00383 2213 
Fe, 5.630 .00391 2201 
Fe: 5.762 .00383 2207 
Nir 5.500 .00400 2200 
Nis 5.517 .00396 2185 
Fe Ni 5.676 .00383 2174 


2 See Hume-Rothery, The Metallic State, p. 51, but note 
that his a’s are temperature coefficients of conductivity. 
However, the similar relation for temperature coefficient 
of resistance is easily derived. It may also be shown that 
the relation a@ppm=ap follows by considering that the 
resistivity of the alloy is made up of two parts, one having 
the resistivity and temperature coefficient of pure zinc 
for these small concentrations) and the other being an 
additional resistance unaffected by temperature. 
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CABLE V. Comparison of various samples of zine 
Fe lines 
p po pu p detected 
E.W.B. (Poppy 6.161 | 5.842 | 1.055 none 
E.W.B. (Way 6.170 | 5.842 1.056, none 
E.W.R. (Poppy 6.218 | 5.882 | 1.057 18! 
B.H. (Way 6.253 | 5.902 | 1.058 37 


two samples differed somewhat, which he at- 
the 


It seemed worth while therefore to check up on 


tributed to a difference in iron content. 
this point. A set of crystals was grown of the 


material of lower resistivity used by Poppy 
(designated E.\W.B.) and their resistivities were 
determined. Spectrographic comparisons were 
made between this material, the B.H. zinc and 
the test sample” previously referred to. 

The results are summarized in Table V. It will 
be noted that the resistivities are qualitatively 
consistent with the difference in iron content. 
Moreover, the E.\W.B. is practically iron-free. 
In the light of the present work the other im- 
purities" in any of the zincs are believed to have 


Also, the 
=.W.B. sample agrees very nicely in poo with 


negligible effect on the resistivity. 


the determination of py) for spectroscopically 
pure zinc,” (09 =5.847 X10-* ohm-cm). Assum- 
ing thus that the values for E.W.B. are correct 
for 100 percent zinc, it is possible to extrapolate 
backwards from the iron resistivity-concentra- 
tion curve (Fig. 2) and determine the actual 
iron content of the B.H. zinc. The result is that 
the B.H. zinc contained 0.0019 percent Fe. 

Professor E. P. T. Tyndall's kind suggestions 
and cooperation in the development of this 
research problem are deeply appreciated by the 
author. 


18 None of Poppy's material of higher resistivity (desig- 
nated E.W.R.) was available. The test sample was said 
by Hanson (reference 1) to contain about the same amount 
of iron as the E.W.R. and it was used therefore in the 
spectrographic comparison. 

14 Lead should perhaps be excepted since it is an impurity 
in all the samples of Table V. However, a second sample 
of B.H. zinc, obtained after the completion of this work 
and differing from the first only in the fact that it had 
less lead (2 lead lines, as against 8 in first sample), had 
precisely the same resistivity. 


16E. P. T. Tyndall and A. G. Hoyem, Phys. Rev. 38, 


820 (1931). The value quoted above is corrected to a 
density of 7.13 g/cm', the authors having computed cross 
section of crystals by weighing a measured length and 
using a density of 7.15 g/cm’. 
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Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. The Board 
of Editors does not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Pressure Broadening of the HCN Band Lines and Inter- 
molecular Forces 


Some time ago, in an investigation of the fine structure 
of the photographic infrared HCN bands it was found by 
Herzberg and Spinks! that the strong line broadening 
observed is due to an exceptionally large pressure broaden- 
ing. It was further concluded that it is connected with a 
strong interaction of the HCN molecules at rather large 
distances and not due to ordinary Lorentz collision 
broadening. Recently however Cornell and Watson? rein- 
vestigated the HCN pressure broadening coming to rather 
different conclusions, namely that (1) the broadening is not 
particularly large but even slightly less than that found by 
Becker for the first harmonic band of HCI so that (2) dis- 
persion forces are the predominating broadening influence 
and there is no evidence for molecular association. Also the 
dependence of the line width on the rotational quantum 
number found qualitatively by Herzberg and Spinks was 
not confirmed by the measurements of Cornell and Watson. 

In order to investigate these differences one of us 
(W. W. W.) 
Herzberg and Spinks. The photograms were then measured 
by both of us. Also, to shed further light on the causes of 


has photometered the spectrograms of 


pressure broadening of the lines of infrared absorption 
bands of various molecules, we have photometered spectro- 
grams taken by Herzberg and Spinks of the 10,370A band 
of C,H,’ the 11,977A band of HCI‘ and the 11,300A band 
of CH,, all taken with atmospheric pressure of absorbing 
gas in a column 400 cm long. The slit width of the spectro- 
graph was 0.05 mm for all the plates and the dispersion 
varied between 5.11A and 5.35A per mm. The photometer 
traces have been made with the 40:1 magnification of the 
registering Koch-Goos microphotometer. 

Unfortunately, because these spectrograms were taken 
with the first Agfa plates sensitive to the 10,000A—12,000A 
region they have a rather heavy background fog, the lack 
of much contrast between absorption lines and continuous 
spectrum rendering the photometering difficult. Also there 
is no possibility of correcting the traces for any non- 
linearity of the emulsion characteristic. Line widths for 
consecutive lines in a band do fluctuate considerably in an 
irregular manner, therefore, and in particular one cannot 
be too certain about the measured relative widths of weaker 
as compared to stronger lines. Nevertheless, in view of the 
fact that all apparatus 
except for the pressure variation for the several HCN 


conditions remained constant 


spectrograms, the comparison of these roughly quantitative 
half-widths is of some significance. 


The average half-widths of all the P-branch lines meas 
ured in the 10,3844 HCN band are approximately 0.57, 
0.84, 0.94 and 1.37 cm™ at vapor pressures of 130, 265, 
360 and 650 mm, respectively. This is a somewhat more 
rapid rate of increase of line width with pressure than that 
found by Cornell and Watson. The explanation of this 
difference is undoubtedly the narrower apparatus line 
width in the case of Herzberg and Spinks’ spectrograms. 
There is evidence here for the nonadditivity of apparatus 
line width and pressure broadening width in the range of 
low gas pressures (under 1 atmosphere). Because of this 
the linearity of the dependence of line width and pressure 
is not as certain as was believed by Cornell and Watson. 
For an accurate investigation of pressure broadening the 
apparatus width should be small compared to the width to 
be measured, a condition which is not always easy to fulfill. 

The width of the HCN lines for atmospheric pressure, 
including apparatus width is about 1.5 cm™ from Herzberg 
and Spinks’ plates as compared to Cornell and Watson’s 
1.35 cm". 

From Herzberg and Spinks’ plates we obtain a half- 
width of about 0.5-0.6 cm™ for the absorption lines of CH,, 
C,H, and HCI at atmospheric pressure. Possibly the HCl 
line width is a trifle higher than that of the others. Con- 
sidering that most of this line width is apparatus width, it is 
seen that the HCN line width is in fact very considerably 
larger than those of CH,, C2:H»s and HCl in agreement with 
Herzberg and Spinks’ original qualitative estimates. Thus 
also Herzberg and Spinks’ original conclusion is confirmed 
that there are particularly large van der Waals forces in 
the case of HCN which perturb the absorbing molecules. 

As to the dependence of the broadening on the rotational 
quantum number the qualitative observation of Herzberg 
and Spinks of a decrease in line width for high J values 
could not definitely be confirmed by the quantitative 
measurements though they seem to indicate a slightly less 
broadening for both lines of high and low J value. Much 
more accurate measurements are necessary to decide this 


point. 
G. HERZBERG 
J. W. T. Sprxks 
W. W. Watson 


Department of Physics, 
University of Saskatchewan 
Sloane Physics Laboratory, 
Yale University, 
November 23, 1936. 


!G. Herzberg and J. W. T. Spinks, Proc. Roy. Soc. 147, 434 (1934 

2S. D. Cornell and W. W. Watson, Phys. Rev. 50, 279 (1936) 

3G. Herzberg and J. W. T. Spinks, Zeits. f. Physik 91, 386 (1934) 
4G. Herzberg and J. W. T. Spinks, Zeits. f. Physik 89, 474 (1934), 
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TO 


The Electroresistive Effect and a Rectifying Property of 
Carborundum Crystals 


The generalized Ohm’s law 


E;=ZRily, i=1, 2,3 1 


k=l 


does not adequately predict the actual current-voltage 
equilibrium states, particularly in the group of crystals 
which rectify alternating currents. Experiments performed 
by the writer show that the electrical resistances of single, 
homogeneous crystals of carborundum and zincite vary 
rapidly with the applied voltage even under isothermal 
conditions and when the voltage gradients are uniform. 
The isothermal variation of the electrical resistance with 
the applied voltage will be called the ‘‘electroresistive 
effect.”’ 

When the effects of the applied magnetic fields can be 
neglected, the current-voltage equilibrium may be approx 
imated by the power series 


Ea ZRakit So SRwlsatS 2S ORitial EE 
k= 1 lm~l 
te oe GRinncli 2. ++- Gat, 2.3 2 
k=1 l=1 ™ La 1 


? 


where /; and F; (?=1, 2, 3) are the cartesian components 


of the current and electric field intensity and Ry, Ri 
Rikim, etc., are components of second, third, fourth, etc. 
order tensors. These tensor components are in general 
functions of temperature and pressure. The terms con 
taining Rix, and Rixim are, respectively, linear and quad 
ratic in the components of E. They describe, respectively, 
the linear and quadratic electroresistive effects. The term 
current-voltage equilibrium is here used since in many 
cases several hours are required for the current to approach 
a steady value following a change in the voltage. 

The number of distinct or independent ‘‘electroresistive 
Rixt, Rikim, 


deduced from Neumann's hypothesis. Thus in the form of 


constants” etc., in any crystal class can be 
carborundum which belongs to the symmetry group Ce 
AE, the 
electroresistive constants for the linear and quadratic 


(Wyckoff’s notation) or (Voigt’s notation 


effects were found to be: 


Linear case 
R333; Rs ™ R322; 
Risi= Rese; Riis= Ree 
(3a) 


Quadratic case 
R3333; Risa = Ree22; Rsiis= 
Rsai = R322; Roa = — Roe = Rau; 
2Reu2= Rii— Re; Rive = Ren; 
Riai2= Rais; Risis = Roses; 
Riiss = Reo3s. 

(10 constants 3b) 





(4 constants) 





All other constants are zero. The integers 1, 2, 3 pertain, 
respectively, to the X, Y and Z axes. The Z axis is chosen 
as the axis of sixfold symmetry and the plane x =0 as one 
of the vertical planes of symmetry. 

Assume that in (2 


terms involving higher powers than 
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the second of the components of E are negligible and 
suppose that a uniform voltage gradient E; is established 
in a rectangular carborundum crystal such that £,; = E,=0. 
Since Ri, = Ra¥R 


of current along Z should obey the equation 


E3= R33] I;E 1;E;?. 4 


and Rj; =0 where «#7, the component 


+R +R 


The resistance of the specimen along the Z direction is 
given by 


R; . E; I; Ry T R333 Fs T R333 E ,, (5) 
oo ee ; 
if the sign of Es is reversed, the magnitude of R; changes 
by an amount depending upon R3;;E3, the term which 
originates from the linear electroresistive effect. Hence 


this theory affords an explanation of the rectifying property 
of carborundum crystals. The portion of rectification due 
to a point and plane is, of course, not taken into account 
by this phenomenological theory. Any term involving odd 
powers of E in (2) will affect the rectifying property. 

Carborundum crystals were found, as predicted, to 
rectify currents caused to flow along Z by applying uniform 
voltage gradients E;sin wf. Preliminary experiments 
indicate that in the case of voltage gradients which are 
greater than one volt/cm higher order tensors than the 
fourth are needed. 

A later paper will contain a more complete account of 
the theory together with a list of the linear and quadratic 
constants for the various crystal classes and experimental 
data upon the intricate electroresistive phenomena in 
large, single crystals of carborundum and zincite. 


HAROLD OSTERBERG 
University of Wisconsin, 
November 24, 1936. 


Exchange Forces Between Elementary Particles 


The Fermi theory of 8-decay introduces an interaction 
term into the Hamiltonian of the elementary particles: 
proton, neutron, electron and neutrino. A consequence of 
this theory is the existence of an exchange force between 
protons and neutrons; this force is caused by the emission 
of light particles by one heavy particle and the absorption 
of these by the other particle, both processes occurring in one 
elementary act. It has been pointed out that only an inter 
action term containing the derivatives of the wave func- 
tions can give rise to a force of sufficient magnitude between 
heavy particles. The interaction term has the form 


Q = g-acP*TWP*( x)? (x) - (0/Ox 


mye 


x) 


0/dx)"¢7(x)+conj. (1 


wv, *, ¥, ¢ are the wave functions of the proton, neutron, 
electron, and neutrino, respectively; x stands for the four 
space-time coordinates; the Greek letters are the spin 
indices; (0/dx)™ is a differential operator of degree m with 
respect to space and time; @®’’’ is a combination of Dirac 
matrices, making the term relativistically invariant; such a 
combination can be constructed in several ways. 

It can be shown that the interaction energy (1) also gives 
rise to exchange forces between other elementary particles. 
Let us consider a neutron at a place I and an electron at a 
place II. Let the neutron decay into a proton, an electron, 
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and a neutrino. In a subsequent process the proton can be 
transmuted back into a neutron at place II by absorbing 
the neutrino emitted and the initial electron, leaving the 
electron emitted in the first process at place I. The net 
result is that the initial neutron and electron have changed 
places. The two processes occurring in one elementary act 
produce an exchange force analogous to that between 
heavy particles. 

In general, exchange forces between two elementar\ 
particles A and B occur whenever B can be created by 
decay of A, or A by a decay of B. Such forces are possible 
between: proton and neutron; proton and positive electron; 
proton and neutrino; neutron and neutrino; neutron and 
negative electron; neutrino and electron (positive and 
negative 

To get the last two of these exchange forces we must 
assume the existence of a transmutation of an electron into 
a neutrino or vice versa with the simultaneous creation of a 
pair of heavy particles. It should be possible for a positive 
electron, for instance, to decay into a neutrino with emission 
of a proton and a neutron. It is easily shown that this kind 
of process is a consequence of the interaction term (1 
Since we assume that the wave functions of all elementary 
particles follow the Dirac equation, we are forced*to con 
sider all negative energy states as occupied in the empty 
space. Then every interaction energy effecting a transmu- 
tation causes also a creation of two particles, since a 
transmutation of a particle in a negative state into another 
particle in a positive state must be considered as a creation 
of two particles. 

A rough evaluation gives the following result: the poten 
tial of the exchange force between two elementary par 
ticles a and + in the distance r is of the following order of 
magnitude: 


s 


Fre ~er(= ‘\’ a(n oS (2 ) hy. 
€ € 


€ 





Here « is the mean energy of @-ray particles and go is a 
dimensionless constant go=g-(e/hc)3*™*"e'~5 X10-" the 
order of magnitude of which is independent of the values 
assumed for m and n. S, and S, are the degrees of differen- 
tiation of the wave function of a and } in (1); p, and pf, are 
the corresponding momenta. We assume as usual that this 
expression only holds for distances r >d, d being of nuclear 
dimensions, and that J®(r)~J#(d) for r<d. Now m=1, 
n=2 is the only combination which can account for the 
observed proton-neutron force and for the asymmetrical 
position of the maximum in the intensity curve of 8-ray 
spectra. We obtain for the exchange forces at zero distances 


(in volts) between 


proton-neutron: ~10* neutron-neg. electron: ~0.5 < 10# 
proton-electron: ~0.5 X 104 neutron-neutrino: ~2 X 107 


proton-neutrino: ~2 X 10~! electron-neutrino: ~10 


(The energies of the light particles are taken as about 10° 
volts.) Thus the exchange forces, except in the case of two 
heavy particles, are too weak to be observed. 

V. F. Wetsskopt 


Institute for Theoretical Physics 
Copenhagen, Denmark, 
November 5, 1936 


¥=J—xdJ/dx to 


mentary functions” 
Kulenkampff-function) on which the analysis formerly was 
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On the Hard Component of Cosmic Rays 


In a former paper' it was shown that the depth-ionization 
curve of cosmic rays measured 80-240 m below sea level 
can easily be represented by an exponential function. The 
newer measurements of Weischedel? confirm this in prin 
ciple. Therefore it is of interest to examine what range 
distribution is obtained by extrapolation of this exponential 


behavior of the curve due to rays coming in from all direc 


J= Joe aad l 


The intensity curve for rays descending vertically through 


converted by the transformation 


For the average range one obtains 


and for the most frequent range 


Y=Joe 4r(1+ ux 2 
and the range distribution p= —dy/dx will be finally 
p(x dx = p?x e # dy 3 
R=2/u } 
5 


Res =I1/u 


1 gives these relations with the value «=1.889x 10 
per meter of water, formerly obtained! for the hard com 
ponent. A range phenomenon, curve II, is quite apparent 
after transformation to vertical incidence. The range dis- 
tribution curve III is very like a Maxwellian one. A direct 
conclusion concerning the energy distribution is not yet 
available. In particular one cannot suppose a loss of energy 
exclusively produced by primary ionization, for this would 
require such small energy values of the incident rays that 
in the equatorial belt only a fraction could reach the earth's 
surface through the blocking effect of the earth's magnetic 
field. And then, contrary to all observations there should 
also result a strong latitude effect of this component 

It was pointed out by Lenz,‘ that the different “‘ele 


exponential function, Gold’s integral, 











meters of water 
Observed intensity of the hard component extrapolated 
II) Intensity after transformation to vertical incidence 
(III) Range distribution 
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based, show a maximum after being multiplied with x. 
Thereupon he bases a method of analysis, forming the 


‘deformed curve’’ x-J(x) by multiplying the measured 


values J by the absorbing depth. In the case of a band 


structure of the radiation this curve includes various 


maxima and so may be decomposed into deformed ele 
the relations above mentioned a 


mentary functions. By 


totally new smgnificance is given to this proceeding. This is 
seen at once by writing Eq. (3) in the form 
p(x) =p*xJ(x)/Jo. 6 


This is exactly the deformed intensity curve. In the case of 
an exponential depth-ionization curve (or one representable 
by a sum of exponential functions) it is seen that the method 
of Lenz‘ signifies physically a decomposition into different 
In the depth-ionization curve—the 


range part 


above 80 m 


groups. 
may be considered to be formed by two ex 
ponential functions. Thus the application of the Lenz- 
method means a decomposition into two range groups, the 
harder one we have treated above more rigorousl\ 

With great pleasure | acknowledge my indebtedness to 
Professor Dr. E. L.. da Fonseca Costa for his helpful 
advice. 

B. Gross 


Instituto Nacional de Technologia 
do Ministerio do Trabalho, 
Rio de Janeiro, 
October 28, 1936 

! B. Gross, Zeits. f. Physik 92, 755 (1934) 

2 F. Weischedel, Zeits. f. Physik 101, 732 (1936). 
} B. Gross, Zeits. f. Physik 83, 217 (1933). 
4‘ E. Lenz, Zeits. f. Physik 83, 193 (1933). 


Direct Effects of Particular Solar Eruptions on Terrestrial 
Phenomena 

In this journal previously, I' called attention to the 

occurrence of widespread sudden fadeouts of high frequency 

radio transmission, and suggested that study be undertaken 

to determine possible relations of the effect to solar, earth- 

I his 


led to effective cooperation in such study, which has been 


current, terrestrial magnetic, and other phenomena 


rewarded by the observation of simultaneous happenings 
on numerous occasions in each of the four fields named. 
The work has corroborated the hypothesis advanced in the 
1935 article that a radio fadeout is due to a sudden solar 
emanation lasting only a few minutes. 

In the two years from July 1934 to June 1936 there were 
94 reported radio fadeouts; 39 of these are known to have 
been accompanied by a simultaneous visible eruption on the 
sun; 15 of them are known to have been accompanied by 
simultaneous fluctuations of terrestrial magnetic elements; 
and for 4 of fluctuations of earth 
currents have been reported. The fluctuations in terrestrial 


them simultaneous 
magnetism and earth currents are abrupt pulses of about 
the same duration as the radio fadeouts, not similar to the 
fluctuations characteristic of magnetic storms. In none of 
these fields has there been an opportunity to make a 
comprehensive study of world-wide data, and there is good 
reason to suppose that the number of coincidences has been 
greater than is known at this time. 

The occurrence of visible effects in the causative solar 
eruptions is probably more or less fortuitous. The cause of 
the terrestrial phenomena is a sudden outburst of highly 
ultra 


radiation of ultraviolet (i.e., 


penetrating solar 
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visible) frequencies. This produces sudden ionization of the 


atmosphere below (or possibly in) the E-laver of the 
ionosphere, where the air density is great enough to insure 
numerous collisions of moving ions and hence rapid 
absorption of the radio wave energy during the few minutes 
it lasts. The degree of absorption is less for the higher 
radiofrequencies than for lower ones. The ionization density 
of the F:-layer is not affected. With the sudden ionization 
below (or in) the E-layer, sudden changes occur in the 
currents flowing there, currents which are the cause of a 


certain part of the earth’s magnetism. The changes in 
these currents directly cause the special magnetic fluctua- 
tions which are observed, and also inductively produce 
changes in earth currents. 

The absorbing ionization and the radio fadeouts are 
produced on the illuminated side of the earth and do not 
They 


occur on the dark side are most intense where the 


sun’s radiation is perpendicular (i.e., around noon and at 


low latitudes). This distribution is also true of the terrestrial 
magnetic variations, and probably of the earth current and 
other phenomena if any. The effect is thus quite distinct 
from ordinary magnetically disturbed periods, sometimes 
storms, whose effects are much more 


called magnetic 


marked at high than low latitudes. Ordinary magnetic 
disturbances and the special effect here reported have no 
necessary connection with each other, except for a general 
tendency to occur more frequently when sunspots are more 
numerous. Ordinary magnetic disturbances occur simul 
taneously with lowered ionization density in the F»-layer, 
as shown by my colleagues. Thus magnetic storms and the 
special magnetic effects simultaneous with radio fadeouts 
are associated with different levels in the atmosphere. It is 
believed that certain rare observations reported in the past 
by astronomers, of magnetic pulses simultaneous with 
visible solar effects, such as the coincidence reported by Pro- 
fessor C. A. Young as occurring on August 3 and 5, 1872, 
were the special type of magnetic effect now found to be 
associated with radio fadeouts. 

The radio fadeouts reached an extraordinary climax in 
May and June, 1936. In the three months, April to June, 
there were 54 reported fadeouts, more than the total 
number reported previously (in the two years the phe 
nomenon has been known). Then in the following three 
months, July to September, there were only 16 fadeouts 
Of the 54 fadeouts in April to June, 21 were accompanied 
by solar eruptions, and 10 were accompanied by the special 
magnetic effect. As there is reason to suppose that radio 
fadeouts (and the associated magnetic and other phe 
are caused by one type of ultraviolet (i.e., ultra 
that 


disturbances, auroras, etc., are caused by a solar emanation 


nomena 


visible) solar radiation, and ordinary magnetic 
of different character, these phenomena present powerful 
means of keeping track of two different and very interesting 
types of events on and in the sun. They deal with ema 
nations which are not only invisible but probably do not 
reach the earth’s surface at all. 

J. H. DELLINGER 


National Bureau of Standards, 
November 24, 1936 


! Dellinger, Phys. Rev. 48, 705 (1935 
? Kirby, Gilliland, Judson, Smith, Reymer, Phys. Rev. 48, 849 (1935) 
and $0, 258 (1936). 
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Investigations of the Deuteron-Deuteron Reaction 


Last year we reported some results on the yield of 
neutrons produced by bombarding different targets con- 
taining deuterium, beryllium or lithium with ion beams of 
deuterium, which were accelerated by a 400 kv transformer 
rectifier outfit.! 

In continuation of these investigations we have studied 
in more detail the neutrons and the protons produced by 
bombarding different deuterium targets with fast deuterons 
(‘‘D-D” reaction). This case is especially interesting on 
account of the apparently very high yield? (about 10~*) and 
because these neutrons have all the same energy”: ? (about 
2X 10° volts). 

The accelerating voltage is measured both with an 
improved generating voltmeter and with a high resistance 
and microammeter.‘ The current carried by the mag- 
netically analyzed ions is recorded by a current integrator. 
In some experiments the molecular ions (mass 4) from a low 
voltage arc are used. Voltages applied between the Faraday 
cage containing the target and adjacent diaphragms and 
an auxiliary small magnetic field prevent any effects from 
secondary electrons. The protons produced in the D-D 
reaction penetrate thin aluminum windows and are counted 
by means of an ionization chamber, linear amplifier and a 
scale-of-eight thyratron counter.’ The neutrons are re- 
corded with the same apparatus, by the use of a paraftin- 
lined chamber for the fast neutrons, a boron-lined chamber 
for slow neutrons, or a simple air chamber (for calculating 
absolute values). The yield of the neutrons was compared 
with that of neutrons from a 15 mg radium-beryllium 
source replacing the target. A comparison of the radium 
source with a radon beryllium source of the same gamma- 
ray activity showed the neutron activity to be the same 
within 10 percent.® 

As solid targets of D-compounds did not last long under 
bombardment, we used liquid D;PO, for longer series of 
measurements, which proved to be very satisfactory. 

The results obtained so far are the following: 

1. The excitation function of the ‘‘thick” targets be 
tween 40 and 100 kv is the same for protons and neutrons 
as found before by Alexopoulos? and agrees approximately 
with that given by Oliphant, Harteck and Rutherford. The 
relative yield is 1 : 4.4: 8 for 100, 200 and 300 kv. No 
resonance has been found between 40 and 100 kv although 
measurements were made every 5 kv. 

2. The increase of slow neutrons absorbable by 0.6 mm 
of cadmium with increasing thicknesses of paraffin cylinders 
surrounding the source and the ionization chamber has 
been measured, both for the D-D neutrons and the Ra-Be 
neutrons. The two such curves coincide above 5 cm of 
paraffin. For smaller thicknesses relatively more slow 


neutrons are observed from the Ra-Be. This effect is 


probably due to the large number of relatively slow 
neutrons originating from the Ra-Be source. 

3. The yield of Cd neutrons produced in 7.5 cm of 
paraffin from a D;PO, target bombarded by 1uA D-ions at 
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100 kv is the same as that of our 15 mC Ra-Be source; but a 
slight change of the geometry of the neutron source inside 
the paraffin has an appreciable effect. The yield of fast 
neutrons for luA D-ions at 100 kv is the same as that of 
44 mC Ra-Be. The main reason for this difference is the 
higher initial proportion of relatively slow neutrons from 
Ra-Be which give no measurable recoil protons in the 
paraffin-lined chamber. 

4. The total number of protons produced in D;PO,, 
measured at 90° from the direction of the ion beam, is 
7.10* sec. at 100 kv and izA on the assumption of iso- 
tropic emission in all directions. Actually appreciably 
more neutrons were found in the direction of the beam than 
at 90°, in accordance with recent reports from the Caven- 
dish laboratory.’ Multiplying the yield by 50/3, the ratio 
of the number of electrons in D;PO,; and the number of 
D-atoms, one gets as yield for the D-D reaction 1 proton 
for 6.10° impinging deuterons. Solid targets—(ND,)2SO,, 
ND,Cl, KOD 
correcting for the stopping power in the same way as fo 
the D;PO, target. Oliphant, Harteck Rutherford 
estimate the yield to be of the order of magnitude of 1 in 
10 at 100 kv. 

5. When using the simple air chamber for counting fast 
neutrons and taking the cross section of 1.810~*4 cm? 
given by Dunning’ for the scattering of neutrons from 
Rn-Be, we find about half as many neutrons as we found 


under the assumption that every recoil 


give within 20 percent the same yield afte: 


and 


protons before, 
O or N atom from scattered neutrons produces a count in 
our ionization chamber.’ So it is safe to say that the two 


reactions produced by the D-D collision give approxi 


mately the same yield. 
We wish to express our appreciation to the Rockefeller 


Foundation for a special grant and to Messrs. Y. Beers, 


J. Giarratana and M. H. Kanner for valuable help in some 
of the experiments. 
R. LADENBURG 
R. B. ROBERTS 


Palmer Physical Laboratory, 
Princeton, New Jersey, 
November 12, 1936. 


! Ladenburg, Roberts and Sampson, Phys. Rev. 48, 467 (1935) 

2 See Oliphant, Harteck and Lord Rutherford, Proc. Roy. Soc. 144, 
612 (1934); Alexopoulos, Helv. Acta 8, 513 and 601 (1935). 

3 Dee, Proc. Roy. Soc. A148, 623 (1935); Bonner and Brubaker, Phys 
Rev. 49, 19 (1936). 

‘ Hafstad, Heydenburg and Tuve, Phys. Rev. 50, 504 (1936) 

5 This device was constructed for these experiments by J. Giarratana 
who developed a special Wynn-Williams circuit, adapted for RCA 
thyratrons. 

* Our method for mixing Ra and Be is described by Erbacher and 
Philipp, Zeits. f. physik. Chemie A176, 175 (1936) who used a similar 
method. 

7 Dunning, Phys. Rev. 45, 586 (1934). 

§ Alexopoulos estimates (reference 2, p. 513) the number of neutrons 
produced by lwA and 100 kv to be 4.10 per second, i.e., about 50 times 
as big as we find, and without correcting for the stopping power of his 
target he gives as yield 1 neutron per 1.6X10® deuterons of 100 kv 
The yield given by Kikuchi, Aoki and Husimi (Proc. Phys. Math. Soc 
Japan 18, 122, 299 (1936)) does approximately agree with our value, 
but the percentage of D-ions in their beam is only estimated and not 
measured by magnetic analysis. 

* Kempton, Brown and Maasdorp, Proc. Roy. Soc. 157, 372 (1936). 
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Erratum: On the Binding of Neutrons and Protons 


J. B. Fisk, Soctety of Fellows, Harvard University, 
L. I. Scurer, Massachusetts Institute of Technology, 
J go 

W. SHockLey, Bell Telephone Laboratories 


! Phys. Rev. 50, 1090 (1936 


In the the manuscript for the above 


Letter to the Editor the wrong drawing was used in making 


preparation ol 
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Fic. 1. The variational binding energy of Ha vs. the “effective radius 
of interaction,”’ ro, for several values of D’/D, the ratio of depths of 
neutron-neutron to neutron-proton holes. 


the cut for Fig. 1. The authors are not responsible for this 
error and the correct cut for their figure is shown here. 


AMERICAN INSTITUTE OF PHysiIcs 
Publications Department 


Erratum: Precision Cosmic-Ray Measurements up to 
Within a Percent or Two of the Top of the Atmosphere 


R. A. MILtrKAN, H,. V. NEHER ANDS. K. Haynes, Norman 
Bridge Laboratory of Physics, Pasadena, California 


Phys. Rev. 50, 992 (1936 


To correct an error in the checking of proofs in this 
office the word “proton,” last line page 995 and line 23 
second column, page 998 should read “photon.” 


AMERICAN INSTITUTE OF PHYsSICs 
Publications Department 


The 8-Ray Spectrum of Li*® 


According to a recent communication of Rumbaugh and 
Hafstad' there is a discrepancy between the mass of Li’ as 
obtained from the reaction Li?+H*—Li’+H! and the re- 
action Li*~+Be*+e~. From the former they conclude that 
Li§>8.022(8) and from the latter? Li’=8.019(3)+0.001. 


This discrepancy can be explained by supposing that the 
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8-ray spectrum of Li® gives rise to an excited Be* nucleus 
so that the reaction reads Li*—>(Be*)’+e~ with the under- 
Be’)’, respectively, for excited 
and normal nuclei, having an energy difference (Be*)’ — Be® 

>0.0035 +0.001 mass units. In order that this explanation 
should work it is necessary to suppose that the disintegra 


standing that Be® stand, 


tion of Li’ into Be*+e™ is ruled out by a selection rule and 
that there is available an excited level of Be* at least 0.002 
mass units above the normal which is not ruled out. Ac 
cording to the calculations of Feenberg and Wigner® there 
exists a 'D level of Be® at approximately 0.003 mass units 
above the normal !S level. The (Be*®)’ nucleus is in the 'D 
level according to this explanation. On the other hand 
their calculations for Li® give a P level which is built on 
doublets in neutrons and on doublets in protons. The ob 
served effect can, therefore, be explained by supposing that 
the normal level of Li’ is a *P: because then‘ the selection 
rule At=0, +1 for total angular momentum will make the 
transition to (Be*)’ possible and will rule out that to Be’. 
An explanation could be attempted by supposing that Li‘ 
is left in an excited state (Li’)’ after the ejection of the pro 
ton. It could then emit a y-ray, become a normal Li*® and 
finally emit a 8-ray. The energy of (Li® 
about 4 or 3 Mev above that of Li’. It should be more 
probable that the protons will be emitted with the full 


’ would have to be 


available energy and that the Li’ nucleus will be left in its 
normal state. One could suppose that Li’ cannot be formed 
on account of a selection rule. The likely selection rules are 
those of the orbital, spin and total angular momenta as well 
as of parity.’ The normal state of Li’ is* an odd *P 5/2 state, 
that of H? an even 4S; and that of Li® an even *P. It is 
possible to have conservation of parity® and of the total 
orbital and spin angular momenta as well as of their re 
sultant by assuming the angular momentum of relative 
one 


motion to be (a) zero initially and one finally or else (b 


initially and zero finally. Case (a) is more probable than 
(b) because the incident deuterons have a lower energy 
than the protons which could be emitted with the forma 
tion of Li’. The fact that the final angular momentum of the 
proton is one and not zero does not decrease the probability 
of the reaction by a large factor because the addition of 
h*? / Mn* to the potential barrier amounts to about 5 Mev at 
a distance of 4 10~" cm. There is, thus, no special reason 
why the formation of Li’ in its normal state should not 
occur. It appears, therefore, that the formation of an ex 
cited Be® is a more natural explanation. Gamma-rays of 
about 3 Mev energy should accompany the reaction on the 
present view. 

E. WIGNER 

Princeton University, 
G. Breit 


University of Wisconsin, 
November 22, 1936. 
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Thermal conductivity of refractories, J. L. Finck 
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MacPherson, M. A. Starr—389(A 
Mobility of ions (see Ions, mobility 
Molecular radii (see Collision diameter 
Molecular structure and constants (see 
molecular and Raman spectra) 
Of S. by electron diffraction, L. R. 
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Relativity 


Conservation of photons, H. P. Robertson—385(A 


Michelson- Morley experiment, W. B. Cartmel—1097(.A 

New relativity, D. G. Bourgin—864 

And nuclear forces, H. Margenau—342 i 
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1186(L) 

Of He, G. H. Dieke—797 

Of OH*, calculation of constants, M. Kovner—188(I 

Of InO, W. W. Watson, A. Shambon—607 

Low electronic states of heteropolar diatomic molecules, 
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X-rays, absorption 

Of compounds of Br, Ru, and Sr, S. T. Stephenson 
387(A), 790 

By elements of high atomic number, M. T. Jones—110 

By elements of low atomic number, T. R. Cuykendall 
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